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ABSTRACT 

There  is  a  growing  consensus  that  die  sound  generated  by  breaking  waves  is 
responsible  for  much  of  fee  ambient  noise  level  in  the  ocean.  While  numerous  field 
measurements  have  shown  a  strong  correlation  between  the  ambient  noise  spectrum 
level  (N)  in  the  range  100  Hz  to  25  kHz  and  wind  speed  in  die  ocean,  very  little  has 
been  done  to  establish  a  comparable  correlation  between  the  ambient  noise  spectrum 
level  and  surface  wave  field  parameters.  The  difficulty  in  establishing  this  relationship 
is  remarkable  given  that  the  frequency  and  intensity  of  wave  breaking  are  dependent 
on  the  characteristics  of  die  wave  field. 

In  Fall  1991,  an  experiment  was  conducted  from  the  research  platform  Flip  130 
kilometers  off  the  coast  of  Oregon,  where  the  ambient  noise  between  2.5  and  25  kHz, 
the  wind  speed,  and  the  sea  surface  elevation  using  wir  wave  gauges  weremfeajainicf.  f**? 

The  correlation  between  N  and  the  root  mean  square  wave  ampfitudedStfas  * 
found  to  be  poor  but  could  be  unproved  if  the  swell  was  filtered  out  fsim^ewfave  ,  ' , 
elevation  time  aeries.  The  influence  of  swell  on  the  value  of  a  was  disprop^rtionate  to  .  , 
the  level  of  ambient  noise  since  its  characteristics  were  not  directly  dwe  to  the  local 
wind-wave  conditions.  Observations  of  the  dependence  of  the  high  ffequepcy  wind 
waves  and  die  directional  wave  spectrum  under  turning  winds  suggested  that  the  high 
frequency  wave  components  responded  more  quickly  to  changes  in  the  wind  speed  and 
wind  direction  than  die  energy-containing  frequencies. 

The  ambient  noise  level  also  correlated  well  with  the  root  mean  square  wave 
slope  $.  This  is  consistent  with  previous  laboratory  measurements  which  showed  that 
die  steepness  of  a  packet  of  waves  correlates  with  the  strength  of  wave  breaking  and 
with  characteristics  of  breaking  waves  mifh  as  loan  of  momentum  fliiv  diggipatinnr 
initial  volume  of  air  entrained,  mixing,  and  sound  generation.  5 

Comparisons  of  surface  wave  dissipation  estimates  using" fiel4. measurements 
and  models  developed  by  Phillips  (1985)  and  Hassehnann  (1974)  show  that  although 
die  two  models  have  very  different  forms,  they  give  values  that  are  comparable  in 
magnitude  Die  relationship  between  die  ambient  noise  level  and  log  of  dissipation  give 
correlation  coefficients  (0.93-0.95)  that  are  comparable  to  those  between  ambient  noise 
and  wind  speed.  Die  mean  square  acoustic  pressure  was  shown  to  vary  with  the 
dissipation,  with  p2  «  D0-^0-8.  The  results  suggest  that  measurements  of  ambient  sound 
may  prove  to  be  useful  in  inferring  surface  wave  dissipation. 

Thesis  Supervisor  W.  Kendall  Melville 
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experiment  (43°  42'  N,  125°  59'  W).  The  site  is  approximately  130  km 
west  of  Oregon  and  600  km  south  of  Vancouver  Island. 

Figure  3.2.  A  schematic  diagram  of  Flip  showing  the  locations  of  the  69 

different  instruments  used  in  the  experiment.  The  capacitance  wire 
wave  gauge  is  15  meters  and  the  resistance  wire  wave  gauge  is 
17  meters  from  the  hull  of  Flip. 

Figure  3.3.  Schematic  diagram  of  the  sampling  and  preprocessing  70 

method  used  in  gathering  and  storing  data. 

Figure  3.4.  A  photograph  of  the  directional  hydrophone  taken  72 

immediately  after  the  platform  'flipped'  from  vertical  to  horizontal 

attitude.  The  photograph  was  taken  from  a  location  near  the  bow 

of  Flip.  The  hydrophone  is  the  dark  circular  disk  at  the  center  of  the 

photograph. 


Figure  3.4.  A  dose-up  photograph  of  the  directional  hydrophone  taken  73 

while  Flip  was  in  port  in  San  Diego.  The  instruments  in  the 
foreground  are  the  SIO  Doppler  sonars  used  to  image  the  sea  surface. 

Figure  3.6.  A  simplified  diagram  of  the  directional  hydrophone  and  the  74 

mounting.  The  dimensions  are  not  to  scale. 

Figure  3.7.  Sensitivity  level  (SL)  of  the  directional  hydrophone  in  78 

dB  re  1  VoltVlpPa2  obtained  from  a  calibration  of  toe  instrument 
at  toe  Naval  Ocean  Systems  Center. 

Figure  3.8  Plot  of  array  gain  (AG)  of  the  directional  hydrophone  as  a  79 

function  of  frequency.  AG  values  computed  from  NOSC  calibration  are 
denoted  by  O. 

Figure  3.9.  Some  typical  ambient  noise  power  spectra  in  dB  re  lpPa2/  81 

1  Hz  at  various  wind  speeds. 

Figure  3.10.  Photograph  of  the  4-wire  resistance  wave  gauge  array  83 
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and  the  capacitance  wire  wave  gauge  (heavy  black  cable)  deployed  at 
the  end  of  the  boom.  The  heavy  white  cable  near  the  center  of  the 
photograph  ia  die  ground  wire. 

Figure  3.11.  Photograph  of  the  bottom  spreader  bar  being  lowered  to  84 

position  prior  to  full  deployment  of  die  wire  wave  gauge  array.  The  bar 
la  made  of  stainless  steel  and  each  arm  is  0.5  meters  long. 

Figure  3.12.  Plot  showing  the  values  for  the  four  different  field  85 

calibrations  of  die  capacitance  wire  wave  gauge  The  dashed  line  is  a 
linear  least  square  fit  of  die  sample  points.  The  equation  for  the  best 
linear  least  squares  fit  is 

a^m)  *  0.541  vc  (Volts)  +  0.800 

where  is  the  calibration  position  and  vc  is  the  mean  voltage  reading  at 

that  position. 

Figure  3.13.  Plots  showing  file  characteristics  of  a  typical  wave  height  89 

power  spectrum.  The  data  were  sampled  at  8  Hz.  The  spectra 

were  computed  from  a  one  hour  wave  gauge  record  with  U10  *  8  m/s. 

a)  Power  spectrum  computed  from  1024-point  FFTs  used 

throughout  this  work,  b)  Power  spectrum  over  the  same  period 

computed  from  2560-point  FFTs  to  resolve  the  lower  frequency 

spectral  peaks. 

Figure  3.14.  Plan  view  schematic  diagram  showing  die  location  and  91 

orientation  of  the  resistance  wave  gauge  array  and  the  coordinate 
systems  used  in  computing  the  wave  slope  s. 

Figure  3.15.  Plan  view  of  schematic  diagram  showing  the  coordinate  96 

system  used  in  computing  die  directional  wave  spectrum  $((0,6). 

Figure  4.1.  Time  series  of  Ujq,  u*  and  the  wind  direction  during  the  101 

NOBS  experiment  Note  that  die  wind  rapidly  changes  direction 
on  Julian  day  278  after  blowing  from  die  North  for  approximately 
six  days. 

Figure  42.  Some  typical  512-point  wave  height  spectra  $(f)  [m2/Hz]  103 

taken  from  one  hour  records  at  the  start  of  each  day  of  the  NOBS 
experiment  a)  269.  b)  270.  c)  271.  d)  272.  e)  273.  f)  274.  g)  275.  h)  276. 
i)  277.  j)  278.  k)  279. 1)280. 

Figure  43.  a)  Root  mean  square  amplitude  s  (m)  and  significant  105 

wave  height  H,  (m)  computed  from  the  full  spectrum,  b)  Mean 
frequency  computed  from  the  full  spectrum  f .  c)  a*  and  H, 
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computed  with  swell  filtered  out  d)  Mean  frequency  with  the  swell 
filtered  out  f^,. 

Figure  4.4.  Time  secies  of  log  $(m2/Hz)  at  selected  frequencies:  106 

0.25  Hz  -  0, 0.5  Hz  -  •,  0.75  Hz  -  1  Hz  - 

Figure  4.5.  Time  series  of  the  RMS  wave  slope  s  (Equation  3.12)  computed  107 

from  die  resistance  wire  wave  gauge  array  data. 

Figure  4.6.  Time  series  of  the  ambient  noise  spectrum  level  N  (dB  re  109 

1  pPa2/l  Hz)  at  selected  WOTAN  frequencies.  Note  the  similarities 
between  these  time  series  and  die  Ujo  tune  series.  In  particular, 
note  the  drop  and  in  N  on  Julian  day  275.  a)  43  kHz.  b)  8.0  kHz. 
c)  140  kHz.  d)  wide  band  [2-14  kHz  or  I]. 

Figure  47.  Directional  wave  spectra  from  0.25  to  1  Hz.  The  sequence  112 

of  plots  demonstrate  die  evolution  of  the  wave  field  in  response  to 
turning  winds  on  Julian  day  278.  The  x-axis  for  the  plots  is  frequency 
in  Hz  and  the  y-axis  is  wave  direction  64.  Dashed  lines  indicate  the 
direction  the  wind  is  blowing  towards  (6wlnd). 

Figure  48.  The  directional  wave  spectra  in  Figure  47  normalized  117 

by  the  spectrum  level  at  the  given  frequency  to  highlight  the 

directional  distribution  of  the  wave  energy  in  response  to  turning 

winds  on  Julian  day  278.  The  x-axis  for  the  plots  is  frequency  in 

Hertz  and  the  y-axis  is  the  wave  direction  84.  Dashed  lines  indicate  the 

direction  the  wind  is  blowing  towards  (6wjn<i). 

Figure  49.  Plot  of  U10  vs.  selected  N  frequencies.  Lines  indicate  124 

regression  lines  obtained  by  various  investigations.  This  study  -  solid 
line,  Evans  et  al  (1984)  -  short  dash,  Leman  et  al  (1984)  -  dash-dot-dot, 

Vagle  et  al  (1990)  -  dash-dot 


Figure  410.  Plot  of  u.  vs.  selected  N  frequencies.  The  figure  shows  a  good  125 

correlation  between  u*  and  N. 

Figure  411.  a)  Plot  of  the  drag  coefficient  Cp  =  (u./Ujq)2  as  a  126 

function  of  (Ujg/Cp).  b)  Plot  relating  the  inverse  wave  age  based  cm  U10 
and  on  u*. 

Figure  412.  Plot  of  N  vs.  the  root  mean  square  wave  amplitude  a.  131 

The  scatter  in  the  data  give  correlation  coefficients  =  0.67  for  all 
the  plots. 

Figure  413.  Plot  of  N  vs.  fl^  The  data  was  high  pass  filtered  to  eliminate  132 

the  contribution  of  swell  energy  to  The  figure  shows  a  significant 
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improvement  in  the  correlation  between  ambient  noise  N  and  the 
RMS  wind  wave  amplitude  a*  (see  Table  43). 

Figure  4.14.  The  high  frequency  surface  wave  energy  and  ambient  133 

noise  N [are  wdtcorrelatecL  a)  Plot  erf  N  vs.  ®(0.5Hz).  b)  Plot  of  N 
vs.  «(0.7SHs).  c)  Plot  erf  N  vs.  «(1.0Hz). 


Figure  4.15.  Plot  of  wave  slope  versus  the  ambient  noise  N  at  selected  137 

frequencies. 

Flgrae  4.16.  Plot  of  wave  slope  s  as  a  function  of  wind  speed  U10,  wind  138 

stress  u>  and  inverse  wave  age  (Ujo/Cp)  and  (u»/Cp). 

Figure  4.17.  The  wind  and  wave  time  series  during  NOBS  can  be  142 

subdivided  into  steady  and  unsteady  conditions.  The  data  in  section  II 
were  considered  steady  while  those  in  sections  I,  III  and  IV  were 
considered  unsteady. 

Figure  4.18.  Plot  of  N(43  kHz)  versus  selected  wind  and  wave  144 

parameters.  The  plots  on  die  left  are  far  data  during  unsteady 
conditions.  The  plots  on  die  right  are  feu  data  during  steady 
conditions. 

Figure  4.19.  Plot  of  U10  versus  selected  wave  parameters.  The  plots  on  145 

the  left  are  for  data  during  unsteady  conditions.  The  plots  on  the  right  are 
for  data  during  steady  conditions. 

Figure  430.  Plot  of  the  non-dimensional  wave  energy  versus  the  inverse  148 

wave  age. 

Figure  431.  Selected  normalized  1-hour  wave  spectra.  The  spectra  were  149 

taken  from  wave  height  records  sampled  at  0000  UT.  The  solid  line 
denotes  Q  for  a  «  0.11  (Equation  417).  a)  Spectra  from  12  day  record, 
b)  Spectra  front  fully  developed  seas. 


Figure  5.1.  Comparisons  of  the  time  series  of  the  log  of  the  total  155 

dissipation  [kg/s3]  computed  using  the  formulas  given  in  Table  5.1. 

a)  Dh  (solid  line)  and  Eh  (short  dash),  b)  Dp  (solid  line)  and  Ep  - 

(short  dash),  c)  Dp  (solid  line)  and  Dh  (short  dash),  b)  Ep  (solid  line) 

and  %  -  (short  dash). 

Figure  53.  Plots  showing  the  spectral  dependence  of  the  KHH  157 

dissipation  dn(f)  and  e^f)  [kg/^Hz].  These  spectral  estimates 
were  computed  from  die  wave  spectra  shown  in  Figure  2.4. 

Numbers  below  die  plots  indicate  die  Julian  day  from  which  the  plots 


15 


List  ef  Figure* 


•retaken. 

Figure  53.  Plot  showing  a  good  correlation  between  0.89<£o£  and  au»g.  161 

The  regression  line  has  a  slope  - 1  and  the  correlation  coefficient 
between  die  two  parameters  is  0.92. 

Figure  5.4.  a)  Rots  showing  die  time  senes  of  die  and  162 

coefficients  0.89 (solid  line),  ccu»g  (short  dash)  (m2/*3].  b)  Plot 

showing  the  relative  magnitude  of  die  difference  between  0.89<£ja£ 
andouug. 

Figure  5.5.  Plots  showing  die  spectral  dependence  of  die  Phillips  165 

dissipation  dp(f)  and  ep(f)  [kg/^Hz].  These  spectral  estimates 
were  computed  from  die  'vave  spectra  shown  in  Figure  2.4. 

Numbers  below  die  plots  indicate  the  Julian  day  from  which  the 
plots  were  taken. 


Figure  5.6.  Plots  of  the  KHH  dissipation  versus  ambient  noise  at  169 

selected  frequencies.  The  dashed  line  indicates  the  regression 
computed  from  die  data  (see  Table  5.1  for  die  coefficients). 

Figure  5.7.  Plots  of  the  reduced  KHH  dissipation  versus  ambient  170 

noise  at  selected  frequencies.  The  dashed  line  indicates  die  regression 
computed  from  die  data  (see  Table  5.1  for  die  coefficients). 

Figure  5.8.  Plots  of  die  Phillips  (1985)  dissipation  Dp  versus  ambient  171 

noise  at  selected  frequencies.  The  dashed  line  indicates  die  regression 
computed  from  the  data  (see  Table  5.1.  for  the  coefficients). 

Figure  5.9.  Plots  of  the  Philips  (1985)  dissipation  Ep  versus  ambient  172 

noise  at  selected  frequencies.  The  dashed  line  indicates  die  regression 
computed  from  die  data  (see  Table  5.2  for  die  coefficients). 

Figure  A.1.  Plot  of  the  Sensitivity  Level  (SL).  203 

Figure  A3.  a)  Plot  of  the  beam  pattern  of  the  hydrophone  at  1  kHz.  204 

b)  Plot  of  the  beam  pattern  of  the  hydrophone  at  3  kHz.  c)  Plot  of 

die  beam  pattern  of  the  hydrophone  at  4  kHz.  d)  Plot  of  the  beam 

pattern  of  die  hydrophone  at  5  kHz.  e)  Plot  of  the  beam  pattern 

of  die  hydrophone  at  10  kHz.  f)  Plot  of  die  beam  pattern  of  the 


hydrophone  at  15  kHz.  g)  Plot  of  die  beam  pattern  of  the  hydrophone 
at  20  kHz. 

Figure  A3.  Geometry  of  Urick's  (If  %)  model  for  the  distribution  of  212 

sound  sources  at  the  ocean  surface. 
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Figure  AA.  Urkk'a  0986)  cos  + model  of  the  vertical  distribution  of  213 

ambient  noise  (dashed  fine)  agrees  weD  with  measurements  of 
Axelrod  et  al  (1965)  at  891  Hz  (thick,  1966). 

Figure  AJL  Geometry  for  the  coordinate  transformation  from  the  true  216 

coordinate  system  to  die  hydrophone  coordinate  system. 

Figure  AA.  Plot  showing  die  array  gain  (AG)  of  die  directional  219 

hydrophone  computed  from  die  available  beam  patterns.  Dipole 
source  assumed  -  O,  omnidirectional  source  arsumed  -  •.  The  line 
indicate*  die  fit  to  the  data  used  in  computing  N(f)  {Equation  3.1]. 
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1  INTRODUCTION  AND  UTEBATUKB  REVIEW 


Since  World  War  II,  there  has  been  oonrideate  research  interest  m  understanding 
the  natural  eouroee  of  underwater  sound  Hie  initial  motivation  of  the  acoustic 
community  was  prfanarilythe  desire  to  chasadotee  the  retain  of  these  noise  sources 
«o  underwater  instruments  could  be  need  mow  effectively.  Underwater  ambient  nofee 
wm  primarily  of  military  and  biological  interest  A»  our  knowledge  of  the 
uMncvraoci  oc  uMflwwcr  souuoi  ^rcw^  oosiio|npiiB9  Mu  •couwaans  D^m  to 
realize  that  ambient  noise  could  be  ueed  to  understand  the  physical  processes  that  are 
responsible  for  generating  underwater  sound  and  ultimately,  increasing  our 
knowledge  of  the  atmosphere  and  the  ocean. 

In  numerous  experiments  in  the  last  decade,  investigators  have  used  passive  acoustic 
measurements  of  underwater  ambient  noiac  to  infer  the  wind  speed  [e.g.,  Shaw,  Watts 
ft  Rossby  (1978)  and  Evans,  Watts,  Halpem  ft  Bouiassa  (1984)],  estimate  the  intensity 
and  distribution  of  breaking  waves  (Fanner  and  Vagle,  1988),  trade  die  spatial  and 
temporal  position  of  breaking  waves  in  die  vicinity  of  die  hydrophones  [e.g..  Ding  ft 
Fanner  (1992),  Crowther  ft  Han&la  (1993)],  infer  the  existence  of  an  acoustic  wave 
guide  formed  by  the  relatively  high  bubble  concentration  near  the  ocean  surface 
(Fanner  ft  Vagle,  1989).  Laboratory  experiments  on  breaking  waves  by  Melville, 
Loewen,  Felizardo,  Jessup  ft  Buckingham  (1988)  and  by  Loewen  ft  Melville  (1991a) 
suggest  that  the  dissipation  of  surface  wave  energy  in  the  ocean  may  also  be  inferred 
from  measurements  of  ambient  noise.  This  hypothesis  is  further  explored  in  this 
work. 


It  is  now  believed  that  die  sound  generated  by  bubble  mechanisms  in  breaking  waves 
is  die  primary  source  of  the  wind-dependent  ambient  noise  in  the  ocean  [Kerman 


(IMS,  1992)].  Yet  while  decade*  of  field  measurements  have  shown  a  well-correlated 


powre  tow  relationship  between  underwater  ambient  notoe  and  wind  speed 
[Knudsen,  Alford  Jc  Boding  (1948),  Urick  (1966)],  very  lew  measurements  have  tried 
to  establish  a  similar  correlation  between  ambient  noiae  and  wave  field  parameters. 
The  lew  measurements  that  have  attempted  to  do  so  show  that  ambient  noise  does 
not  correlate  weB  with  the  significant  wave  height  [PenhaQow  &  Dietz  (1964),  Perrone 
(1969),  Fanner  6c  Lemon  (1984)].  This  to  remarkable  given  that  breaking  waves  are 
more  directly  related  to  the  characteristics  of  the  wave  field  rather  than  die  wind.  We 
believe  that  bus  could  be  due  to  the  fact  that  swell  and  the  longer  wave  components 
which  do  not  break  contribute  a  significant  proportion  to  the  variance  of  the  wave 
field.  The  laboratory  measurements  of  Melville  and  his  coworkers  [e.g.,  Melville  et  al 
(1988),  Loewen  4c  Melville  (1991a),  Lamarre  4c  Melville  (1991)]  suggest  that  both  the 
sound  generation  and  the  dissipation  m  breaking  waves  are  coupled  by  the  air 
entrainment  process  and  that  this  to  why  laboratory  experiments  have  shown  that 
these  two  parameters  correlate  with  each  other. 

The  aim  of  this  thesis  to  to  investigate  whether  comparable  correlations  hold  in  the 
ocean.  To  this  end,  we  have  conducted  an  experiment  measuring  ambient  noise, 
wind,  and  die  surface  wave  field  from  a  floating  platform  (RP  Flip)  off  fire  coast  of 
Oregon.  These  measurements  will  be  used  to  investigate  die  correlation  between 
ambient  sound  and  a  number  of  surface  wave  parameters,  especially  surface  wave 
energy  dissipation  which  to  most  directly  related  to  wave  breaking.  In  addition,  we 
will  use  these  measurements  to  obtain  empirical  relationships  between  these 
parameters. 

Before  proceeding  to  a  discussion  of  the  field  experiment  and  the  data,  we  will,  in  this 
chapter,  review  file  available  literature  on  ambient  noise  measurements  and  their 
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correlation  with  wind  and  wave  parameter.  We  will  then  discuss  the  mechanisms  for 
sound  generation  associated  with  breaking  wave*.  Later,  we  will  examine  the  results 
of  some  recent  laboratory  experiments  on  breaking  waves.  Finally,  we  will  discuss  the 
literatim  an  surface  wave  energy  dissipation. 

U  Ambient  noise  and  the  wind  speed. 

Pioneering  underwater  sound  measurements  by  Knudsen,  Alford  &  Ending  (1948) 
during  World  War  II  showed  that  the  noise  spectrum  level  N(f)  in  the  100  Hz  -  25  kHz 
range  correlates  with  sea  state.  The  study  resulted  in  what  is  now  commonly  referred 
to  as  the  Knudsen  curves',  a  family  of  curves  showing  a  slope  of  approximately  -17  to 
-20  dB  per  decade  within  this  frequency  range  as  a  function  erf  wind  speed.  The 
Knudsen  curves  have  been  verified  by  other  investigators  over  the  last  few  decades 
[see  Wens  (1962)  and  Urick  (1986)  for  a  review].  Figure  1.1  reproduces  a  figure  from 
Wenz  (1962)  which  summarizes  foe  sources  and  sound  spectrum  levels  of  the 
prevailing  ambient  noise  in  the  ocean.  The  Knudsen  curves  in  this  figure  are 
expressed  in  terms  of  sea  state  which  was  commonly  used  at  die  time.  Table  1.1  gives 
die  equivalent  wind  speed  for  a  given  sea  state.  Although  this  figure  is  three  decades 
old,  it  is  still  useful  in  providing  a  general  idea  of  the  sound  levels  and  sources  in  the 
ocean. 

Functional  expressions  between  ambient  noise  and  wind  speed  were  first  proposed  by 
Crouch  8c  Burt  (1972)  using  field  measurements  by  Piggott  (1964)  and  Perrone  (1969). 
They  found  that  this  relationship  can  be  described  by  die  empirical  relation 
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Figure  l.l.  Sources  and  spectrum  level  of  deep  water  ambient  noise  (Wen z,  1962). 
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Tabic  1.1.  Approximate  relationship  between  sea  state  and  wind  speed  [adapted 
from  Wens  (1962)]. 


Beaufort 
sea  state 

Mean  wind  speed 
(m/s) 

Sea  criteria 

0 

<0.5 

Mirror-like 

1 

1.0 

Ripples 

2 

2.5 

Small  wavelets 

3 

45 

Large  wavelets,  scattered  whitecaps 

4 

7.0 

Small  waves,  frequent  whitecaps 

5 

10.0 

Moderate  waves,  many  whitecaps 

6 

12.5 
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15.5 

Heaped-up  sea,  blown  spray,  streaks 

8 

19.0 

Moderately  high  long  waves,  spindrift 

N(f)=N,(f)+20n(f)logU  (1.1)» 


where  f  is  the  frequency  at  which  the  measurements  were  made,  Nj(f)  is  the  noise 
spectrum  level  at  1-knot  wind  speed,  U  is  the  wind  speed  in  Inots,  and  n(f)  is  a 
coefficient  of  order  one.  N  is  the  wind  dependent  noise  spectrum  level  (dB  re  1 
mbar2/!  Hz) 


N(f)  =  10logTi^  . 
.PSrf  . 


0-2) 


where  p2(f)  is  the  ambient  noise  power  spectrum.  Crouch  &  Burt  (1972)  used  a 
reference  pressure  p^  =  1  mbar2  which  is  different  from  the  current  System 
International  (SI)  unit  p^f  =  1  pPa2  more  commonly  used  today.  Note  that  Equation 

1  All  common  logarithms  used  in  this  work  are  denoted  by  'log'  while  natural  logarithms  are 
denoted  by  'In’. 
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WOTAN  (Un)  m/fc 


Figure  1.2.  A  plot  comparing  the  wind  speed  measured  by  a  vector-averaging  wind 
recorder  (VAWR)  type  anemometer  3.5  meters  above  the  sea  surface  and  by  a 
WOTAN  at  4.3  kHz  ambient  noise  frequency  (Evans  et  al,  1984). 
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Figure  13.  Plot  of  log  U  versus  the  ambient  noise  level  at  8  and  14.5  kHz.  Note  that 
when  U  exceeds  10  meters/second,  the  departure  from  the  linear  relationship 
between  N  and  log  U  becomes  significant  (Fanner  &  Lemon,  1984). 
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1. 1  is  not  in  dimensionless  form.  Crouch  k  Burt  obtained  values  of  n(f)  ranging  from 
0.8  to  2.8.  By  fitting  data  obtained  by  other  investigators  [Cato  (1976),  Shooter  8c 
Gentry  (1981)/  and  Burgess  k  Kewley  (1983)],  Urick  (1986)  obtained  values  of  n 
ranging  from  0.6  to  1.5.  The  differences  in  n  may  be  due  to  toe  differences  in  the  way 
toe  investigators  sampled  the  data,  the  site  dependence  of  toe  U  vs.  N  relationship, 
and  toe  differences  in  the  sensitivity  of  N  at  each  frequency  to  changes  in  wind  speed. 
The  degree  of  scatter  in  toe  data  appears  to  be  larger  at  toe  lower  frequencies 

[f  ■  0(100  Hz)]  and  the  wind  dependence  at  the  higher  frequencies  [f  >  0(1  kHz)] 
appears  to  be  stranger. 

In  toe  last  decade,  Ok  relationship  between  ambient  noise  and  wind  speed  has  been 
used  as  a  method  of  predicting  wind  speed  from  ambient  noise.  Field  measurements 
by  Shaw,  Watts  k  Rossby  (1978)  showed  that  this  relationship,  described  by  an 
expression  of  toe  form 

l</gU  =  m(f)N(f)+n(f)  (1.3) 

can  be  used  as  a  means  of  measuring  the  wind  speed  using  hydrophones.  In  Equation 
13,  U  is  the  wind  speed  in  meters  per  second,  f  is  the  frequency  in  Hertz,  m(f)  and 
n(f)  are  coefficients,  and  the  noise  spectrum  level  N(f)  is  in  dB  re  1  pPa2/l  Hz.  This 
technique,  commonly  referred  to  as  Weather  Observations  Through  Ambient  Noise  or 
WOT  AN  (Evans,  Watts,  Halpem  k  Bourassa,  1984),  has  been  used  in  several 
measurements  over  the  last  decade  [Evans  et  al  (1984),  Lemon,  Fanner  k  Watts  (1984) 
and  Vagle,  Large  k  Farmer  (1990)].  The  development  of  WOT  AN  instruments  is 
described  by  Hill  (1984).  In  the  subsequent  discussion  in  this  work,  equations  relating 
U  and  N  in  toe  form  of  Equation  1.3  will  be  referred  to  as  WOT AN  equations’.  Figure 

1.2  shows  a  figure  taken  from  Evans  et  al  (1984)  which  compares  wind  speed 
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Table  1J.  Some  coefficients  lor  flu  WOT  AN  eiyrmion  and  die  equivalent  power 
lawa  obtained  from  previous  field  measurements. 

2 

logU  -  m(f)N(f)+ n(f)  4-  -  mU" 
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8.0 

0.0419 
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1.5 

2.4 

14.5 

0.0434 

-1.217 

0.64 
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Lemon.  Farmer  &  Watts  (19841 

site  A 

4.3 

0.0336 

-1.045 

1.29 

3.0 

8.0 

■ES53X 

-0.925 

0.64 

3.0 

14.5 

0.0341 
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0.38 
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Lemon.  Farmer  &  Watts  (1984) 

siteB 

4.3 

WKSSSM 

-1.64 

4.6 

2.2 

14.5 

0.0465 

-1.32 

0.69 

2.2 

Vagje,  Large  &  Farmer  (1990) 

1  8.0 

0.0378 

-1.006 

0.46 

2.7 

measurements  made  by  WOTANs  and  by  vector  averaging  wind  recorders  (VAWR) 
3.5  meters  above  the  sea  surface.  Hie  plot  shows  good  agreement  between  the  two. 

Table  1.2  summarizes  die  coefficients  m  and  n  in  Equation  1.3  obtained  by  several 
investigators.  Their  studies  confirmed  the  validity  of  Equation  1.3  but  showed  that  m 
and  n  are  site  and  frequency  dependent.  The  table  shows  that  although  the  regression 
of  log  U  and  N  is  linear,  there  can  be  significant  differences  in  the  equation  of  the 
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regression  lines  among  different  sites.  In  particular,  the  ambient  noise  measurements 
reported  by  Lemon  et  al  (1984)  and  given  in  Table  1.2  were  both  in  Queen  Charlotte 
Sound  on  the  west  coast  of  Canada  yet  their  WOT  AN  Equation  coefficients  are  very 
different  The  water  depth  at  the  measurement  sites  are  comparable  (265  meters  at  Site 
A  and  287  meters  at  site  B).  Site  A  is  doee  to  the  shore  while  site  B  is  at  the  center  of 
the  channel.  This  suggests  that  die  surrounding  topography  has  a  significant  impact 
on  either  die  noise  generating  mechanisms  or  on  how  die  sound  reaches  die 
hydrophone. 

It  is  important  to  point  out  however  that  although  the  WOT  AN  coefficients 
summarized  in  Table  1.2  differ  from  site  to  site,  they  do  give  wind  speed  estimates 
that  do  not  differ  significantly  from  each  other  within  the  4  to  12  m/s  range.  This  is 
demonstrated  by  Table  1.3  which  shows  a  ±0.5  m/s  difference  when  U  =  5  m/s  and 
±1.5  m/s  when  U  *  10  m/s.  If  we  restrict  ourselves  to  the  4.3  kHz  data,  the  difference 
at  10  m/s  is  less  than  ±1.0  m/s.  The  error  is  more  substantial  at  higher  U  and  higher 
frequencies. 


Sound  absorbed  by  bubbles  generated  by  breaking  waves  has  been  shown  to  cause 
significant  departures  from  linearity  at  higher  ambient  noise  frequencies  and  higher 
wind  speeds.  Figure  1.3  shows  a  plot  of  log  U  versus  N  at  8  and  14.5  kHz  measured 
by  Fanner  &  Lemon  (1984)  at  Queen  Charlotte  Sound.  The  figure  shows  that  when  die 
wind  speed  U  exceeds  approximately  8  meters/second,  there  is  a  clear  departure  from 
a  linear  relationship  between  log  U  and  N.  Fanner  A  Lemon  (1984)  found  that  this 
discrepancy  is  more  pronounced  at  higher  frequencies  and  at  higher  wind  speeds. 
They  suggested  that  at  higher  wind  speeds,  waves  break  more  frequently  and  more 
intensely.  These  breaking  waves  generate  bubbles  and  the  concentration  of  these 
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Table  1.5.  Some  typkal  wfakd  speeds  predicted  by  the  WOT  AN  equation  and 
coefficients  given  la  Table  13. 


N(43kHz) 
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LFW  Site  A 
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48 

48 

46 

6.4 

56 

10.4 

8.4 

12.9 

60 

15.3 

11.4 
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1  ■■■■■■■■■ 

LFW  Site  A 

LFW  Site  B 

44 

49 

4 

53 

51 

9.9 

15 

113  ! 

55 

148 

11.6 

173 

LFW  -  Lemon,  Fanner  fc  Watts  (1984) 
N(f)  is  in  dB  re  1  pPia2/!  Hz 


bubbles  increases  with  increasing  wind  speed.  These  bubbles  persist  near  the  ocean 
surface  and  absorb  the  sound  generated  by  breaking  waves. 


Vagle  et  al  (1990)  examined  whether  WOT  AN  coefficients  obtained  at  me  site  can  be 
used  to  predict  die  wind  speed  at  other  sites.  They  compared  the  WOT  AN  coefficients 
obtained  during  die  Frontal  Air-Sea  Interaction  Experiment  (FASINEX)  at  a  deep 
water  site  500  km  southwest  of  Bermuda  with  the  wind  and  ambient  noise  data 
obtained  at  three  shallow  water  and  me  deep  water  site  They  concluded  that  the 
FASINEX  WOT AN  equations  for  ambient  noise  less  than  8  kHz  can  reliably  predict 
the  wind  speed  from  the  ambient  noise  at  another  deep  water  site  as  long  as  the 
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instruments  are  carefully  calibrated  and  extraneous  noise  is  screened  out  They  also 
concluded  that  the  variations  in  the  WOTAN  equations  at  sites  dose  to  the  shore  such 
as  those  reported  by  Lemon  et  si  (1964)  were  due  to  the  effect  of  wind  direction  on 
letch  (and  consequently  the  number  of  breaking  waves),  the  effect  of  the  bottom,  die 
geometry  of  die  surrounding  topography,  and  the  influence  of  industrial,  shipping 
and  other  man-made  noise. 


Although  the  WOTAN  expression  is  dimensionally  inhomogeneous,  we  can 
nonetheless  rearrange  Equation  1.2  to  relate  U  and  N  in  terms  of  a  power  law  of  the 
form 


(1.4) 

Pnf 


where  m  » lO-n/lOm  and  »  *  l/(10m).  Based  on  the  values  of  Table  1.1,  we  see  that 

p2- 1/2-3.  (1.5) 


UZ  Ambient  noise  and  the  wave  field. 


Although  foe  linear  relationship  between  N(f)  and  log  U  is  well  established  by 
decades  of  measurements,  die  physics  behind  this  correlation  cannot  be  gleaned  from 
die  WOTAN  equations  alone.  Theoretical  studies  and  laboratory  experiments  cm 
bubbles  and  bubble  clouds,  which  are  reviewed  in  the  next  section,  have  led  to  a 
growing  consensus  that  the  sound  generated  by  breaking  waves  is  responsible  for  the 
wind  dependence  of  ambient  sound  in  die  Knudsen  range.  However,  attempts  at 
establishing  a  direct  and  predictable  relationship  between  ambient  noise  and  surface 
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Table  1.4.  Correlation  coefficients  between  N,  H,  and  U.  (Penhallow  Jc  Diets,  1964) 


steady  winds 
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0.66 

0.79 

0.42 

1  u 

0.88 

0.67 

0.47 

wave  field  parameters  have  been  few  and  largely  unsuccessful.  Perrone  (1969)  found 
that  the  correlation  coefficient 


<-l  M 


between  U  and  N  is  0.80  while  die  correlation  coefficient  between  die  significant  wave 
height  H,  and  N  is  only  0.60  in  the  200  Hz  -  2.8  kHz  frequency  band.  Fanner  &  Lemon 
(1984)  obtained  comparable  correlation  coefficients  for  Hg  and  ambient  noise  at  43, 8.0 
and  14.5  kHz.  By  classifying  measured  wind  velocities  as  steady,  increasing  and 
decreasing,  Penhallow  &  Dietz  (1964)  found  that  although  U  tracked  ambient  noise 
well  under  all  these  conditions,  Hg  correlated  well  rally  under  steady  conditions.  Table 
1.2  summarizes  their  results.  [Both  Perrone  (1969)  and  Penhallow  &  Dietz  (1964) 
computed  the  linear  correlation  coefficient  between  N(f)  and  U  instead  of  log  U.] 
Perhaps  due  to  die  relatively  poor  correlation  between  N  and  Hg,  none  of  these 
studies  published  the  regression  coefficients  between  the  two  variables. 
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The  poor  correlation  between  N  and  Hg  may  appear  surprising  given  that  the  ambient 
noise  source*  are  expected  to  be  direct  products  of  the  wave  field  itself.  One  should 
then  expect  that  the  correlation  with  wave  parameters  should  be  better  than  the  wind 
speed  correlation.  This  apparent  discrepancy  may  partly  be  explained  by  the 
relatively  large  distance  between  the  wave,  wind  snd  N  instruments  in  die  above 
experiments.  Perrone's  (1969)  wave  and  wind  measurements  were  made  30  miles 
from  the  hydrophone  site.  Fanner  k  Lemon  (1984)  deployed  their  wave  buoy  'near' 
their  WOTAN  although  it  is  not  dear  what  the  exact  distance  is.  Penhallow  4c  Dietz 
(1964),  who  made  measurements  dose  to  shore,  had  the  anemometers  and  wave 
gauges  installed  off  a  dock  white  the  hydrophone  was  at  12  meters  depth,  300  meters 
offshore. 


Perhaps  the  more  plausible  explanation  is  that  the  presence  of  swell  increases  Hg 
without  proportionally  increasing  the  incidence  of  wave  breaking.  Since  swell  is 
generated  by  storms  far  from  die  observation  site,  it  is  possible,  particularly  in  the 
open  ocean  with  large  fetch  lengths,  for  Hg  to  be  large  even  on  a  calm  day  with  low 
wind  speeds.  Same  recent  studies  [(Thorpe,  1992)  and  (Melville,  1992)]  suggest  that  it 
is  die  wind  waves  at  frequencies  greater  than  die  wind  sea  spectral  peak  that  are 
breaking.  If  their  models  are  correct,  then  wave  components,  including  wind  waves 
and  swell,  which  do  not  break  make  a  significant  contribution  to  die  variance  of  the 
sea  surface  displacement  but  not  to  the  noise.  If  wave  breaking  is  indeed  primarily 
responsible  for  the  ambient  noise  in  die  ocean,  then  correlations  with  N  should  be 
made  with  wave  parameters  that  quantify  wave  breaking. 
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Figure  1.4.  Ambient  noise  spectrum  level  during  a  breaking  event  detected  by  a 
hydrophone  14  meters  below  sea  level.  The  background  spectrum  immediately  before 
breaking  shown  in  the  first  plot  is  reproduced  by  the  dashed  lines  in  the  second  and 
third  plots.  The  90%  confidence  level  is  indicated  by  the  vertical  bars  in  each  plot 
(Fanner  k  Vagle,  1989). 
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In  the  last  decade,  there  have  been  a  significant  number  erf  investigations  that  have 
looked  kilo  the  acoustic  properties  of  breaking  waves.  There  appears  to  be  a  growing 
consensus  that  the  prevailing  noiee  level  in  the  Kmidsen  range  is  primarily  due  to  the 
wove  breaking  mechanieat*  associated  with  bubble  creation,  coalescence,  splitting  and 
collective  bubble  cloud  oscillation. 

Feimer  4c  Vagle  (1989)  described  experiment  off  Vancouver  Island  in  140  m  of  water 
and  at  a  tile  500  km  southwest  of  Bermuda  in  4000  m  water.  By  comparing  video 
recordings  with  simultaneous  acoustic  measurements,  they  demonstrated  that 
breaking  waves  can  generate  ambient  noise  at  frequencies  as  low  as  100  Hz.  Figure 
1.4,  taken  from  their  paper,  shows  that  at  frequencies  within  the  Kmidsen  range, 
breaking  waves  can  generate  sound  levels  significantly  higher  than  the  background 
level  HoOett  (1989)  reported  measurements  of  breaking  wave  generated  noise  as  low 
as  30  Hz  in  the  Mediterranean  Sea. 

By  comparing  die  WOT  AN  measurements  in  Georges  Strait,  British  Columbia  at  4.3 
kHz  with  a  model  of  sound  generation  by  randomly  distributed  breaking  waves. 
Farmer  4c  Vagle  (1988)  showed  that  changes  in  die  dominant  wave  period,  the  mean 
breaker  spacing,  and  die  mean  breaking  wave  intensity  can  be  inferred  from  changes 
in  die  ambient  noise  spectrum  level  Their  model  also  suggests  that  the  underwater 
ambient  noise  level  is  proportional  to  die  frequency  and  intensity  of  wave  breaking  in 
die  vicinity  of  die  hydrophone.  The  importance  of  their  work  is  that  it  suggests  that 
information  about  the  characteristics  of  the  wavefield  and  die  distribution  of  breaking 
waves  can  be  obtained  from  measurements  of  ambient  sound. 
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Several  theories  an  how  breaking  wave*  generate  sound  have  been  proposed  [tee 
Kennan  (1964)  for  a  review  of  early  theories}.  These  include  splashing  water  droplets 
[(Frans;  1999)  and  Wilson  (1963)},  cavitation  noise  (Purduev,  1966),  nonlinear 
interaction  of  water  waves  [Longuet-Higgins  (1951),  Marsh  (1963),  and  Kuo  (1968)], 
and  atmospheric  turbulence  (Fitzpatrick  4c  Strasberg,  1956).  Kennan  (1984)  reviewed 
these  theories  and  discussed  some  of  their  major  shortcomings. 

Although  die  oscillation  of  bubbles  generated  by  breaking  waves  has  been  considered 
as  a  plausible  ambient  noise  mechanism  for  several  decades  (Wenz,  1962),  more  recent 
laboratory  experiments  and  theoretical  modeling  have  suggested  that  this  is  the 
dominant  mechanism  responsible  for  the  characteristics  of  wind-dependent  ambient 
noise  at  higher  frequencies  (f  >  500  Hz).  In  addition,  other  recent  investigations 
reviewed  later  in  this  section  suggest  die  possibility  that  die  collective  oscillation  of  a 
breaking  wave  bubble  doud  could  be  responsible  for  a  significant  proportion  of  the 
noise  spectrum  level  below  500  Hz.  An  extensive  survey  of  the  current  research  an 
wave  breaking  sound  sources  and  mechanisms  can  be  found  in  the  proceedings  of  the 
1987  NATO  workshop  in  Lend,  Italy  (Kennan,  1988)  and  die  1990  Conference  on  the 
Natural  Physical  Sources  of  Underwater  Sound  in  Cambridge,  UK  (Kennan,  1992). 

In  a  now  classic  paper,  Minnaert  (1933)  studied  die  sound  generated  by  die  resonant 
volumetric  oscillation  of  a  newly  created  bubble  in  water.  He  showed  that  if  adiabatic 
compression  is  assumed,  then  die  resonant  frequency  4  of  this  mode  of  oscillation  is 
given  by 


4 


(1.6) 
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where  r0  is  the  mean  radius  of  the  bubble,  y  «  cr/cv  is  the  ratio  of  the  specific  heats  of 
the  entrained  gas,  P  is  the  ambient  pressure  at  the  bubble  and  pw  is  the  density  of 
water.  For  an  ideal  gag  under  adiabatic  conditions,  y-  1.4. 

In  linear  wave  theory,  the  fluctuations  in  the  sound  pressure  p  is  related  to  the 
velocity  of  the  fluid  particles  v  by  the  equation  (Urick,  1975,  p.  12) 

P-pev  (1.7) 

where  p  is  the  density  of  the  fluid  and  c  is  the  propagation  velocity  of  the  pressure 
wave.  The  value  of  the  specific  acoustic  impedance  pc  in  sea  water  is  1.5  x  10s  g/cnA 
and  42  g/cm^  for  air.  The  large  difference  in  pc  between  air  and  water  implies  that 
for  an  underwater  sound  wave,  the  pressure  is  essentially  zero  (p  *  0)  at  the  air-water 
interface  (pressure  release  condition). 

The  resonant  frequency  of  a  bubble  near  the  water  surface  is  modified  by  the  pressure 
release  condition  at  the  air-water  interface.  For  a  wave  incident  on  the  air-sea 
interface,  the  phase  of  the  reflected  acoustic  energy  is  shifted  by  180°  to  satisfy  the 
pressure  release  condition.  A  bubble  oscillating  near  an  acoustically  smooth  interface 
as  a  monopole  nt  a  depth  L/2  below  the  water  surface  creates  an  image  monopole 
source  L/2  above  the  water  surface  oscillating  out  of  phase  with  the  real  bubble.  In 
short,  the  far  Add  behavior  of  a  bubble  near  the  sea  surface  can  be  represented  by  a 
dipole  whose  resonant  frequency  is  higher  than  that  of  the  bubble  resonating  as  a 
monopole.  If  the  bubble  radius  rQ  «  L,  then  the  resonant  frequency  of  the  dipole  f^ 
can  be  approximated  by  (Longuet-Higgins,  1990) 


(18) 


Figure  1.5.  The  sound  radiated  by  a  typical  oscillating  bubble  generated  by  a 
laboratory  spilling  breaker  can  be  described  as  an  exponentially  damped  sinusoid 
(Medwin  &  Beaky,  1989). 
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From  observations  of  quasi-steady  spilling  wave  breakers  in  the  laboratory  using  a 
high  speed  movie  camera  synchronized  with  die  soundtrack.  Banner  A  Cato  (1988) 
found  that  the  noise  bunts  detected  by  a  hydrophone  corresponded  to  the  bubble 
formation  at  the  leading  edge  of  the  spill,  the  impact  and  coalescence  of  bubbles,  and 
foe  splitting  of  bubbles.  They  also  suggested  that  foe  range  of  sound  frequencies  they 
detected  corresponded  to  the  resonant  frequencies  of  the  range  of  bubble  sizes  they 
observed.  A  more  careful  experiment  by  Pumphrey  A  Ffowcs  Williams  (1990)  who 
studied  die  free  oscillations  of  individual  bubbles  generated  by  water  spilling  from  a 
trough  and  by  air  forced  through  a  submerged  hypodermic  needle  showed  results 
that  are  consistent  with  Banner  A  Cato's  (1988)  experiments.  Laboratory 
measurements  of  spilling  breakers  in  an  anechoic  tank  (Med win  A  Beaky,  1989) 
showed  that  the  oscillation  of  these  bubbles  ran  bp  characterized  as  damped  sinusoids 
(Figure  1.5)  which  suggests  that  foe  bubbles  were  forced  initially  and  then  oscillated 
freely-  Medwin  A  Beaky  (1989)  then  proposed  that  the  sound  from  a  breaking  wave  is 
due  to  the  incoherent  addition  of  the  sound  generated  by  foe  individual  bubbles 
generated  by  foe  breaking  wave. 

Using  Medwin  A  Daniel's  (1990)  measurements  of  bubble  size  distribution  in 
laboratory  breaking  waves  and  by  assuming  that  foe  oscillating  bubbles  radiate 
acoustic  pressure  as  a  dipole,  Loewen  A  Melville  (1991b)  proposed  a  semi-empirical 
sound  source  model  for  gently  spilling  waves  that  closely  predicted  the  spectral 
magnitude  and  shape  of  foe  measured  sound  spectrum.  The  largest  bubble  Medwin  A 
Daniel  (1990)  observed  had  a  7.4  mm  radius  and  a  440  Hz  resonant  frequency. 

Loewen  A  Melville  (1991b)  speculated  that  bigger  bubbles  may  be  generated  by 
breaking  waves  in  the  ocean  and  these  can  radiate  sound  at  even  lower  frequencies. 
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For  a  large  bubble  cloud  generated  by  intense  breaking  waves,  it  may  be  appropriate 
to  treat  it  is  a  fluid  continuum  whose  density  is  comparable  to  the  surrounding  fluid 
but  whose  compressibility  is  significantly  higher.  Carey  k  Browning  (1988)  suggested 
that  the  oscillation  of  such  a  bubble  cloud,  whose  dimensions  are  significantly  larger 
than  the  largest  individual  bubbles  in  the  ocean,  might  account  for  the  low  frequency 
wind-dependent  noise.  In  a  subsequent  paper,  Carey  8c  Fitzgerald  (1992)  showed  that 
a  spherical  collection  of  bubbles  can  oscillate  at  a  resonant  frequency 
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(1.9) 


where  r0  is  the  radius  of  die  bubble  cloud,  y  *  1  is  the  ratio  of  specific  heats  under 
isothermal  conditions,  P  is  the  ambient  pressure,  pw  is  the  density  of  surrounding 
fluid,  and/is  the  mean  void  fraction  of  the  bubble  cloud.  Note  that  this  equation  is 
similar  to  the  Minnaert  (1933)  formula  for  the  resonant  oscillation  of  spherical  bubbles 
(Equation  1.6).  For  a  spherical  bubble  cloud  in  water  of  radius  rG  -  0.1  meters, 

/=  10  %  air,  and  an  ambient  pressure  p  =  101  kPa,  Equation  1.9  predicts  a  resonant 
frequency  fQ  *  90  Hz.  If  the  bubble  cloud  has  a  radius  rG  =  1  meter,/  =  1  %  air. 
Equation  1.9  predicts  a  resonant  frequency  f0  -  30  Hz. 

A  more  general  theory  for  sound  generation  by  collective  oscillations  of  bubble  clouds 
was  proposed  by  Prosperetti  (1988)  and  Lu,  Prosperetti  k  Yoon  (1990).  Luetal  (1990) 
outlined  their  approach  and  used  it  to  predict  the  resonant  collective  oscillation 
frequencies  of  bubbly  regions  of  various  geometries.  In  particular,  they  showed  that 
their  method  successfully  predicted  the  frequency  of  the  resonant  oscillations  of  a 
cylindrical  bubble  column  generated  in  the  laboratory  (Yoon,  Crum,  Prosperetti  k  Lu, 
1991). 
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A  6-0.40 


Figure  1.6.  Hot  showing  the  c£/served  peak  frequency  of  the  pressure  spectrum 
from  breaking  waves  and  the  predicted  resonant  frequency  f r  for  a  cylindrical  bubble 
cloud  computed  using  the  method  described  by  Lu  et  al  (1990).  [Loewen  &  Melville, 
1992]. 
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Laboratory  experiments  using  bubble  columns  (Yoon  et  al,  1991)  and  bubble  plumes 
(Kolaini,  Roy  4c  Crum,  1991)  suggest  that  this  mechanism  might  explain  the  low 
frequency  pressure  fluctuations  corresponding  to  modes  of  the  bubble  plumes 
detected  by  a  hydrophone.  Kolaini  et  al  (1991)  created  bubble  plumes  by  dropping  a 
cylindrical  volume  of  water  into  a  tank  of  water.  They  found  that  substructures 
consisting  of  a  sphere  of  bubbly  water  of  high  void  fraction  break  off  from  die  main 
plume  and  oscillate  at  frequencies  as  low  as  tens  of  Hertz.  The  resonant  frequencies 
they  observed  matched  the  lowest  eigenfrequendes  of  die  substructures  if  a  void 
fraction  of  40%  within  die  plumes  was  assumed. 


Recently,  Loewen  4c  Melville  (1992)  [see  also  Loewen  (1991)]  observed  low  frequency 
peaks  of  die  order  of  tens  erf  Hertz  in  die  noise  spectrum  of  two-dimensional  and 
three-dimensional  plunging  waves  in  die  laboratory.  By  using  the  measurements  of 
void  fraction  of  breaking  wave  bubble  clouds  (Lamarre  4c  Melville,  1992)  and  the 
theory  of  collective  bubble  cloud  oscillations  by  Lu  et  al  (1990),  they  found  that  the 
low  frequency  peak  in  the  spectrum  of  pressure  fluctuations  picked  up  by  a 
hydrophone  matches  the  resonant  frequency  of  collectively  oscillating  bubble  clouds 
predicted  by  Lu  et  al  (1990)  [see  Figure  1.6]. 


1.4  Breaking  wave  sound  and  dissipation. 

In  the  last  few  years,  our  research  group  has  conducted  a  series  of  laboratory 
experiments  to  establish  quantitative  correlations  between  the  acoustics  and  the 
dynamical  characteristics  of  breaking  waves.  The  relationship  between  the  energy 
dissipated  by  breaking  waves  and  the  sound  radiated  was  initially  examined  by 
Melville,  Loewen,  Felizardo,  Jessup  k  Buckingh*  t  (1988)  for  controlled  breaking 
waves  in  the  laboratory.  This  experiment,  as  well  as  the  subsequent  investigation  by 
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Figure  1.8.  The  fraction  of  energy  dissipated  by  a  breaking  event  D  is  proportional  to 
die  initial  volume  per  unit  width  of  air  entrained  by  a  breaking  wave  Vc.  The  figure 
shows  laboratory  results  from  three  different  wave  packets  of  various  wave  slope 
(adapted  from  Loewen,  1991). 
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Loewen  Ac  Melville  (1991*)  using  •  tagger  set  of  wave  packets,  established  flat  over  a 
limited  range  of  parameters,  the  acoustic  energy  radiated  by  a  breaking  event  is 
proportional  to  the  surface  wave  energy  dissipated  by  the  event  (Figure  1.7).  In 
addition,  Loewen  4c  Melville  (1991a)  showed  that  approximately  10"®  of  the  dissipated 
energy  is  converted  to  acoustic  energy  and  that  the  energy  of  the  radiated  sound  is 
proportional  to  the  duration  of  the  event  By  relating  the  electrical  conductivity  of  an 
air-water  mixture  to  die  void  fraction  of  the  medium,  Lamarre  4c  MdviUe  (1991) 
mapped  die  spatial  and  temporal  void  fraction  distribution  in  laboratory  breaking 
waves.  By  studying  die  evolution  of  die  bubble  cloud  generated  by  the  breaking 
wave,  they  found  that  a  significant  proportion  of  the  dissipated  energy  (up  to  50%)  is 
due  to  work  done  by  the  liquid  in  entraining  air  against  buoyancy.  Using  data  from 
similar  experiments,  Loewen  (1991)  (see  also  Loewen  and  Melville  (1992)]  showed 
riwt  the  energy  dissipated  by  a  breaking  event  is  proportional  to  die  initial  volume  of 
air  entrained  (Figure  1.8).  This  is  consistent  with  die  Loewen  4c  Melville  (1991b)  dipole 
sound  source  model  which  suggested  that  for  similar  bubble  size  distributions,  die 
total  acoustic  energy  radiated  is  proportional  to  die  volume  of  air  entrained. 

These  laboratory  studies  clearly  suggest  that  sound  generation  and  the  dissipation  of 
imytmniral  wave  energy  are  coupled  through  die  air  entrainment  process.  We  believe 
that  this  is  die  reason  why  the  two  correlate  well  in  the  laboratory  and  we  expect  that 
these  two  parameters  should  also  correlate  well  in  die  field. 

1.5  Wind  wave  dissipation. 

Direct  measurements  of  dissipation  due  to  wave  breaking  in  the  field  are  not  yet 
feasible  although  some  progress  is  being  made  to  achieve  this  (Agrawal  et  al,  1992). 
However,  dissipation  can  be  estimated  using  the  energy  budgets  arising  from  wind- 


wave  evolution  model*.  If  tenns  other  than  disjpatian  can  be  effectively  measured  or 
modeled,  then  the  dissipation  can  be  inferred  through  an  energy  balance. 

We  can  describe  the  distribution  of  the  sea  surface  elevation  a(x)  using  the 
wavenumber  spectrum  (Philos,  1985) 

Y(k)-^  £a(x)a(**r)e"**d*  (1.10) 

which  is  given  in  terms  of  the  covariance  of  the  wave  height  between  two  points  at  x 
and  x+r.  The  evolution  of  the  surface  wave  field  can  then  be  represented  by  the 
energy  transfer  equation  (Phillips,  1985) 

P.t5^-P.*”+P.<(e,+v)  V.1-— Vj.^+wOO-dft*  0.11) 

where  pw  is  the  density  of  water,  is  the  group  velocity  and  V  is  the  current  velocity. 
The  three  forcing  tenns  on  the  right  hand  side  are  die  net  spectral  flux  from  resonant 
wave-wave  interactions  -Vk  •  T(k),  die  spectral  wind  farcing  term  w(k),  and  the 

spectral  dissipation  term  d(k).  The  net  spectral  flux  is  computed  using  the  collision 
integral  relating  the  component  wavenumbers  (Hassehnann,  1962, 1963) 

•  T(k)  *  P*JJJ  Q^{[$(l*)+ S(fci)]S(ka)S(ks)  -[Sflsj) + S(k3)]S(k)S(k,)} 

x80t  k]  -  k2  -  k})8(a> + 00]  -  a>2  -  (Dsldkj  dka  dk3 

(1.12) 
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where  Sit  the  action  spectral  density  S(k)  -  (g/tt^fk)  and  Q  is  a  complicated 
function  of  order  k*.  This  expression  is  dourly  too  complex  to  evaluate  analytically 
and  la  computed  numerically  when  used  by  wave  evolution  models  (WAMDIG, 
1988).  Sefl  4c  Hasactmtnn  (1972),  Longuet-Higgins  C 1978 )  aid  Fox  (1976)  have 
developed  a  narrow  bandwidth  approximation  to  the  net  spectral  flux  tom.  They 
found  that  die  interactions  are  primarily  local  for  die  wind  wave  components  whose 
amplitudes  and  slopes  are  small  Phillips  (1985)  then  argued  that  Equation  1.12 
behaves  as 


-Vk  TW-pQ^OOkVco 

-pgk^Ot) 


(1.13) 


By  combining  measurements  of  atmospheric  pressure  fluctuations  over  surface  waves 
in  die  Bight  of  Abaco  and  Miles'  (1952)  wind  input  model,  Snyder,  Dobson,  Elliott  & 
Long  (1981)  proposed  that 


w(k)  «  0.25pa  gmax£o,^“^cos8  - 1  jjco'F(k) 


(114) 


where  pa  is  the  density  of  air,  U5  is  the  wind  speed  5  meters  above  the  mean  sea  level, 
c  is  die  phase  speed  of  die  surface  wave  frequency  component,  and  0  is  die  angle 
between  die  wind  and  wave  direction. 


A  comparable  expression  was  preposed  by  Plant  (1982)  from  a  survey  of  wind  input 
measurements  of  Kawai  (1979),  Snyder  et  al  (1981),  Shemdin  &  Hsu  (1967),  Larson  & 
Wright  (1975),  Plant  4c  Wright  (1977),  Wu,  Hsu  4c  Street  (1979),  and  Stewart  4c  Teague 
(1980).  He  showed  that  the  various  wind  input  measurements  satisfy  the  equation 
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w(k)  -  (a04  ±0.02)^  jcos6^  *  j  «*(k) 

where  u»  is  the  friction  velocity  of  the  air  flow  over  the  waves  (Figure  1.9). 

White  the  characteristic*  of  both  the  spectral  flux  and  wind  input  teems  are  clearly 
complex,  it  is  generally  believed  that  the  expression  for  the  dissipation  d  is  the  least 
understood  component  of  the  energy  transfer  equation.  Because  of  the  difficulty  in 
quantifying  this  in  the  field,  there  have  been  no  direct  field  measurements  of  the  total 
surface  wave  dissipation.  Only  recently  have  investigators  begun  to  make  field 
measurements  of  fee  volumetric  dissipation  rate  in  the  water  column  adjacent  to  the 
air-water  interface. 

Agrawal  et  al  (1992)  measured  die  surface  velocity  field  under  conditions  of  strong 
breaking  from  a  tower  in  Lake  Ontario.  They  found  that  die  rate  of  dissipation  of 
kinetic  energy  under  those  conditions  can  exceed  the  values  predicted  by  assuming  a 
velocity  distribution  based  on  a  fiat,  constant  stress,  turbulent  boundary  layer  or  'wall- 
layer'  theory  (Tennekes  4c  Lumfey,  1972)  by  two  orders  of  magnitude.  They  concluded 
that  die  near-surface  layer  is  a  region  of  enhanced  dissipation  and  that  studies  of 
upper  mixed  layer  dynamics  have  to  incorporate  the  influence  of  wave  breaking 
explicitly. 


14  Breaking  wavs  dissipation  models. 


In  principle,  the  total  wave  energy  dissipation  can  be  computed  from  the  sum  of  the 
energy  dissipated  by  all  file  breaking  events  for  an  area  of  the  sea  surface  over  a  given 
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Figure  1.9.  A  plot  of  the  wind  induced  wave  growth  w/'PCk)  as  a  function  of 
frequency  at  u* »  0.45  m/s.  Data  were  compiled  from  seven  different  investigations 
by  Plant  (1982).  Dashed  lines  show  limits  of  Equation  1.14.  Data  from  Snyder  etal 
(1981)  is  indicated  by  shaded  area  (Plant,  1982). 
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length  of  time.  The  difficulty  in  this  method  Is  in  coming  up  with  a  technique  that  can 
consistently  identify  all  the  breaking  events  over  the  given  area  and  in  determining 
how  much  energy  is  dissipated  by  each  of  the  observed  events.  However,  this 
approach  can  be  approximated  by  assuming  that  the  distribution  of  breaking  waves  in 
the  ocean  is  statistically  stationary,  and  that  the  wind  and  wave  conditions  are  quasi¬ 
steady  so  that  the  percentage  of  wave  crests  that  break  at  a  given  location  over  a  finite 
period  of  time  is  proportional  to  the  percentage  of  breaking  wave  crests  over  the  sea 
surface  at  any  given  time. 

fat  a  recent  paper,  Thorpe  (1992)  suggested  that  die  fraction  of  breaking  wave  crests  £ 
is  related  to  the  inverse  wave  age 

\3 

(1.16) 


5«(4±2)xl(T^ 


( 


This  equation  was  obtained  front  an  empirical  fit  to  die  data  from  lake  measurements 
of  Thorpe  6c  Humphries  (1980),  Weissman,  Atakturk  &  Katsaros  (1984),  Atakturk  & 
Katsaros  (1991)  and  Thorpe  (1990).  We  see  from  Figure  1.10  that  this  compilation  also 
includes  two  points  from  Longuet-Higgins  6c  Smith's  (1983)  tower  measurements  off 
the  coast  of  Netherlands  whose  data  actually  fall  above  the  upper  limit  of  Equation 
1.16.  There  is  also  enough  scatter  in  the  data  points  that  they  look  less  correlated 
without  the  dashed  lines. 

fat  any  case,  he  assumed  that  on  average,  die  energy  dissipated  per  unit  length  of  a 
breaker  can  be  estimated  from  laboratory  measurements  of  dissipation  by  steady 
breaking  waves  (Duncan,  1981) 
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Figure  1.10.  Log-log  plot  of  die  fraction  of  breaking  waves  £  versus  the  inverse  wave 
age  (Uio/Cp).  Data  was  taken  from  various  field  investigations:  Longuet-Higgins  & 
Smith  (1983)  -  ■,  Weissman  et  al  (1984)  WA+K,  Thorpe  &  Humphries  (1980)  -  O, 
Katsaros  &  Atakturk  (1991)  -  •,  and  Thorpe  (1990)  -  □  (Thorpe,  1992). 
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d,*(0.044±0.008)Aw4/g  (117> 

where  q,  is  the  speed  of  propagation  of  the  wave  crest  Therefore,  the  total  energy 
dissipated  per  unit  surface  area  of  the  ocean  is 

Dr-«U/*p 

.(W±U)X10-*^.M  (1.18) 

where  Ap  is  the  wavelength  of  the  dominant  wave.  If  deep  water  linear  dispersion  is 
assumed,  cp2  =  gAp/2*  and 


Dr 


*(3.0±L8)xlO“5pwU|o 


(119) 


We  can,  in  principle,  measure  the  variables  in  Equation  1.19  to  compute  Dp.  However, 
unlike  the  other  variables  an  the  left  hand  side  of  that  equation,  q,  is  difficult  to 
measure.  Although  Ding  &  Farmer  (1992)  have  recently  developed  a  method  for 
tracking  the  motion  of  breaking  waves  using  the  phase  differences  in  the  acoustic 
signal  received  by  an  array  of  hydrophones,  the  value  of  q,  is  difficult  obtain  using 
commonly  available  meteorological  and  oceanographic  instruments. 

Alternatively,  Thorpe  (1992)  suggested  that  if  q,  is  a  fraction  of  Cp  (or  if  the  wave 
frequency  components  higher  than  the  spectral  peak  break) ,  we  can  obtain  q,  as  a 
function  of  Cp  if  we  know  the  total  energy  dissipation  Dp  in  Equation  1.19.  He 
assumed  that  the  volumetric  dissipation  in  the  upper  5  meters  is  comparable  to  Oakey 
&  Elliot's  (1982)  measurements  of  the  mean  volumetric  dissipation  in  the  20-meter 
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upper  mixed  layer  5-10  meters  below  the  sea  surface.  Thorpe  (1992)  found  that  if  he 
assumed  that  (c^/Cp)5  *  1  in  Equation  2.30,  then  Dp  exceeds  Oakey  k  Elliot's  (1982) 
dissipation  estimates  by  a  factor  of  1000.  Hence,  for  die  two  to  be  comparable, 

(Cb/Cp)5  - 10-3 

Cb-0.25cp.  (1.20) 

The  measurements  of  Agrawal  et  al  (1992)  showed  that  the  dissipation  just  beneath 
die  water  surface  to  a  depth  of  die  order  of  die  waveheight  can  be  up  to  two  orders  of 
magnitude  larger  than  that  predicted  by  a  turbulent  boundary  layer  near  a  rigid  wall 
Melville  (1992)  proposed  these  results  suggest  die  existence  of  a  shallow  layer  of 
enhanced  dissipation  just  beneath  the  sea  surface.  He  then  showed  that  by  combining 
Phillips'  (1985)  dissipation  model  and  measurements  of  turbulence  and  mixing  by 
laboratory  breaking  waves  (Rapp  k  Melville,  1990),  the  rate  of  enhanced  dissipation 
and  die  thickness  of  this  layer  can  be  estimated,  and  that  these  estimates  are 
consistent  with  Agrawal  et  al's  (1992)  results. 

In  comparing  his  model  with  that  of  Thorpe  (1992),  Melville  (1992)  argued  that  by 
using  Duncan's  (1981)  dissipation  measurements  of  steady  laboratory  breaking 
waves,  Thorpe's  (1992)  dissipation  formula  Dp  (Equation  1.19)  overestimates  the  total 
energy  dissipated  by  individual  breaking  events.  He  showed  that  laboratory 
measurements  (Loewen  k  Melville,  1991a)  suggest  that  the  energy  dissipated  by 
unsteady  breaking  waves  is  given  by 


du=(0.0032-0.016)pwcj/g. 


(1.21) 
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He  suggested  diet  since  spiffing  waves  constitute  most  of  the  breakers  in  deep  water, 
the  value  of  d^  that  should  be  used  should  be  closer  to  die  lower  limit  of  Equation 
1.21,  approximately  one  order  of  magnitude  smaller  than  Equation  1.17.  Hence, 
Melville  (1992)  concluded  that  Dj>  (Equation  1.19)  overestimates  the  total  dissipation 
of  individual  breaking  waves  by  one  order  of  magnitude. 

Melville  (1992)  also  disagreed  with  Thorpe  (1992)  that  die  total  dissipation  should  be 
comparable  to  die  upper  mixed  layer  dissipation  measured  by  Oakey  &  Elliott  (1982). 
Based  on  his  estimates  of  enhanced  dissipation  in  the  shallow  layer  beneath  die  sea 
surface,  Melville  (1992)  suggested  that  die  dissipation  rate  is  two  orders  of  magnitude 
larger  than  Oakey  &  Elliot's  (1982)  measurements  over  a  depth  one  order  of 
magnitude  smaller  than  die  upper  mixed  layer. 

Melville  (1992)  proposed  that  since  Equation  1.19  overestimates  the  total  individual 
breaking  wave  dissipation  by  a  factor  of  10,  and  since  using  Oakey  &  Elliot’s  (1982) 
results  underestimates  the  surface  wave  dissipation  by  a  factor  of  10,  then  the  ratio 
(Cb/Cp)5  in  Equation  1.20  is  should  be  0.1.  Consequently,  Melville  (1992)  concluded 
that 


Cb-0.63cp.  (1.22) 

This  value  is  comparable  to  the  results  of  Ding  (1993)  who  found  that  the  speed  of 
breaking  waves  they  observed  is  0.4-0.7cp  (Figure  1.11). 

Both  Thorpe  (1992)  and  Melville  (1992)  propose  methods  for  estimating  the  total 
dissipation  from  measurements  of  die  total  dissipation  by  breaking  waves.  However, 
their  approach  suffers  from  two  major  drawbacks.  First,  the  accuracy  of  estimates 
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Figure  1.11.  Normalized  speed  of  a  breaking  event  (q,/cp)  versus  the  phase  speed  cp 
(Ding,  1993). 
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using  their  models  is  premised  on  the  validity  of  the  empirical  relationship  between  \ 
and  (U10/Cp)  given  by  Equation  1.16.  Given  the  scatter  in  Figure  1.10,  there  is  a  need 
to  obtain  more  breaking  wave  surveys  to  adequately  justify  that  Equation  1.16 
accurately  describes  this  relationship.  Second,  since  the  ratio  (q,/Cp)  is  raised  to  the 
5th  power,  then  small  errors  in  the  estimate  of  this  value  will  result  in  large 
differences  in  die  total  dissipation  estimate.  For  example,  die  upper  and  lower  limits 
of  Ding's  (1993)  measurements  of  q,  gives  dissipation  estimates  that  differ  by  one 
order  of  magnitude. 


1.7  Outline  of  thesis. 


The  results  of  the  laboratory  experiments  of  Melville  and  his  coworkers  suggest  that 
die  ambient  noise  level  in  die  ocean  is  related  to  die  magnitude  of  the  wave  field 
characteristics  that  correlate  with  the  frequency  and  die  strength  of  breaking.  In 
particular,  die  laboratory  experiments  suggest  that  wave  energy  dissipation  and 
ambient  noise  generation  are  coupled  together  by  the  air  entrainment  process  in 
breaking  waves.  While  the  dissipation  of  surface  wave  energy  can  be  quantified  in  the 
laboratory  under  well-controlled  conditions,  as  the  lake  experiment  by  Agrawal  et  al 
(1992)  indicates,  this  cannot  be  easily  accomplished  in  the  field.  Hence,  we  need  to 
resort  to  models  to  provide  us  with  estimates  of  dissipation  in  the  ocean.  In  Chapter  2, 
we  will  discuss  the  two  major  dissipation  models  in  the  literature.  We  will  then  derive 
methods  of  estimating  dissipation  from  integral  characteristics  of  the  wave  field  based 
on  the  assumptions  of  the  two  models.  In  Chapter  3,  we  will  describe  an  experiment 
we  conducted  cm  a  research  platform  off  the  coast  of  Oregon  where  we  made 
simultaneous  measurements  of  ambient  noise,  wind  speed  and  wind  direction, 
directional  wave  spectra,  and  video  recordings  of  the  occurrence  of  breaking  waves. 

In  Chapter  4,  we  will  discuss  the  results  of  the  measurements  and  examine  the 
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characteristics  of  the  wave  spectrum  and  how  ambient  noise  relates  to  different  wave 
field  parameters.  In  Chapter  5,  we  will  use  die  wave  measurements  to  compute 
estimates  of  the  dissipation  using  the  dissipation  models  described  in  Chapter  2.  We 
wifi  then  show  that  die  relationship  between  wave  energy  dissipation  and  ambient 
noise  can  also  be  described  by  well-correlated  power  law  expressions.  These  results 
are  then  summarized  in  Chapter  6  where  the  consequences  of  these  results  for  the 
remote  sensing  of  the  upper  ocean  using  ambient  noise  are  discussed. 
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2  SURFACE  WAVE  DISSIPATION  MODELS 


Since  reliable  methods  far  computing  dissipation  directly  from  breaking  waves  in  the 
field  are  still  currently  unavailable,  we  adopt  the  alternative  approach  of  estimating 
dissipation  from  die  energy  budget  of  a  wave  evolution  model  This  method  was 
introduced  in  Section  1.5  and  is  discussed  in  detail  here  In  this  chapter,  we  examine 
how  the  assumption  of  equilibrium  of  die  wind  input,  dissipation  and  nonlinear 
spectral  flux  can  lead  to  expressions  for  the  spectral  dissipation.  We  then  show  that 
die  total  dissipation  can  be  estimated  by  substituting  simple  idealized  forms  of  die 
wave  spectrum  into  the  spectral  dissipation  equation. 

2-1  Equilibrium  spectral  slopes. 

In  an  earlier  paper,  Phillips  (1958)  argued  that  in  a  fully  developed  sea,  die 
characteristics  of  die  wind  generated  waves  are  described  by  an  equilibrium 
determined  by  wave  breaking  but  independent  of  time  and  fetch.  Using  dimensional 
arguments,  he  concluded  that  the  equilibrium  wave  height  frequency  spectrum  is 

®(co)«ing2cif3,  (2.1) 


where  m  is  a  non-dimensional  constant  and 


•(«) «  ~  Ja(*.  t)a(x,  t + x)t~im6x . 


(2.2) 


He  also  concluded  that  die  wavenumber  spectrum  is 
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¥<k)~F(0)lf 


(23) 


where  F(0)  if  the  directional  wavenumber  distribution.  Observational  evidence  from 
subsequent  experiments  |e.g.,  Toba  (1973),  Kawai,  Okuda  &  Toba  (1977),  Forristal 
(1981),  Donelan,  Hamilton  k  Hui  (1984)]  suggest  that  the  form  of  the  wave  spectrum 
in  fact  depends  on  fetch.  These  investigations  show  that  the  form  of  the  wind  wave 
spectrum  is  closer  to  Toba's  (1973)  expression 


«&(a>)  *  ccmgaT4,  (2.4) 

which  he  obtained  from  wind-wave  tunnel  experiments  and  dimensional  analysis. 

In  toe  last  chapter,  we  discussed  how  toe  energy  transfer  equation  (Equation  1.11) 
describes  toe  evolution  of  the  surface  wave  field.  Phillips  (1985)  argued  that  lor  the 
higher  frequency  components  of  toe  wave  field  (wind  waves),  die  temporal  and 
spatial  scales  for  growth  and  decay  are  significantly  longer  than  the  internal  wind 
wave  time  and  length  scales.  From  this,  he  assumed  that  these  waves  are  essentially 
in  equilibrium 


Pwg-“*-Vk-  T(k) + W(k)  -  d(k) s °,  (2.5) 

w 

and  that  the  three  source  terms  are  proportional  and  of  comparable  order 


-VkT(k)~w(k)~d(k). 


(2.6) 
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Using  ftae  assumptions  and  combining  the  spectral  flux  form  (Equation  1.13)  and 
the  Plant  (1982)  wind  input  expression  (Equation  1.15),  he  showed  that  die 
equilibrium  forma  for  die  wavenumber  and  wave  frequency  spectra  are 

^bOO  -  u*"1/V7/2  (2.7) 

^(aO-amgm-4  (2.8) 

where  a  -0.11  (the  Toba  coefficient).  The  value  of  fi  is  obtained  from  die  equation 
P  -  a/4I(p)  which  relates  die  magnitudes  of  the  frequency  and  the  wavenumber 
spectra,  and  die  surface  wave  spreading  function 

f«/2 

I(p)«  I  cos'ede.  (2.9) 

J- un 


For  p  «  0.5, 1(p) «  2.4  and  P  =  0.0115. 


12  Phillis  diacipation  modeL 


Using  the  equilibrium  arguments  described  in  the  last  section,  Phillips  (1985) 
evaluated  d(k)  from  Equation  2.5  using  the  proposed  forms  of  die  expressions  for  die 
spectral  flux  Vk  T(k)  and  w(k).  He  showed  that 

W-PwTJ^fk4^3^)]3  (2.10) 

d(k)  =  PwYgk,,F3(k)  (111) 
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dtgffttf  2 

where  y  is  a  numerical  constant  Since  w(k)  and  d(k)  are  of  comparable  magnitude, 
the  value  of  7  can  be  estimated  by  equating  dswapation  with  the  wind  input  term, 

OP-  -  0.04.  Hence,  7  "305.  To  express  d(k)  in  terms  of  $  and  co,  we  substitute  Equation 
2.7  into  Equation  2.10,  use  the  linear  dispersion  relation  co2  «  gk,  and  get 

d-pwfTk~*[poot'eju«"1/2k1/2]3 

'p*‘Yk"fe7aa4“"I 

<Km)-p»nfk-<^^^«>E(co)l.  (2.12) 

The  total  Phillips  dissipation  in  the  equilibrium  range  is 


To  obtain  an  estimate  of  the  total  dissipation  based  an  Phillips  (1985)  assumptions 
using  tiie  measured  wave  spectrum,  we  use  0(<d)  instead  of  Qe(go)  and  substitute 
Equation  2.12  into  Equation  2.13,  we  get  an  expression  for  the  dissipation  in  the  wind 
wave  range  from  field  measurements  of  the  wave  spectrum  based  on  the  Phillips 
(1985)  assumptions 


Dp  =  IQ&Bxi  X  f  co1I<I>3(co)dco. 

KStlCrtf  K 


(2.M) 
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The  lower  Mmit  of  the  equilibrium  range  is  assumed  to  be  the  spectral  wind  wave 
peak.  The  integral  m  Equation  2.14  will  converge  as  long  as the  spectral  dope  erf  O(co) 
it  steeper  than  or4  when  cd  -t  Note  that  our  use  of  $(a>)  instead  of  Ogfoo)  is  a 
modification  of  Phillips'  (1965)  explicit  assumption  of  an  equilibrium  spectrum. 
However,  Equation  2. 14  is  justified  if  it  can  be  demonstrated  that  $(<»)  it  reasonably 
approximated  by  Phillips  (1985)  assumed  equilibrium  form  of  fee  wave  frequency 
spectrum  (Equation  18). 

If  we  instead  combine  Equations  2.7, 2.11  and  2.13,  we  get  Phillips'  (1985)  reduced 

-S» — *  —  -**  —  — ^ «  - 

TiWynriii  wyiiiMy* 


(2.15) 


«2yP3l(3p)pwu? 


(2.1 0 


Phiffips'  (1965)  assumptions  lead  to  a  logarithmic  increase  in  the  dissipation  with 
increasing  wavenumber.  Hence,  there  is  a  need  to  specify  an  upper  spectral  limit  to 
his  dissipation  expression.  Phillips  (1985)  assumed  that  die  lower  limit  of  the 
equilibrium  range  to  be  die  spectral  peak,  ko  *  g  /  .  He  then  argued  that  the  upper 

limit  is  determined  by  die  presence  of  wind  drift  c^,  which  suppresses  die  formation 
of  high  frequency  waves  traveling  at  the  same  velocity  (Banner  &  Phillips,  1984). 
Using  die  work  of  Keulegan  (1951)  whose  field  measurements  suggest  that  *  u»,  he 
chose  die  upper  limit  to  be  kj  *=rg/u?  where  r  is  a  constant  of  order  1. 
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-2*>I (3p)p.o.s 


(2.17) 


Hsssehninn  (1974)  argued  that  although  wave  breaking  is  locally  a  strong  nonlinear 
process,  it  it  in  general  *wcak  kt  the  mean1  and  its  effect  on  the  dissipation  d(co)  ia  a 
quasHmef  function  of  $(<n)  and  a  damping  coefficient  proportional  to  die  square  of 
the  frequency  a*,  Le, 


<Ka>-iK»fa(<0).  (2.18) 

He  suggested  that  the  value  of  the  coefficient  t)  be  obtained  from  an  equilibrium 
balance  of  the  three  source  terms  on  die  right  hand  side  of  die  energy  transfer 
equation.  Komen,  Hanehnaxm  Ac  Hassebnaim  (1984)  [herein  referred  to  as  KHH] 
propoaed  an  expression  for  the  coefficient  q  which  they  obtained  by  performing 
numerical  simulations  on  the  forms  erf  tj  to  reproduce  the  characteristics  of  fetch- 
limited  wave  growth.  They  used  a  numerical  evaluation  of  die  full  Hassehnann  (1962, 
1963)  spectral  flux  equation  (Equation  1.12)  for  Vk  TOO  and  the  Snyder  et  al  (1981) 
wind  input  model  for  w(k).  A  paper  by  die  Wave  Model  Development  and 
Implementation  Group  (WAMDIG,  1988)  describes  a  numerical  implementation  of  a 
wave  development  model  using  a  numerically  stable  form  of  the  KHH  dissipation 
expression.  While  the  WAMDIG  (1988)  expression  is  slightly  different  from  die  KHH 
formula,  it  is  based  on  the  same  assumptions  and  gives  dissipation  estimates  that  are 


CkapiarZ 

of  the  same  order  of  magnitude  as  those  given  by  the  KHH  formula.  From  KHH,  the 
total  dissipation  due  to  wave  breaking  is 

I>H«Pwt  fc^CO/SfCo/^K^S^CDldtt  (2.19) 


where  Cff  •  333  x  10^,  5  ■  ec^/g2  is  a  measure  of  the  wave  steepness, 

SpM -4.57x10-3  is  the  value  of  a  fora  fuUy  developed  wind  sea  based  on  the  Pierson- 
Moskowitz  spectrum,  pw  « 1025  kg/m3  is  the  density  of  sea  water,  e*J^O(co)da)  is 

the  total  surface  wave  energy,  and 

<o  *  J^4a>dOD/£,Qd(D  (2.20) 

is  the  mean  frequency  of  the  wave  spectrum. 

Although  both  KHH  ami  Phillips  (1985)  derived  dissipation  expressions  by  balancing 
the  contributions  of  the  three  source  terms,  they  arrived  at  substantially  different 
expressions  for  D  because  of  the  differences  in  their  assumed  spectral  dependence  of 
d(a>),  differences  in  die  source  functions  for  w(k)  and  Vk  T(k)  used,  and  differences 
in  their  metivxi  of  balancing  the  source  terms. 

2.4  Reduced  Hassehnann  dissipation  estimate. 

In  many  instances,  the  full  wave  spectrum  is  not  available.  Instead,  the  summary 
wave  data  is  limited  to  a  few  variables,  e.g.,  root  mean  square  wave  amplitude  a  and 
the  mean  frequency  of  the  wave  spectrum  F.  We  can  use  these  parameters  to  obtain  a 
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dissipation  estimate  by  assuming  that  the  frequency  spectra  can  be  approximated  by  a 
simple  power  law  in  terms  of  to.  This  was  essentially  Phillips'  (1985)  approach  when 
he  derived  Equation  2.16  from  Equation  2.11.  However,  his  rationale  was  based 
primarily  on  the  assumption  of  the  existence  of  equilibrium  spectra  and  is  valid  only 
inthewind  wave  region  of  the  spectrum. 

While  die  basis  of  die  KHH  dissipation  expression  is  independent  of  die  existence  and 
form  of  an  equilibrium  spectrum,  based  on  the  discussions  in  Section  2. 1 ,  we 
nevertheless  assume  an  idealized  wave  spectrum  of  die  form 


0*4(00)=  miug  to-4 

*0 


OOiOOp 

a>«0p 


(2.21) 


where  m  is  a  constant  Both  the  frequency  and  wavenumber  spectrum  satisfy 
Paraeval's  theorem 


a*  *  j*  'F(k)dkIdka  =  Jo(co)da>. 


(2.22) 


Substituting  Equation  2.21  into  Equation  2.22  and  integrating  gives 


a2*  jmiugcOp3 


m 


u*8 


(123) 


Combining  Equations  121  and  123  gives 


$*(<»)  *3fl*a£pf4. 


00  2  (Dp. 


(2.24) 


The  relationship  between  oo  and  00p  is  obtained  by  substituting  Equation  2.24  into  the 
definition  of  a>  (Equation  2.20).  For  a  spectrum  with  an  (O'4  tail 

co  =  1.5tOp.  (2.25) 


Substituting  Equations  2.24  and  2.25  into  Equation  2.19  gives  the  equation  for  the 
reduced  KHH  dissipation 


*0.889pwgcw(a/o^M)2fl2w2  I  a>“2 doo. 


(2.26) 


Integrating  Equation  2.26  gives 


Eh 


« 1.333pwg  Cft 


cfa. 


(2.27) 


In  this  chapter,  we  derived  four  spectral  estimates  of  the  dissipation  based  on  the 
wave  evolution  models.  To  evaluate  the  total  dissipation  based  on  these  models,  we 
need  measurements  of  both  the  wave  spectrum  and  the  wind  speed.  In  the  next 
chapter,  we  will  describe  an  experiment  whose  measurements  will  be  used  to 
compute  the  surface  wave  energy  dissipation  using  these  models.  In  Chapter  5,  we 
will  then  examine  the  spectral  characteristics  of  these  estimates  and  compare  these 
results  with  ambient  noise  measurements. 
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9  EXPERIMENT 


In  the  literature  review,  we  saw  that  although  there  have  been  numerous  studies  that 
dearly  show  that  there  is  a  well-correlated  relationship  between  wind  speed  and  the 
prevailing  underwater  ambient  noise  in  the  Knudsen  range,  the  mechanics  behind 
this  relationship  is  not  well  understood.  It  is  believed  that  the  sound  generated  by 
breaking  waves  are  responsible  to  the  wind  dependent  noise  We  pointed  out  that 
although  laboratory  experiments  have  shown  that  wave  field  parameters  that 
characterise  wave  breaking  correlate  with  the  acoustic  energy  radiated,  there  have 
been  no  measurements  that  show  that  this  relationship  is  true  in  the  field  as  well. 

The  main  objectives  of  our  experiment  are:  to  correlate  the  prevailing  noise  spectrum 
level  with  wave  field  parameters,  particularly  wave  height,  wave  slope  and  wave 
spectra;  to  relate  these  results  to  measurements  of  the  wind  speed  dependence  of 
ambient  noise;  to  examine  die  spectral  characteristics  of  surface  wave  dissipation 
models  and  to  study  how  dissipation  correlates  with  ambient  noise. 

3.1  NOBS  Description. 

The  main  thrust  of  the  NOBS  (Noise  on  Basalts  and  Sediments)  experiment  was  to 
study  the  influence  of  the  surface  wave  spectrum  on  the  very  low  frequency  [0(1  Hz)] 
ambient  noise  on  the  sea  floor.  The  principal  investigator  was  John  Hildebrand  of  the 
Soipps  Institution  of  Oceanography  (SIO)  of  the  University  of  California  San  Diego. 
The  Upper  Ocean  Physics  Group  of  Robert  Pinkel  (SIO)  and  Jerome  Smith  (SIO) 
conducted  meteorological  measurements,  as  well  as  measurements  of  the  surface 
wave  directional  spectra  using  an  acoustic  Doppler  sonar  system  (Smith,  1989)  and  a 
wire  wave  gauge  array. 
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TabU  3.1.  Instruments  used  In  Ail  experiment 


Instrument 

Institution 

Sample 
rate  (Hz) 

Type 

Directional 

hydrophone 

MIT 

51,200 

Fabricated  at  MIT 
using  ITC-8181A 

Anemometer 

SIO 

8 

Weathertronks 

Model  2020/2030 

Air  temperature 

SIO 

8 

wefttnertronKs 

Model  4470 

Barometer 

SIO 

8 

Weathertronks 

Model  7105-A 

Accelerometer 

SIO 

8 

Fabricated  at  SIO 

Heading 

Flip 

8 

Flip  instrumentation 

Capacitance  wire 
wavegauge 

MIT 

8 

Fabricated  at  MIT 

4-wire  resistance 
wavegauge  array 

SIO 

8 

Fabricated  at  SIO 

Video  camera 

MIT 

N.A. 

NECTI-23A 

Although  our  experiment  was  unrelated  to  the  main  objective  of  NOBS,  observations 
of  the  wave  field  and  the  meteorological  parameters  were  common  to  both  the  ocean 
bottom  studies  of  Hildebrand  and  to  ours.  Table  3.1  summarizes  the  major 
instruments  used  in  our  study. 


The  research  platform  Flip,  a  manned  spar  buoy,  was  towed  from  San  Diego, 
California  and  moored  in  3000  meters  of  water  approximately  130  km  off  the  coast  of 


Figure  3.1.  A  map  showing  the  location  of  Flip  during  the  NOBS  experiment  (43°  42* 
N,  125°  59'  W).  The  site  is  approximately  130  km  west  of  Oregon  and  600  km  south  of 
Vancouver  Island. 
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Figure  33.  Schematic  diagram  of  the  sampling  and  preprocessing  method  used  in 
gathering  and  storing  data. 
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Reedsport,  Oregon  in  September,  1991  (Figure  3.1).  Flip  was  deployed  in  a  3-pt 
mooring  configuration.  The  length  from  the  bow  to  the  stem  of  Flip  is  approximately 
110  meters  and  the  cylindrical  hull  has  a  maximum  diameter  of  6  meters.  At  vertical 
attitude,  90  meters  of  Flip  is  submerged  under  water.  The  section  of  the  hull  which 
intersect*  the  mean  water  tine  it  approximately  3.8  meter*  in  diameter  (Bronson  fc 
Qoaten,  1965).  The  surface  wave*  reflected  from  die  hull  of  the  platform  could  be 
observed  by  eye  but  these  did  not  appear  to  have  a  noticeable  effect  on  the  wave  field 
near  die  capacitance  and  resistance  wavegauges  15  and  17  m  away.  Figure  3.2  shows 
a  diagram  of  Flip  and  die  locations  of  the  different  instruments  used.  A  schematic 
diagram  of  die  data  sampling  and  preprocessing  method  is  shown  in  Figure  33.  This 
is  described  in  detail  in  the  subsequent  sections. 

3.2  Directional  hydrophone. 

The  hydrophone  was  mounted  an  Flip  33  meters  below  sea  leveL  Since  one  erf  die 
initial  objectives  of  this  experiment  was  to  track  breaking  waves  using  the  beam 
pattern  of  a  directional  hydrophone,  and  because  there  was  some  concern  that  wave 
splashing  on  the  hull  of  Flip  might  dominate  the  acoustic  signal,  it  was  considered 
desirable  to  introduce  directionality  to  die  hydrophone.  This  was  accomplished  by 
mounting  an  ITC-8181A  omnidirectional  hydrophone  on  die  focus  of  a  42-inch 
parabolic  dish  (Figures  3.4, 3.5  and  3.6).  The  parabola  was  filled  with  foam  to  provide 
a  perfectly  reflecting  surface  and  improve  the  directional  characteristics  of  die 
hydrophone  The  resulting  beam  pattern  was  obtained  from  a  calibration  of  die 
directional  hydrophone  at  the  Naval  Ocean  Systems  Center  (NOSC)  Transducer 
Calibration  Facility  in  San  Diego,  CA.  The  results  of  die  NOSC  calibration  are  given  in 
Appendix  A.1. 


Figure  3.4.  A  photograph  of  toe  directional  hydrophone  taken  immediately  after  the 
platform  'flipped'  from  vertical  to  horizontal  attitude.  The  photograph  was  taken  from 
a  location  near  the  bow  of  Flip.  The  hydrophone  is  the  dark  circular  disk  at  the  center 
of  the  photograph. 


Figure  3.5.  A  dose-up  photograph  of  the  directional  hydrophone  while  Flip  was  in 
port  in  San  Diego.  The  instruments  in  die  foreground  are  die  SIO  Doppler  sonars  used 
to  image  the  sea  surface. 
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The  configuration  resulted  in  *  3-dB  main  lobe  with  a  2P  inducted  angle  at  3  kHz  and 
6*  el  15  kHz.  The  centerline  oi  the  parabola  was  oriented  upward  at  a  45*  angle  from 
the  kvtgitodinal  axis  of  Ffy.  At  3  kHz,  the  main  lobe  creates  an  elliptical  footprint  an 
the  gee  surface  whose  major  and  minor  axes  are  30  meters  and  15  meten  long.  At  25 
kHz,  the  major  and  minor  axes  of  the  effiptical  footprint  are  33  meters  and  7  meters 
long. 

A  video  camera  mounted  on  the  crow's  nett  28  m  above  see  level  was  oriented 
towards  the  footprint  of  the  hydrophone  mein  lobe.  At  that  location,  the  camera  can 
view  a  5.2  m  x  48  m  section  of  the  sea  surface.  Four-hour  videos  of  the  ocean  surface 
were  taken  every  day  starting  at  2  PM  Pacific  Daylight  Time  (PUT)  or  2100  UT.  The 
video  image  was  embedded  with  the  time  code  to  synchronize  the  video  data  with  the 
acoustic  record.  The  video  data  was  then  recorded  using  a  Panasonic  AG-2500  video 
cassette  recorder. 

321  Dmt*  »mmpUmg  and  preprocessing 

Figure  33  shown  a  schematic  diagr*  of  the  hydrophone  data  conditioning  and 
sampling  method  used.  The  hydro  < >^ie  signal  was  high  passed  at  23  kHz  using  a 
Krohn-Hite  3202R  filter,  amplified  50  dB  using  a  Wilcoxon  Research  AM-5  Low  Noise 
Amplifier,  and  low-pass  filtered  at  27  kHz,  sampled  &  preprocessed  using  die 
procedure  described  in  die  next  paragraph  using  an  ALR  486/25  PC  equipped  with  a 
Spectrum  Signal  Processing  TMS320C30  Real  Time  Board.  This  add-in  board  is  an 
independent  computer  that  has  a  specialized  digital  signal  processing  unit.  It 
communicates  with  die  host  PC  through  the  16-bit  ISA  bus.  The  data  was  stored  on 
650  MB  Sony  EDM-1DA1  magneto-optical  disks.  Selected  ambient  noise  data  were 
recorded  in  DAT  cassettes  using  a  Panasonic  SV-255  DAT  recorder. 


Prior  to  the  start  of  each  10-minute  sampling  segment,  the  tune  code  generated  by  a 
Datum  9900  Time  Code  Generator  was  sampled  from  the  PC  aerial  port  and  recorded. 
The  hydrophone  data  was  then  sampled  at  51.2  kHz  using  the  TMS320C30  Real  Time 
Board.  The  power  spectrum  was  computed  and  avenged  from  nine  512-point  FFTs 
and  then  tent  to  the  RAM  of  the  host  PC  every  0.1  seconds  (approximately).  This 
method  generated  6000  256-point  real  power  spectra  during  each  10-minute  sampling 
segment  which  were  stored  in  2-byte  integer  records.  The  time  code  was  sampled  and 
recorded  at  die  end  of  the  segment  prior  to  storing  die  data  in  the  optical  drive. 

3JL2  Data  procuring. 

The  ambient  noise  spectrum  level  N(f)  [dB  re  1  uPaVl  Hz]  for  an  equivalent 
omnidirectional  hydrophone  was  computed  from  die  raw  power  spectrum  R(f)  using 
the  formula 


N(f)-R(f)-SA+SL(f)+AG(f)  (3.1) 

where  R  is  in  dB  re  1  Volt2/!  Hz,  SA  is  the  signal  amplification  (50  dB),  and  SL  is  the 
hydrophone  sensitivity  level  in  dB  re  1  pPaVl  Volt2  (Figure  3.7).  The  array  gain  AG 
in  dB  was  computed  using  die  formula  (Dyer,  1989) 


AG  *10log 


(3.2) 
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where  Q  is  the  solid  angle.  In  computing  AG,  radial  symmetry  was  assumed  in  both 
the  hydrophone  pattern  B  and  the  sound  source  pattern  S.  A  compilation  of  field 
measurements  of  the  vertical  directionality  of  ambient  noise  (Urick,  1966)  suggest  that 
the  sound  source  pattern  can  be  modeled  as  S(+)  -  cos2  4  where  +  is  the  angle  from 
the  vertical.  The  details  on  how  AG  was  computed  is  given  m  Appendix  A.1.  Figure 
3.8  shows  a  plot  of  AG  computed  from  die  beam  pattern  of  the  hydrophone. 

The  interaction  between  the  hydrophone  and  several  of  the  instruments  on  Flip  led  to 
the  contamination  of  some  of  the  hydrophone  signaL  The  main  noise  sources  were 
characterised  as  follows: 


1.  Unidentified  noise  sources  covering  a  broad  band  of  frequencies  above  14  kHz 
made  acoustic  infbnnatian  in  this  portion  of  die  spectrum  inaccessible  in  a 
consistent  manner.  The  acoustic  data  above  142  kHz  was  considered  to  be 
unusable. 


2.  FUp  operated  an  acoustic  beacon  which  transmitted  a  12-second  4  kHz 
frequency  tone  every  2  min  40  seconds.  Higher  order  harmonics  at  8  and  12  kHz 
of  this  tone  were  also  sampled  by  the  hydrophone  The  data  was  corrected  by 
neglecting  the  power  spectra  from  the  acoustic  data  that  were  sampled  while  the 
sonar  beacon  was  in  operation. 

3.  The  acoustic  Doppler  sonar  mounted  on  the  hull  of  Flip  generated  10- 
millisecond  noise  pulses  every  0.625  seconds.  While  die  pulses  were  at  195  kHz, 
they  produced  Inroad  band  electrical  noise  which  contaminated  acoustic 
frequencies  well  within  our  hydrophones  range  This  has  resulted  in  the 
elimination  of  approximately  one-fifth  of  all  acoustic  power  spectra. 
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Figure  3.7.  Sensitivity  level  (SL)  of  the  directional  hydrophone  in  dB  re  1  VoltVliiPa2 
obtained  from  a  calibration  of  the  instrument  at  die  Naval  Ocean  Systems  Center. 
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Figure  3.8.  Plot  of  array  gain  (AG)  of  the  directional  hydrophone  as  a  function  of 
frequency.  AG  values  computed  from  NOSC  calibration  are  denoted  by  O. 
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Figure  3.9  shows  some  typical  unsmoothed  power  spectra  of  die  ambient  noise  [dB  re 
1  iiPaVl  Hz]  for  wind  speeds  between  4-9  meters/second.  The  spectra  were  averaged 
over  one  hour.  The  broad  peak  at  2.5  kHz  is  due  to  the  high  pass  filter  cutoff  at  that 
frequency. 


The  data  from  environmental  sensors,  the  wave  gauges,  the  compass  and  the 
accelerometers  were  sampled  at  S  Hz  by  a  Macintosh  Ilfx  equipped  with  an  NB-MIO- 
16  analog-digital  converter  (National  Instruments,  Austin,  TX)  (see  Figure  3.3].  The 
data,  which  were  stored  on  optical  disks,  were  provided  to  us  by  Jerome  Smith  of  SIO. 

Hie  wind  speed  U  was  measured  by  a  Weathertronics  (Qualimetrics,  Sacramento, 

CA)  Model  2020  3-cup  anemometer  while  the  wind  direction  was  measured  by  a 
Model  2030  vane.  Both  instruments  were  mounted  on  a  crossbar  positioned  at  28  m 
above  sea  level  These  measurements  were  verified  during  die  experiment  by 
comparing  diem  with  readings  from  die  Flip  anemometer  three  times  per  day.  One 
hour  averages  of  U  were  later  computed  and  reduced  to  die  wind  speed  at  10  meters 
elevation  U10.  In  computing  U1(>  we  used  die  implicit  wave  age  dependent  method 
for  computing  die  friction  velocity  u*  and  die  wind  speed  profile  Uz  described  by 
Maat,  Kraan  &  Oost  (1991).  Details  are  given  in  Section  4.1.1. 

Air  temperature  measurements  were  sampled  using  a  Weathertronics  Model  4470 
Platinum  Resistance  Sensor.  The  barometric  pressure  was  measured  using  a 
Weathertronics  Model  7105-A  analog  barometer.  Although  instruments  for  measuring 
die  water  temperature  and  relative  humidity  were  provided,  these  devices  foiled  to 
function  properly. 
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Figure  3.9.  Some  typical  ambient  noise  power  spectra  in  dB  re  lnPa2/l  Hz  at  various 
wind  speeds. 
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Two  independent  measurements  of  the  surface  waves  were  used.  An  array  consisting 
of  four  Nkhrome  resistance  wire  wavegauges  designed  and  built  at  the  Scripps 
Institution  of  Oceanography  (SIO)  was  deployed  on  a  boom  8  m  above  the  mean  sea 
level  and  17  meters  from  the  hull  The  wires  were  suspended  from  the  ends  of  a  cross 
cuch  that  the  wires  are  0.5  m  from  the  center.  The  ground  wire  was  suspended  from 
the  center  of  the  crow.  A  spreader  was  placed  at  the  lower  end  of  the  wires  to  keep 
them  apart  The  wave  slope  and  directional  wave  spectra  calculations  presented  in 
this  study  use  this  data.  Figure  3.10  shows  a  photograph  of  die  4-wire  wave  gauge 
array  mounted  an  tiie  boa  taken  from  the  deck  of  the  electronics  laboratory  of  Flip. 
Figure  3.11  shows  the  spreader  as  it  is  being  lowered  prior  to  deploying  the  wave 
gauge  army  at  tiie  end  of  the  boom. 

A  capacitance-type  Tantalum  wire  wavegauge  built  at  MIT  based  on  a  design 
developed  at  the  Applied  Physics  Laboratory,  Johns  Hopkins  University  was  also 
deployed  along  tiie  boom  15  meters  from  tiie  huIlofF/p.  It  was  calibrated  in  situ  four 
times  during  tiie  12-day  observation  period.  The  field  calibration  of  tiie  capacitance 
wavegauge  was  accomplished  by  raising  or  lowering  tiie  wire  at  0.5  meter  intervals, 
sampling  the  voltage  for  50  seconds,  and  computing  the  mean  voltage  at  each 
position.  The  resulting  elevation  vs.  voltage  calibration  curve  was  computed  using  a 
linear  least  squares  fit  of  tiie  sample  points  (Figure  3.12). 

An  IMET  Positional  Sensor  Package  located  in  the  working  laboratory  recorded  the 
apparent  acceleration  due  to  the  motion  of  Flip.  It  was  developed  by  the  Upper  Ocean 
Physics  Group  at  SIO  and  it  measures  acceleration  and  tilt  components  in  the 


Figure  3.10.  Photograph  of  the  4-wire  resistance  wave  gauge  array  and  the  capacitance 
wire  wave  gauge  (heavy  black  cable)  deployed  at  the  end  of  the  boom.  The  heavy 
white  cable  near  the  center  of  the  photograpl  the  ground  wire. 
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Figure  3.11.  Photograph  of  die  bottom  spreader  bar  being  lowered  to  position  prior  to 
full  deployment  of  the  wire  wave  gauge  array.  The  bar  is  made  of  stainless  steel  and 
each  arm  is  0.5  meters  long. 
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Figure  3.12.  Plot  showing  the  values  for  the  four  different  field  calibrations  of  the 
capacitance  wire  wave  gauge  The  dashed  line  is  a  linear  least  square  fit  to  the  sample 
points.  The  equation  for  the  best  linear  least  squares  fit  is 

a,,  (m)  *  0.541vc  (Volts)  +  0.800 

where  a,,  is  the  calibration  position  and  vc  is  the  mean  voltage  reading  at  that  position. 


85 


Cktf ittfS 

horizontal  and  vertical  (x,y,z)  directions.  The  signal  from  die  wavegauges,  die 
accelerometers  and  Flip's  compass  were  low-passed  at  4  Hz  and  sampled  at  8  Hz. 


3.5  Wave  spectra. 


The  power  spectra  of  die  time  series  of  die  sea  surface  elevation  a(t)  obtained  from  the 
capacitance  wire  wave  gauge  were  used  to  examine  die  evolution  of  the  wave  field 
during  die  experiment  If  we  assume  that  the  spatial  and  temporal  probability 
distribution  function  of  die  wave  height  is  stationary  over  a  finite  time  period  T,  then 
die  finite  Fourier  transform  pair  (Bendat  &  Piersol,  1986,  p.  130) 


A(f)=  faWe^dt 
Jo 


i 

a(t)*- 1  A(f)e’j2Kftdf 
*Jo 


(3.3) 


can  be  used  to  compute  die  autocorrelation  (or  power)  spectrum  of  die  wave  field 


*(f,T)-iA(f,T)A*(f,T)  (3.4) 


where  oo  -  2xf-  2x/T. 

We  approximated  die  continuous  spectrum  with  the  equivalent  discrete  power 
spectrum  computed  from  the  discretely  sampled  data.  This  was  accomplished  by 
subdividing  the  time  series  into  finite  overlapping  records  of  length  N  with  the 
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overlaps  N/2  records  long.  The  N-point  records  were  then  windowed  using  an  N- 
potnr  Hanning  taper 


w,«0.5ll-cos(2jo7N)]  (3.5) 

to  reduce  siddobe  leakage.  Windowing  the  data  introduces  a  reduction  in  die  total 
energy  in  the  time  series  which  may  be  corrected  by  multiplying  die  Fourier 
transform  with  die  square  root  of  the  variance  of  die  window.  For  a  Hanning 
window,  this  scale  factor  is^§  (Bendat  k  Piersol,  p.  396).  The  discrete  Fourier 

transform  of  die  windowed  time  series  is  therefore  given  by 


N- 1 

A(4)=At^|^^iiwiejM"/w  0.0 

«-0 

where  J = V— 1,  At  is  die  sample  time  step,  and  n  *  -  y y .  The  discrete  Fourier 

transform  was  computed  using  the  Danielson-Lanczos  Fast  Fourier  transform  (FFT) 
algorithm  described  in  Press  et  al  (1986,  p.  390ff).  The  wave  height  power  spectrum 
$(f)  computed  from  A  is  then  given  by 

jA(OA*«.)  (37) 

where  die  asterisk  (*)  denotes  complex  conjugation. 

In  this  experiment,  we  used  a  value  of  N  =  1024.  For  a  sampling  interval  At  -  0.125  s, 
this  gives  a  time  window  T  =  NAt  =  128  s.  To  reduce  die  variance,  die  power  spectrum 
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in  Equation  3.7  was  computed  from  an  avenge  of  K  »  59  power  spectra  of  the  wave 
gauge  time  series.  The  1024-point  data  segments  were  overlapped  512  points  and  the 
resulting  power  spectrum  represented  approximately  an  hour  of  die  data.  The  one¬ 
sided  power  spectrum  which  will  be  used  throughout  this  work  is 

»-Mh  (3-8) 

Figure  3.13  shows  the  characteristics  of  a  typical  wave  height  power  spectrum 
observed  during  the  experiment  The  plot  was  obtained  by  computing  and  averaging 
1024-point  discrete  Fourier  transforms  on  a  1-hour  capacitance  wave  gauge  record. 
Three  major  peaks  are  observed.  Given  the  location  and  fetch  lengths  at  the  site,  a 
significant  peak  due  to  swell  is  present  along  with  the  wind  sea  peak.  The  peak  below 
0.05  Hz  is  due  to  the  27-sec  heave  and  48-second  pitch-roll  resonance  of  Flip  which  are 
unresolved.  The  presence  of  these  peaks  is  consistent  with  previous  studies  of  Flip's 
response  to  waves  (Rudnick,  1967).  While  the  energy  spectrum  of  Flip’s  motion  due  to 
waves  has  components  of  higher  frequency,  their  effect  on  die  measured  wave  energy 
greater  than  0.05  Hz  is  small 


3.6  Wave  slope. 

An  estimate  of  the  root  mean  square  (RMS)  wave  slope  s  was  used  to  parameterize 
the  steepness  of  die  wave  field.  Figure  3.14  shows  a  plan  view  schematic  diagram  of 
die  location  and  orientation  of  die  4-wire  resistance  wave  gauge  array.  We  defined  the 
East  direction  to  be  die  true  x-axis  and  the  North  direction  to  be  die  true  y-axis.  A  local 
coordinate  system  based  on  the  geometry  of  Flip  was  also  defined  such  that  die 
direction  towards  the  keel  was  the  Flip  x-axis  xp.  The  orientation  of  the  Flip  y-axis  yp 
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Figure  3.13.  Plots  showing  the  characteristics  of  a  typical  wave  height  power 
spectrum.  The  data  were  sampled  at  8  Hz.  The  spectra  were  computed  from  a  one 
hour  wave  gauge  record  with  Ujq  *  8  m/s.  a)  Power  spectrum  computed  from  1024- 
point  FFTs  used  in  throughout  this  work,  b)  Power  spectrum  over  the  same  period 
computed  from  2560-point  FFTs  to  resolve  the  lower  frequency  spectral  peaks. 
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follows  from  the  right-hand  rule.  A  third  coordinate  system  based  on  the  wire  wave 
gauge  array  was  also  defined  such  that  the  array  x-axis  xw  bisects  the  angle  formed  by 
the  arms  from  which  wires  1  and  4  were  suspended  This  axis  is  also  perpendicular  to 
the  boom  and  the  array  y-axis  ywiscobnear  with  the  boom. 


The  orientation  of  the  boom  and  of  the  wave  gauge  array  are  given  in  Figure  3.14.  The 
angle  between  the  keel  of  Ffijp  and  the  magnetic  North  6p  was  measured  by  a  compass 
on  Flip.  The  signal  was  sampled  at  8  Hz  together  with  the  wire  wavegauge  signal. 
Prom  Figure  3.14,  the  angle  between  the  true  coordinate  system  and  the  array 
coordinate  system  Is 


e-eF-fnsP+sa0) 

e=eF-209°  (3.9) 

where  0  is  defined  to  be  positive  in  the  clockwise  direction.  The  instantaneous  slope  of 
the  sea  surface  s(t)  *=  (s^)  can  be  estimated  from  the  simultaneous  sea  surface 
elevation  a^t)  measurements  of  the  four  wire  wave  gauges,  i=l  A3,4.  In  the  local 
coordinate  system  (xw/yw),  the  slope  is  computed  from  the  mean  of  the  finite 
difference  of  the  wave  elevation  between  two  wires 


(t)  .[«,(«)-m»1+M>-m«>1 

2**  (3.10) 
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where  Ax  ■  Ay  -  0.7071  m.  The  wave  slope  in  the  true  coordinate  system  is  therefore 


•a«(0-iI«(Ocosfl(t)-sJy(t)stafl(t)  ^ 

*a,(t)  *  •2y(0«lnfl(t)+*Jl(t)co*6(t) 

The  one  hour  variance  of  the  wave  slope  components,  o2( s^)  and  o^s^y),  where 


were  computed  from  Equation  3.11.  The  RMS  wave  slope  s  is  therefore 


(3.12) 


3.7  Directional  wave  height  spectra. 

3.7.1  Theory. 

We  used  the  Maximum  Likelihood  Method  in  computing  directional  spectra  using  the 
4-wire  wave  gauge  array.  The  method  was  first  used  by  Capon  (1969)  in  computing 
die  wavenumber  spectra  of  seismic  waves.  Oakley  &  Lozow  (1977)  showed  that  in  a 
nearly  noise  free  environment,  a  sparse  wave  gauge  array  of  at  least  three  wires  can 
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theoretically  resolve  the  direction  of  waves  whose  wavelengths  are  several  orders  of 
magnitude  larger  than  the  size  of  the  array.  However,  the  presence  of  noise  in  real 
arrays  sets  a  finite  limit  to  the  resolving  power  of  fire  array.  We  present  a  derivation 
of  fire  Maximum  likelihood  technique  for  obtaining  five  directional  wave  height 
spectrum  tyfojb)  from  an  array  of  N  wire  sensors  based  on  an  approach  described  by 
Jeffreys,  Wareham,  Ramsden  &  Platts  (1981 ).  This  section  is  essentially  a  review  of 
their  paper. 


Let  A|(<d)  be  the  Fourier  transform  of  fire  wave  elevation  time  series  for  wave  gauge  i 
at  position 


3*r,(cos8„sind))  (3.13) 

where  rt  and  6j  are  the  distance  and  angle  from  the  origin  to  the  sensor  i,  i*l,...,N.  The 
NxN  cross  spectral  density  matrix  C  for  sensors  i  and  j  is  then  given  by 


Aj(q>)  Aj(q>) 
|Ai(ca)|  |Aj(oo)| 


(3.14) 


The  overbar  in  file  above  equation  indicates  that  a  sufficient  number  of  Fourier 
transforms  should  be  averaged  to  ensure  that  the  cross  spectral  density  estimate 
Cjj(co)  is  statistically  reliable.  The  asterisk  *  indicates  complex  conjugation.  It  is 

common  to  normalize  file  Fourier  transforms  by  their  moduli  to  avoid  errors  due  to 
differences  in  wave  gauge  calibration. 
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If  the  directional  wave  height  spectrum  <b((o,6)  is  expressed  as  the  sum  of  M  plane 
waves  of  power  propagating  along  directions  6m,  then  the  true 

cross  power  spectral  density  function  is  then 


u 

C(o)  -  y>x(tD,e,)<IKa),e,)jt*T(<a,e,).  (3.15) 

*(a>.e«)«e'*-r'  i=l,...,N.  (3.16) 

The  ith  component  of  x  is  the  complex  phase  lag  between  the  ith  wave  gauge  and  the 
origin  for  a  wave  of  frequency  oo  propagating  along  die  direction  0m  while  k  is  the 
wavenumber  vector  whose  dements  are 

«-^(cos0m,sin6m).  (3.17) 

We  implicitly  assumed  the  linear  deep  water  dispersion  relation  co2=gk. 

Oven  the  matrix  C(co)  estimated  using  Equation  3.14,  we  can  invert  Equation  3.15  to 
obtain  an  estimate  of  die  power  spectrum  9(ca,6d)  of  the  wave  propagating  along 
direction  6d 

4K<».8d)«B*JC(a>)Bd 

M 

4<o.emf  *((o,eJ 

mrnl 

M 

»|b*J  *(a>.ed)|2d<<o.0d)+  xfoe^Vco.eJ. 

n>M 
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To  evaluate  $(a>,ed),  we  need  to  determine  the  vector  Ej  such  that  the  contribution 
from  ^((0,6^  is  undistorted  while  the  contributions  from  die  other  directional 
components  ate  minimized.  We  minimize  $(00^4)  with  respect  to  Ed  subject  to 

|E-Jx(m,ed)|-1. 

This  leads  to  an  estimate  of  die  power  in  die  wave  frequency  ©  propagating  along  0d 
(Capon,  I960) 


4>(©.0d)*  [K*T(©,fid)Cr1(©)*(©,0d)]"  .  (318) 

The  directional  spectrum  was  scaled  such  that  die  sum  of  the  wave  energy  at 
frequency  oo  for  all  the  wave  component  directions  is  equal  to  die  spectral  energy  level 
of  die  capacitance  wire  wave  gauge 


M 

(3I9) 

nisi 

Problems  regarding  the  numerical  stability  of  the  C  matrix  can  arise  when  its  inverse 
is  computed  in  Equation  3.18.  If  there  are  fewer  wave  frequencies  than  wave  gauged 
then  C  will  be  overdetermined  and  consequently  singular.  It  can  be  made  positive 
definite  and  its  stability  enhanced  by  introducing  a  small  positive  perturbation  R  « 1 
such  that 
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C*«(l-R)C+ja  (3.20) 

where  I  is  an  NxN  identity  matrix.  This  method  is  equivalent  to  introducing  a  small 
uncorrelated  noise  whose  power  is  R  times  the  wave  energy  into  the  wave  gauges. 

3.7J2  Implementation. 

We  now  outline  the  basic  steps  in  implementing  a  practical  algorithm  for  computing 
the  directional  spectrum  using  die  maximum  likelihood  method  front  the  time  series 
of  die  surface  wave  elevation  measurements  aj(t)  at  wave  gauge  i  for  the  4-wire 
resistance  wave  gauge  array  shown  in  Figure  3.15. 

1.  Compute  the  complex  phase  lags  jq  (Equation  3.16)  for  each  wave  gauge  using  the 
array  coordinate  system  shown  in  Figure  3.15. 

2.  Subdivide  the  time  series  data  aj(t)  into  sample  windows  of  length  In  this 

experiment,  we  used  Nw  - 1024  with  a  512-point  overlap  between  sample 
windows.  Since  die  data  was  sampled  at  8  Hz,  each  window  is  128  seconds  long. 

* 

3.  Window  die  data  to  reduce  sidelobe  leakage.  A  Hanning  window  (Equation  3.5) 
was  used. 

4.  Compute  the  discrete  Fourier  transform  A(fn)  of  the  windowed  time  series  using 
the  Fast  Fourier  Transform  method. 
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5.  Compute  the  4x4  cross  spectral  density  matrix  C(a>)  (Equation  3.14).  Each  C^(cd) 
estimate  was  based  on  32  minutes  of  data  or  29  spectral  averages. 


6.  Perturb  die  C(a>)  matrix  by  introducing  a  small  noise  R  (Equation  3.20).  R  «  10*6 
was  used. 


7.  Compute  die  inverse  of  C(a>). 

8.  Subdivide  die  0  space  into  M  directions.  Compute  die  Maximum  Likelihood 
estimate  of  die  wave  spectrum  4(00,64)  in  direction  64  (Equation  3.18)  few  all  die  M 
directions.  We  used  M  -  72  in  this  experiment  resulting  in  a  A0  =  5°. 

9.  Compute  the  power  spectral  density  of  the  capacitance  wire  wave  gauge  <b(  tn)^  ■ 

10.  Rescale  the  4  computed  in  step  8  so  that  the  sum  of  all  wave  directions  at  a  given 

frequency  oo  is  equal  to  4(co)  (Equation  3.19). 

11.  Convert  die  directional  spectrum  from  the  wave  gauge  array  coordinate  system  to 

die  trut  coordinate  system.  From  Figure  3.15,  the  azimuth  6  of  a  given  direction  is 

0-0,-119*+^  (3.21) 

where  84  is  the  angle  of  that  direction  from  the  array  x-axis  xw/  and  6p  is  the  azimuth 

of  Flip's  keel. 
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4  OBSERVATIONS 


4,1  CSmmmI  coadWoM 

The  observation  period  of  this  experiment  covered  approximately  twelve  days 
between  September  26, 1991  -  October  8, 1991  Cuban  day  269  -  281).  We  began 
recording  the  acoustic  measurements  while  the  winds  blowmg  from  the  south  were 
dying  down.  After  cme  and  a  half  days  of  relative  calm,  the  wind  blew  initially  from 
die  north  (fetch  length:  600  km)  with  die  wind  apeed  increasing  over  a  day.  It  stayed 
in  a  quasi-steady  condition  over  a  period  of  six  days.  Although  die  variability  in  the 
wind  direction  was  minimal  during  die  quasi-steady  period,  die  hourly-averaged 
wind  speed  U  fluctuated  between  6  -12  meters/second.  The  wind  then  changed 
direction  and  blew  from  die  south  (fetch  length:  unlimited)  and  then  died  down  over 
the  course  of  a  day  and  a  half.  The  sudden  change  in  wind  direction  was  accompanied 
by  a  rapid  decline  and  increase  in  the  wind  speed  over  a  period  of  six  hours. 


Wind  speed  U10  and  friction  velocity  u*. 


Although  our  wind  speed  measurements  were  made  at  28  meters  elevation  (U^,  we 
reduced  this  value  to  the  wind  speed  at  10  metes  consistent  with  meteorological 
convention.  It  is  common  to  assume  a  logarithmic  wind  speed  profile  near  die  sea 
surface  of  the  form  (Large  &  Pond,  1981) 

<4i> 
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where  U,  is  the  wind  speed  at  elevation  z,  u.  is  the  friction  velocity,  x  »  0.40  is  the  von 
Karman  constant,  zq  is  the  sea  surface  roughness  length  scale  and  £  is  a  term  that 
accounts  for  the  effect  of  atmospheric  stability  on  die  wind  speed  profile.  Since  u»  was 
not  directly  measured  in  die  experiment,  its  value  was  inferred  from  U^- 

To  evaluate  H  using  the  bulk  aerodynamic  method  (Large  k  Pond,  1981), 
measurements  of  die  sea  surface  temperature  and  the  humidity  are  required. 
Unfortunately,  the  instruments  for  measuring  both  parameters  failed  to  operate 
during  the  experiment  so  we  were  unable  to  evaluate  this  term.  Large  fc  Pond  (1981) 
showed  that  assuming  neutral  stability  (E  *  0)  introduces  an  uncertainty  of  less  than 
4%  in  die  value  of  u*  and  Ur  We  neglected  the  effect  of  utmospherk  stability  cm  the 
wind  speed  profile  and  used  the  equation 

Uxa— ln_l.  (4.2) 

K  Zq 

Maat,  Kraan  &  Oost  (1991)  suggested  that  zq  is  a  function  of  the  wave  age  (Cp/u.) 
where  Cp  is  die  phase  speed  of  the  peak  frequency  of  the  wave  spectrum  computed 
from  the  wave  gauge  record.  Using  an  empirical  fit  to  field  data  obtained  during  die 
1986  Humidity  Exchange  Over  the  Sea  Main  Experiment  [HEXMAXJ,  they  proposed 
that 


*0  =  0-8^ 
g 


(4.3) 


From  die  two  equations  (Equations  4.2  and  4.3),  we  can  determine  the  value  of  the 
two  unknowns  u»  and  zq.  We  computed  u*  and  zq  by  assuming  a  value  for  u*  in 
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270  272  274  276  278  280 

Julian  day 


Figure  4.1.  Time  series  of  Ujo  u»  and  the  wind  direction  during  the  NOBS 
experiment  Note  that  die  wind  rapidly  changes  direction  on  Julian  day  278  after 
blowing  from  the  North  for  approximately  six  days. 
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Equation  43  and  then  cubatituting  the  computed  value  of  zq  into  Equation  4.2  to  get  a 
new  u*  estimate.  The  process  was  repeated  until  the  difference  between  the  estimated 
and  computed  values  for  u»  differed  by  less  than  0.001  m/s.  The  value  of  U10  was 
then  computed  from  Equation  4.2.  Figure  41  shows  time  series  of  hour-averaged  U^y 
u»  and  die  wind  direction  (hiring  the  observation  period. 


As  we  noted  in  Section  3.5,  die  wave  gauge  record  was  contaminated  by  low 
frequency  energy  due  to  the  motion  of  the  platform.  Figure  3.12  shows  that  a 
significant  proportion  of  this  energy  can  be  eliminated  by  specifying  a  high  pass 
frequency  cutoff  at  0.05  Hz,  and  that  the  contribution  erf  the  platform  motion  to  the 
wave  gauge  record  above  0.05  Hz  is  negligible  compared  to  die  wave  energy.  Daily 
plots  of  die  one  hour  wave  height  spectrum  for  the  time  series  beginning  at  0000  UT 
throughout  the  observation  period  are  shown  in  Figure  4.2. 

The  time  series  erf  the  integral  characteristics  of  die  wave  spectrum  are  shown  in 
Hgure  4.3.  The  root-mean  square  wave  amplitude  a  and  die  mean  frequency  f  of  the 
wave  spectrum  were  defined  as 
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Figure  4.3.  a)  Root  mean  square  amplitude  a  (m)  and  significant  wave  height  H,  (m) 
computed  from  the  full  spectrum,  b)  Mean  frequency  computed  from  die  full 
spectrum  f  (Hz),  c)  Up,  (m)  and  H,  (m)  computed  with  swell  filtered  out.  d)  Mean 
frequency  with  swell  filtered  out  f*  (Hz). 


105 


(ui)  sH  (Ui)  «h 


log  $ 


Figure  4.4.  Time  series  of  log®  (m2/Hz)  at  selected  frequencies:  0.25  Hz  -  O, 0.5  Hz  - 
•,  0.75  Hz  -  a,  1  Hz  -  ■. 
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Figure  4.5.  lime  series  of  wave  slope  s  (Equation  3.12)  computed  from  the  resistance 
wire  wave  gauge  array  data. 
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In  Figures  4.3a  and  4.3b,  the  lower  limit  of  integration  in  Equations  4.4  and  45  was  set 
at  fo-0.05  Hz.  The  values  of  a  and  F  in  both  plots  include  spectral  energy  from  both 
wind  waves  and  swell.  In  Figures  43c  and  4.3d,  the  lower  limit  of  integration  used 
was  die  frequency  between  die  swell  and  the  wind  sea  peak  in  which  0(f)  was  a 
minimum.  This  was  chosen  so  that  the  energy  of  die  swell  does  not  affect  die  value  of 
Op  and  Tw.  The  plot  shows  that  die  wind  wave  spectrum  Fw  ranged  from  0.15  to 
0.4  Hz  during  NOBS. 

Figures  4.3a  and  43c  also  shows  the  wave  height  data  in  terms  of  the  significant  wave 
height  H„  defined  as  the  mean  of  the  highest  one-third  of  waves.  For  a  Rayleigh- 
distributed  wave  (Kinsman,  1984  p-  390), 

Hf  =  Z83a.  (46) 

Figure  44  shows  time  series  of  the  spectral  level  of  <X>  at  0.25, 0.5, 0.75  and  1.0  Hz.  The 
plots  in  Figure  44  are  a  measure  of  the  variability  of  the  wave  energy  above  the  wind 
peak. 

The  method  for  computing  the  root  mean  square  wave  slope  s  using  the  4-wire  wave 
gauge  array  was  described  in  Section  3.6.  The  time  series  of  $  is  shown  in  Figure  4.5. 


44  Ambient  noise 

Figure  4.6  presents  the  time  series  of  the  hour-averaged  noise  spectrum  level  N  at 
selected  frequencies  [see  Figure  3.11  for  plots  of  typical  N(f)].  We  chose  to  display  N 
at  4.3, 8.0,  and  14.0  kHz  since  these  are  the  frequencies  WOTANs  commonly  used  (see 
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Figure  4.6.  Time  series  of  the  ambient  noise  spectrum  level  N  (dB  re  1  uPa2/l  Hz)  at 
selected  WOT AN  frequencies.  Note  the  similarities  between  these  time  series  and  the 
U]q  time  series.  In  particular,  note  the  rapid  drop  and  increase  in  N  on  Julian  day  278. 
a)  4.3  kHz.  b)  8.0  kHz.  c)  14.0  kHz.  d)  wide  band  [2.5-14  kHz  or  X]. 
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Table  1.2).  The  total  integrated  acoustic  energy  between  2.5-14  kHz  denoted  by£is 
also  presented.  The  gaps  in  the  acoustic  time  series  are  due  to  records  eliminated 
because  some  of  the  devices  used  in  deploying  instruments  unrelated  to  our 
experiment  generated  electrical  noise  that  contaminated  the  full  acoustic  spectrum. 
Figures  4.1  and  4.5  show  that  the  fluctuation  in  the  wind  speed  Ujo  and  u.  is  clearly 
reflected  in  N.  The  high  correlation  between  Ujq  and  N  is  consistent  with  the 
investigations  we  reviewed  in  Section  1.1.  In  addition,  good  correlations  between  N 
and  properties  of  the  wave  height  spectrum  O,  particularly  the  high  frequency 

wave  spectrum  components  in  Figure  4.4  and  wave  slope  s  in  Figure  4.5,  and  U10  and 
N  is  also  apparent  in  Figures  4.3  -  4.6.  This  will  be  clearer  in  Section  4.6  when  we 
examine  the  empirical  relationships  between  N  and  environmental  parameters. 

4.5  Directional  wave  height  spectra. 

The  method  for  computing  the  directional  wave  height  spectra  was  outlined  in 
Section  3.7.  The  maximum  horizontal  dimension  of  the  array  is  1  meter.  The  smallest 
wavelength  the  array  can  resolve  is  of  this  order.  Consequently,  the  highest  wave 
frequency  that  can  be  resolved  is  approximately  1  Hz. 

Figure  4.1  shows  that  although  the  wind  speed  has  fluctuated  between  5  to  11 
meters /second,  the  wind  direction  has  been  quasi-steady  (from  approximately  30°) 
for  6  days  prior  to  Julian  day  278.  The  wind  then  began  to  change  direction  and  die 
down  with  U10  decreasing  to  4  meters /second  while  turning.  After  several  hours,  the 
wind  had  completely  turned  and  was  blowing  from  the  South  (180°). 

The  evolution  of  the  wave  field  during  this  event  is  shown  by  a  series  of  contour  plots 
of  directional  wave  spectra  <b(W)  in  Figure  4.7.  The  plots  have  a  directional  resolution 
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of  5*  and  a  frequency  resolution  of 0.002  Hz.  They  show  the  directional  distribution  of 
the  wave  energy  within  the  frequency  range  0.25  to  1  Hz.  Since  the  wave  energy 
spectrum  4>(f)  decays  as  a  function  of  frequency,  only  the  wave  energy  at  die  lower 
frequency  limit  are  seen  in  these  plots.  In  addition,  since  the  directional  wave  spectra 
plots  show  die  direction  die  waves  are  moving  to  while  the  anemometer  gives  the 
direction  the  wind  is  blowing  from,  we  added  180°  to  die  wind  direction  to  facilitate 
comparison.  The  figure  lists  the  start  time,  die  wind  speed  and  the  direction  the  wind 
is  Mowing  towards  for  each  plot 

The  wind  direction  6,^  and  wind  speed  are  relatively  steady  in  die  first  six  plots. 
The  wind  begins  to  turn  and  decay  in  Figure  4.7g  with  Ujq  reaching  a  minimum  in 
Figure  4.7L  During  this  period,  die  mean  wave  direction  in  the  energy  containing 
portion  of  the  spectrum  and  die  wind  direction  are  oriented  in  different  directions.  In 
Figure  4.7n,  we  see  some  wave  energy  at  approximately  350°  and  0.45  Hz.  This 
matches  die  wind  direction  for  that  plot  which  is  at  354°.  In  the  next  plot,  the  energy 
at  this  direction  has  increased  and  begun  to  migrate  towards  the  lower  frequencies.  In 
die  remaining  plots,  die  wave  energy  centered  around  the  original  direction  (200°)  has 
decayed  considerably  while  the  wave  energy  parallel  to  die  wind  has  already  grown 
significandy. 

While  these  plots  are  informative,  the  information  content  is  biased  towards  the  lower 
frequencies  because  die  energy  level  there  is  higher  than  at  the  higher  frequencies.  It 
is  possible  to  reduce  this  bias  and  observe  die  high  frequency  directional  behavior  of 
the  wave  field  by  eliminating  the  influence  of  the  spectral  energy  level  on  the  plot. 

The  same  plots  in  Figure  4.7  were  plotted  in  Figure  4.8.  However,  instead  of  equating 
the  total  wave  energy  at  each  frequency  with  the  spectral  energy  level  of  the 
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Figure  4.7.  Directional  wave  spectra  from  0.25  to  1  Hz.  The  sequence  of  plots 
demonstrate  die  evolution  of  die  wavefield  in  response  to  turning  winds  on  Julian 
day  278.  Hie  x-axis  for  die  plots  is  frequency  in  Hertz  and  die  y-axis  is  wave  direction 


0^.  Dashed  lines  indicate  the  direction  the  wind  is  blowing  towards  (6, 


Start  Time  (UT)  U10(m/s)  0*^ 


8.0  211° 


8.1  211° 


0511  T  7.7  T  211° 


&2  212° 


°D>i6SiU9m 


Start  Time  (UT) 


U10(m/s) 

4.6 


115 


Figure  4.8.  The  directional  wave  spectra  in  Figure  4.7  normalized  by  the  spectrum 
level  at  the  given  frequency  to  highlight  die  directional  distribution  of  the  wave 
energy  in  response  to  turning  winds  on  Julian  day  278.  The  x-axis  for  the  plots  is 
frequency  in  Hertz  and  die  y-axis  is  die  wave  direction  64.  Dashed  lines  indicate  the 
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capacitance  wire  wave  gauge  ( Step  10  in  Section  3.7.2  and  Equation  3.19),  we  instead 
chose  to  normalize  <D(f,8)  such  that 

M 

^«Kf.eJ-i.  (4.7) 

m-1 

Figure  4.8  shows  contour  plots  of  how  die  wave  energy  is  distributed  as  a  function  of 
frequency  subject  to  Equation  4.7.  The  response  of  die  high  frequency  wave 
components  to  die  change  in  wind  direction  is  initiated  in  Figure  4.8i-l  when  the  wind 
speed  decreases  to  approximately  4.0  meters/second.  Prior  to  this,  the  wind  has  been 
blowing  along  220°  for  several  days  and  die  wave  direction  is  correspondingly 
centered  around  that  direction  as  well.  In  Figure  4.8m,  we  observe  a  region  of  high 
energy  centered  at  350°  for  f  >  0.65  Hz.  Note  that  this  is  not  observed  in  Figure  4.7m 
because  of  the  spectral  energy  distribution.  Tie  new  wave  direction  matches  die 
increasing  wind  speed  at  that  direction  (354°).  The  plot  also  suggests  that  die  high 
frequency  wave  components  respond  quickly  to  changes  in  wind  speed  and  wind 
direction.  As  die  wind  speed  increases  and  time  progresses,  this  high  energy  region 
grows  and  migrates  towards  the  lower  frequencies  until  the  wave  direction  is  aligned 
with  the  wind  direction  across  all  frequencies  (Figure  4.8r-t). 

4.6  Discussion 

4.6.1  Ambient  noise,  UJ0  and  «•. 


Figure  4.9  shows  plots  of  log  Ujo  versus  die  spectrum  level  at  the  different  ambient 
noise  frequencies  N(f)  given  in  Figure  4.6  showing  a  linear  relationship  between  the 
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Tabic  4.1.  Some  coefficients  for  the  WOT  AN  expression  and  the  equivalent  power 
laws  computed  from  NOBS  data  (ci  Table  li  and  Figure  4.9). 


logU,0  ■  mN+  n  -4-»mUjo 
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2.7 

two  parameters.  We  also  plotted  the  linear  regression  equations  from  the  different 
WOT  AN  measurements  given  in  Table  1.2  and  from  our  data.  Regression  lines  in  the 
form  of  the  WOT  AN  equation  (Equation  1.3)  were  computed  from  the  data  using  an 
orthogonal  least  squares  fit  (Casella,  1990,  p.  584)  since  both  log  U10  and  N(f)  are 
random  variables.  A  summary  of  tite  linear  regression  coefficients  and  the  power  law 
coefficients  between  the  two  parameters  is  given  in  Table  4.1.  The  equation  of  the 
regression  lines  for  our  drta  are  comparable  to  those  found  by  previous  investigators. 
Our  results  suggest  an  ambient  sound  pressure  level  dependence 

p2-Ufo4_a-°.  (48) 

Figure  4.10  shows  four  plots  of  N  versus  log  u»  sampled  and  computed  from  the  wind 
speed  measurements  in  this  experiment.  The  data  is  well  correlated.  It  suggests  that  a 
single  linear  fit  to  the  data  is  adequate  and  that  the  power  law  expression 
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and  N  and  the  equivalent  power 
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(4.9) 


best  describes  the  data  (see  Table  42).  Our  results  contradict  Kerman  (1984)  who 
argued  based  on  empirical  data  that  p2  ~  u.*-5  when  u*  >  u*c  *  0.23  m/s,  and  p2  -  u»3 
when  u»  <  u*c  However,  the  scatter  in  Kerman's  (1984)  data  at  u.  <  u.c  did  not  clearly 
show  that  there  are  indeed  two  subranges. 

Oven  that  u»  was  computed  primarily  from  Uz  in  both  Kerman  (1984)  and  in  this 
work  (Equations  42  and  43)  and  is  only  weakly  dependent  on  Cp,  the  relationship 
between  u*  and  N  is  essentially  similar  to  that  between  N  and  U10. 

In  Section  41,  we  described  an  implicit  method  using  Maat  et  al's  (1991)  equations  for 
the  relationship  between  the  value  of  u*,  Uz  and  the  wave  age.  It  is  desirable  to  obtain 
an  explicit  and  dimensionally  consistent  expression  for  u*  as  a  function  of  Ujq  that  is 
equivalent  to  their  method.  This  can  easily  be  done  from  our  results. 


127 


Figure  4.11a  shows  a  log-log  plot  of  die  drag  coefficient  C0  « (u./  VI0)2  versus  the 
inverse  wave  age  (U10/Cp)  showing  an  increase  of  the  sea  surface  roughness  with 
increasing  wind.  In  this  figure  we  see  that  die  relationship  between  die  two  variables 
can  be  approximated  by  a  power  law 


U 

By  rearranging  die  variables  in  Equation  4.10,  we  can  obtain  an  expression  relating 
U](y  u»  and  Cp.  Figure  4.11b  plots  die  two  versions  of  die  inverse  wave  age  expression, 
(Ujo/Cp)  mid  (u*/ Cp),  against  each  other.  The  regression  line  on  die  plot  translates  to  a 
power  law  expression 


\0j6 

j  ' 


(4.10) 


^=3.85x10-*^^  , 
u.  *3.85x10-2U^. 


(4.11) 


4.6.2  Ambient  noise  and  wave  parameters. 

The  previous  field  experiments  reviewed  in  Section  1.2  [Penhallow  k  Dietz  (1964), 
Perrone  (1969),  Lemon  k  Fanner  (1984)]  examined  die  relationship  between  ambient 
noise  and  die  significant  wave  height  (or  alternatively,  the  root  mean  square 
amplitude  a)  and  found  that  these  two  parameters  poorly  correlated.  The  plots  in 
Figure  4.12  which  compares  our  measurements  of  N  and  a  is  consistent  with  their 
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results.  We  suggested  that  the  presence  of  acoustically  inactive  swell  affects  the 
relationship  between  N  and  s.  While  the  presence  of  swell  significantly  increases  *,  it 
does  not  proportionally  increase  the  steepness  of  fee  waves  nor  does  it  proportionally 


To  reduce  the  influence  of  the  swell  energy  on  the  variance  of  die  wave  field,  we 
filtered  die  swell  out  of  the  wave  height  reconi  and  computed  the  root  mean  square 
amplitude  of  the  data  (see  Section  43).  The  resulting  estimate,  whose  time  series  is 
plotted  in  Figure  43c,  was  plotted  against  N  (Figure  413).  In  this  figure,  we  see  that 
although  die  scatter  in  the  data  is  greater  dun  that  between  Ui0  and  N,  the 
improvement  in  the  correlation  between  a  and  N  is  significant  Table  43  shows  that 
die  power  law  relationship  is  given  by 


p2 -<£*-**  (412) 


We  can  better  understand  die  relationship  between  the  wind  speed  and  wind  waves 
by  examining  the  behavior  of  the  high  frequency  wave  components.  In  Figure  4.4,  we 


Table  U.  Linear  regression  coefficients  relating  and  N,  and  die  equivalent  power 
laws  computed  from  NOBS  data. 


logo*  ■  mN+  n  p2  =  nu£ 
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8.0 

0.89 
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3.3 

2.0 

140 

0.77 

0.0547 
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1.8 
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0.0494 
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showed  the  time  series  of  the  wove  spectrum  4(f)  far  f  -  0.25, 0.5. 0.75,  and  1 .0  Hz. 
Wenoted  that  the  4(f)  time  series  doaely  match  the  shape  of  both  the  U10  and  N  time 
series.  The  behavior  of  the  4(f)  resembles  Uiq  and  N  better  than  die  s*  time  series. 
Since  the  wind  wave  mean  frequency  J  is  between  0.15  and  0.4  Hz  (Figure  42d),  then 
frequencies  greater  than  0.4  Hz  are  within  the  wind  wave  frequency  range. 

In  the  last  section,  we  observed  from  die  directional  wave  spectrum  plots  (Figure  4.7 
and  4.8)  that  die  higher  frequency  wind  wave  components  (f  >  OS  Hz)  respond 
quickly  to  changes  in  wind  speed  and  direction.  Figure  414  shows  plots  oflog  4  at 
0  J,  0.75,  and  1.0  Hz  versus  N.  They  show  a  well  correlated  relationship  between  these 
frequencies  and  N  except  the  lower  4(0.5  Hz)  values  (low  wind  speed, 

4(0 J5  Hz)  <  >1.0)  where  die  scatter  and  nonlinearity  in  the  data  is  more  considerable. 
However,  die  scatter  in  rest  of  die  4(0.5  Hz)  data  is  considerably  smaller. 

Our  results  show  that  a  and  N  are  poorly  correlated  partly  because  of  the  presence  of 
swell  The  characteristics  of  swell  are  not  directly  determined  by  the  local  wind-wave 
relationship  and  its  influence  on  the  local  noise  generation  mechanisms  is  weak.  By 
neglecting  swell  and  concentrating  on  the  wind  sea  portion  of  the  wave  spectrum,  we 
can  obtain  a  better  correlated  power  law  relationship  between  N  and  a.  This  is  clearly 
shown  by  improved  correlation  between  die  energy  at  die  higher  frequency 
components  of  4  and  N  (Figure  414). 

Experiments  by  Rapp  k  Melville  (1991)  and  Loewen  k  Melville  (1991a)  have 
suggested  that  the  wave  slope  of  a  packet  of  M  waves 
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N(  14kHz)  N(4.3kHz) 


a 


a 


Figure  412.  Plot  of  N  vs.  the  root  mean  square  wave  amplitude  a.  The  scatter  in  the 
data  give  correlation  coefficients  «  0.67  for  all  the  plots. 


Figure  413.  Wot  of  N  vs,  a*.  The  surface  elevation  time  series  was  high  pass  filtered 
to  eliminate  the  contribute  of  swell  energy  to  a.  The  figure  shows  that  compared  to 
Figure  4.12,  we  see  a  significant  improvement  in  file  correlation  between  ambient 
noise  N  and  file  root  mean  square  wind  wave  amplitude  a„  (see  Table  4.3). 
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Figure  4.14.  c)  Plot  of  N  vs.  $(1.0Hz). 
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Chapter  4 


can  be  used  to  measure  the  strength  of  wave  breaking  in  the  laboratory.  Their 


experiments  also  suggest  that  Sp  correlates  with  several  important  dynamical 


properties  of  wave  breaking,  eg.,  loss  erf  excess  momentum  flux,  mixing,  dissipation, 
noise  generation.  Their  results  suggest  that  the  overall  steepness  of  the  wave  field  is 
an  important  indicator  of  the  characteristics  of  the  breaking  wave.  We  can  show  that 
tfie  steepness  of  die  wave  Add  correlates  with  some  dynamic  characteristics  of  wave 
breaking  in  die  ocean  as  well 


The  time  series  of  the  RMS  wave  slope  s  was  given  in  Figure  4.5.  A  comparison  of  this 
plot  with  those  of  the  U10  and  N  time  series  (Figures  4.1  and  4.5)  shows  that  the 
general  characteristics  of  those  plots  can  also  be  observed  here.  In  Figure  4.17  we 
presents  four  plots  of  log  s  versus  N.  The  plots  show  that,  consistent  with  the 
laboratory  results  of  Melville  et  al  (1988)  and  Loewen  k  Melville  (1991a),  the 
steepness  of  die  wave  field  correlates  well  with  the  ambient  noise.  The  characteristics 
of  this  relationship  are  summarized  in  Table  4.4  which  gives  correlation  coefficients 


N 

y,U-*Xy<-y) 
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Tabic  4.4.  Power  law  coefficients  for  flic  ambient  sound  pressure  level  and  wave 
slope  s,  and  die  corresponding  correlation  coefficients  between  N  and  log  s. 
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3.1 

0.85 

greater  than  0.85  between  N  and  log  s  except  for  the  14  kHz  noise.  The  plots  of  wave 
slope  versus  N  give  a  power  law  relationship 

p2_^-14  (4.13) 

We  excluded  the  N(14  kHz)  regression  because  of  the  significantly  poorer  correlation 
coefficient  with  this  ambient  noise  frequency.  As  we  mentioned  in  Section  1.1,  field 
experiments  [Fanner  &  Lemon  (1984),  Vagle  et  al  (1990)]  have  suggested  that  the 
influence  of  sound  absorption  by  bubbles  can  significantly  alter  the  ambient  noise 
levels  above  8  kHz. 

Figure  4.16  shows  four  plots  relating  s  to  u*  and  wave  age  parameters  based  on 
both  Ujo  and  u*.  While  most  of  the  data  points  in  the  figure  are  well  correlated,  the 
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scatter  in  the  data  is  more  significant  than  in  Figure  415.  The  scatter  for  data  points 
for  U10  <  5.5  m/s  (u*  <  0l2  m/s]  appears  to  be  greater  than  for  data  points  above  it 

4.63  Ambient  Moist  and  wind-wave  development. 

In  making  simultaneous  wind,  wave  and  ambient  noise  measurements,  Penhallow  & 
Dietz  (1964)  subdivided  their  data  into  portions  where  die  wind  was  steady,  increasing 
and  decreasing.  Their  results  show  (see  Section  1 .2)  that  die  linear  correlation 
coefficient  between  the  ambient  noise  at  630  Hz  and  the  wind  and  wave  variables  is 
higher  when  die  wind  is  steady  and  degrades  considerably  when  the  wind  is 
(hanging.  They  theorized  that  when  die  winds  are  steady  and  the  waves  are  fully 
developed,  there  is  an  equilibrium  between  die  rate  of  energy  input  by  the  wind  to 
the  waves  and  die  rate  of  energy  dissipated  by  waves  due  to  wave  breaking.  It  can 
then  be  reasoned  that  since  die  wind  and  ambient  noise  relationship  is  well-defined 
and  well-correlated,  then  die  wave  conditions  that  satisfy  this  equilibrium  should  be 
well-correlated  with  ambient  noise  as  well.  Unfortunately,  there  are  no  plots  or 
figures  in  their  paper  with  which  we  can  examine  the  data  nor  was  their  research 
followed  up  in  the  literature.  As  we  mentioned  in  our  review,  die  research  in  wave 
and  ambient  noise  relationship  is  relatively  sparse. 

We  decided  to  further  explore  the  ideas  in  Penhallow  &  Dietz  (1964)  by  subdividing 
die  available  NOBS  data  and  performing  a  similar  comparison.  While  the  boundary 
between  the  different  sections  of  die  data  were  determined  by  inspection,  they  were 
chosen  on  die  basis  of  predetermined  and  well-defined  characteristics  of  the  data. 
Based  on  an  examination  of  the  NOBS  data,  the  time  series  was  subdivided  into  four 
general  conditions  (Figure  4.17), 
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Table  4JL  Ifamr  wiwhtkw  wfflcUnti  between  N(UkHi),  and  tht  log  of  several 
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I.  Decaying  and  growing  wind  speed  and  wave  amplitude. 

II.  Fully  developed  quasi-steady  wind  and  wave  condition. 

III.  Turning  winds. 

IV.  Decaying  wind  and  wave  condition. 

Far  our  purposes,  we  considered  portion  II  of  the  record  as  steady,  and  portions  I,  III 
and  IV  unsteady.  We  then  computed  die  correlation  coefficients  for  the  steady,  unsteady 
and  die  full  time  sales  records  for  die  wind  and  wave  variables,  and  N(4.3  kHz).  The 
results  are  summarized  in  Table  45. 

Unlike  Penhallow  &  Dietz,  our  results  show  that  the  correlation  coefficients  for 
unsteady  data  is  higher  than  those  for  steady  data.  This  is  observed  in  both  the  wind 
and  wave  data.  The  difference  between  the  correlation  coefficient  of  the  full  data  set 
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Julian  day 

Figure  4.17.  Hie  wind  and  wave  time  series  during  NOBS  can  be  subdivided  into 
steady  and  unsteady  conditions.  The  data  in  section  II  were  considered  steady  while 
those  in  sections  I,  III  and  IV  woe  considered  unsteady. 
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and  the  unsteady  data  is  negligible  The  reduction  in  the  correlation  coefficient  (hiring 
steady  conditions  is  more  significant  for  wave  parameters. 

The  plots  ki  Figures  4.18  and  4.19  show  that  although  the  scatter  in  the  data  is 
comparable  for  the  steady  and  unsteady  cases,  the  range  of  the  unsteady  data  is  bigger 
and  it  is  comparable  to  die  full  data  range  Consequently,  die  correlation  coefficient  of 
die  unsteady  data  is  comparable  to  the  full  data  and  is  significantly  larger  than  that  of 
die  steady  data.  It  is  not  dear  why  the  PenhaDowAc  Dietz  (1964)  results  are  different 
from  ours  since  they  did  not  provide  enough  details  in  their  paper  for  us  to  examine 
their  data.  However,  given  the  dosenesa  of  the  correlation  coefficients  between  the 
unsteady  and  full  cases,  we  bebeve  our  results  are  robust  and  the  conclusions  in  this 
work  are  not  sensitive  to  the  quantity  of  the  data  we  analyzed. 

4.6.4  Wove  growth  and  wove  age. 

We  can  relate  die  wind  speed  Ujq  to  the  wave  energy  by  adopting  Kitaigorodskii's 
(1962)  non-dimensionalization  of  the  wave  and  wind  variables.  This  approach,  which 
was  used  by  Hasselmann  et  al  (1973)  to  analyze  fetch-limited  wave  growth  during  the 
Joint  North  Sea  Wave  Program  (JONSWAP),  relates  die  non-dimensional  wave 
variance 

(4.14) 

Ujo 

with  die  wave  age  (Cp/Ujo).  The  parameter  e  is  the  variance  of  the  surface  wave  field. 
Kitaigorodskii's  (1962)  non-dimensionalization  was  also  used  in  some  of  die  more 
recent  investigations  [Donelan,  Hamilton  &  Hui  (1985),  Dobson,  Perrie  &  Toulany 


unsteady  steady 


0.1  1  0.1  1 
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log$(0.75Hz)  log$(0.75Hz) 


Hgnn  4.18.  Plot  of  N(4.3  kHz)  versus  selected  wind  and  wave  parameters.  The  plots 
on  die  left  are  for  data  during  unsteady  conditions.  The  plots  on  the  right  are  for  data 
during  steady  conditions. 
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Figure  4.19.  Plot  of  U^o  versus  selected  wave  parameters.  The  plots  on  the  left  are  for 
data  during  unsteady  conditions.  The  plots  on  the  right  are  for  data  during  steady 
conditions. 
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(1969),  Donrian  et  al  (1992)]  on  fetch  limited  wave  growth  in  both  lake  and  ocean 
environments* 


Our  wave  measurements  are  compiksicd  by  the  presence  of  swell  whose 
characteristics  are  unrelated  to  the  local  wind  condition*.  Dobson  et  el  (1989),  whose 
measurements  olf  the  coest  of  Nova  Scotia  during  the  Canadian  Atlantic  Storms 
Program  (CASP)  were  eleo  contaminated  by  swell,  suggested  that  the  frequency 
separating  swell  and  the  wind  sea  portion  of  the  spectrum  f,  is  0.03  Hz  less  the  wind 
sea  peak  frequency.  The  total  energy  of  the  wind  sea  spectrum  is  therefore 


e* 


pCOd*. 


(415) 


Their  approach  is  similar  to  jut  method  of  separating  the  wind  wave  spectrum  from 
the  swell  components  (Section  4.3).  However,  we  defined  to  be  the  frequency  of  the 
minimum  value  of  4»(0  between  the  swell  and  wind  sea  peaks  and  its  value  was 
obtained  by  inspection. 


figure  420  shows  a  well  correlated  log-log  plot  of  e  and  the  inverse  wave  age 
computed  using  the  NOBS  data.  The  correlation  coefficient  between  the  two 
parameters  is  0.97.  Using  linear  orthogonal  regression,  we  found  that  the  relationship 
is  best  described  by  the  expression 


146 


Table  4.6  compare*  our  results  with  the  expersnenls  mentioned  earlier  in  this  section. 
Our  mute  are  dote  to  those  of  the  previous  experiments,  particularly  to  that  of 
Damian  el  al  (1992).  The  JONSWAP  result  did  not  separate  wind  seas  from  swell.  The 
agreement  between  our  result*  and  the  results  of  the  experiments  given  in  Table  4.6  is 
psrtkulariy  remarkable  eince  the  other  experiments  were  made  under  fetch-limited 
conditions.  Our  results  suggest  that  within  experimental  error,  the  empirical 
relationship  between  the  wind  and  wave  parameters  expressed  by  Equation  4.16  is 
valid  over  a  wide  range  of  fetch  lengths. 

4.6.5  Wave  spectrum  dtmr*cterittict 

Using  previous  observational  evidence,  Phillips  (1985)  suggested  that  the  portion  of 
wave  height  spectrum  Q(a>)  above  the  spectral  peak  can  be  described  by  the  equation 

4>(cd)  -  auLgto-4  (4.17) 


Table  4.6.  Coefficients  of  die  non-dimensional  wind  and  wave  growth  expression. 
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20 

where  a  is  the  Tabs  coefficient  Following  Phillips,  we  nondimensionaUzed  $(<o)  such 
that  for  an  ideal  spectrum  of  the  form  of  Equation  4.17, 


(4.18) 


The  term  inside  die  parenthesis  is  the  nondimensionalized  radian  frequency  to.  The 
plot  of  the  log  of  equation  4.18  and  log  (u*co/g)  should  have  a  straight  line  with 
slope*  -4.  Figure  4.21a  shows  the  12  spectra  in  Figure  4.2  which  were  taken  from  one 
hour  samples  of  the  wave  gauge  record  beginning  at  0000  UT  each  day  during  the 
observation  period.  The  solid  line  is  for  the  idealized  4(a>)  whose  a  *  0.11  (Equation 
4.8).  The  plots  show  a  dose  agreement  between  Equation  4.17  and  the  measured  Q(a>) 
in  the  frequency  range  between  a>p  and  approximately  Z5cop.  For  co  >  Up,  the  slope  of 
foe  wave  spectrum  appears  to  be  closer  to  <S>(g>)  ~  co*7.  Figure  4.21b  shows  typical 
wave  spectra  under  folly  developed  conditions  (Julian  days  273-278).  For  to  >  ay  the 
plots  collapse  into  one  spectrum. 

Phillips  (1985)  also  predicted  that  foe  upper  limit  of  foe  or4  slope  is  determined  by  foe 
presence  of  wind  drift  c^-u*  which  would  inhibit  the  creation  of  waves  frequencies 
whose  phase  speed  is  comparable  to  that  value  (see  Section  2.2).  He  proposed  that  the 
wavenumber  of  this  upper  limit  is 


*1  =  ©? /g  =  rg/u2 


(4.19) 


where  r  is  of  order  one.  He  subsequently  estimated  that  r  *  0.16  for  typical  wind  wave 
conditions.  Rearranging  Equation  4.19  gives 
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which  is  approximately  3  times  the  value  suggested  by  Figure  4.21.  This  suggests  that 
die  formation  of  high  frequency  wind  waves  is  inhibited  by  same  mechanism  whose 
phase  speed  is  lower  than  u».  It  is  currently  not  known  what  this  mechanism  is.  It  is 
also  possible  that  Keulegan's  (1951)  measurements  of  die  wind  drift  in  a  closed 
channel  [which  Phillips  (1985)  used]  overestimates  die  value  of  c^  in  the  open  ocean. 


S  SPECTRAL  DISSIPATION  ESTIMATES 


in  Chapter  2,  we  described  the  basis  of  the  Phillips  (1985)  and  the  KHH  surface  wave 
dissipation  models.  We  showed  (hat  die  two  models  differ  in  form,  in  their 
dependence  on  the  wave  spectrum,  and  in  die  method  used  to  evaluate  their 
coefficients.  While  Phillips  (1985)  is  based  on  die  assumption  of  a  spatial  and  temporal 
equilibrium  balance  among  die  source  terms  in  die  energy  transfer  equation,  KHH 
tuned  the  dissipation  parameters  with  die  wind  input  and  nonlinear  flux  source 
functions  to  reproduce  die  behavior  of  fetch-limited  wind  wave  growth.  We  also 
showed  that  we  can  derive  ample,  explicit  dissipation  equations  by  making  general 
assumptions  on  die  characteristics  of  the  wave  spectrum  and  substituting  die  wave 
spectrum  equation  into  the  full  spectral  dissipation  expressions. 

Table  5.1  summarizes  the  four  spectral  dissipation  estimates  discussed  in  Chapter  2. 

In  this  work,  d(a>)  denotes  die  spectral  dissipation  computed  from  die  measured  wave 
height  spectrum  while  e(co)  denotes  die  spectral  dissipation  computed  from  wind  and 
integral  wave  parameters  (reduced  estimate).  The  total  dissipation  estimates  are 
therefore 
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(5.1) 


In  the  subsequent  sections  of  this  chapter,  we  will  examine  the  spectral  behavior  of 
these  model  and  compare  these  estimates  to  ambient  noise  measurements. 
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TaMa  SJL  Dissipation  estimates  (see  Chapter  2  foe  lull  details.) 


1.  spectral  dissipation  equations 


<V(“)-'SP**'W<") 

eP(<»)  *  4tP3I(3p)pw  u2  oo_l 
dH(a>)-pwgcH((a/5)2(a/otm)2SwdKci>) 

cH(a>)*0.89pwgcw(a/o^)2 


2.  total  dissipation  equations 
Dp  m  lOptyw 8  y  fm1Iq)3(a))dto 

16CK P)f  J., 

Ep  «2T(P3I(3p)pwu2li^r(cp/u,f  j 
Dh*Pw« 

EH-1.33pwgcw(a/aPM)2  afo. 


5.1  Experimental  results. 

Although  tiie  four  models  (Dp,  Ep,  Dh»  Eh)  are  based  on  different  assumptions 
regarding  the  spectral  description  of  dissipation,  the  form  of  the  wave  height 
spectrum,  and  the  existence  of  an  equilibrium  wave  spectrum,  these  equations  give 
comparable  estimates  of  the  total  dissipation  for  most  of  the  measured  wind  and  wave 
condition.  This  is  observed  in  Figure  5.1  which  compares  the  time  series  of  the  four 
dissipation  estimates  during  the  NOBS  observation  period. 
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Figure  5.1*  shows  very  good  agreement  between  the  full  and  reduced  KHH 
dissipation  estimates  Dj*  and  Eh-  The  two  time  series  are  virtually  indistinguishable 
through  moat  of  the  observation  period.  Figure  5.1b  shows  a  comparison  between  the 
time  series  of  Dp  and  Ep.  Here  we  see  that  while  there  is  good  agreement  between 
these  two  estimates  over  most  of  the  data,  there  can  be  significant  differences  between 
them  especially  (hiring  periods  of  wind  and  wave  growth  and  decay.  These 
differences  are  explained  further  in  Section  5.1.2.  Figure  5.1c  compares  die  KHH  and 
Phillips  dissipation  estimated  from  the  measured  spectra  Dp  and  Dh.  In  this  figure, 
we  see  that  Dp  is  consistently  higher  by  a  small  factor  than  Djp  This  difference  is 
nearly  constant  This  is  because  the  magnitude  of  the  coefficient  of  Dp  is  slightly 
larger  than  that  of  Dh-  Note  however  that  this  difference  is  small  compared  to  the 
dynamic  range  of  die  data. 

These  time  series  comparisons  (Figures  5.1a-d)  show  that  the  different  estimates  agree 
well  under  steady  wind  and  wave  conditions  but  may  vary  significantly  during 
periods  of  growth  or  decay.  This  result  is  not  surprising  since  die  coefficients  of  the 
different  dissipation  expressions  were  chosen  such  that  the  dissipation,  the  wind 
input  and  the  nonlinear  flux  expressions  are  of  comparable  order.  Note  that  even 
though  the  influence  of  wind  speed  on  die  Snyder  et  al  (1981)  wind  input  formula 
(used  by  KHH]  (Equation  1.14)  is  sHghfly  different  from  the  Plant  (1982)  wind  input 
formula  [used  by  Phillips  (1985)]  (Equation  1.15),  these  two  formulas  give  values  that 
are  comparable  to  each  other  in  the  range  f  <  1  Hz  (Figure  1.9).  Hence,  the  Phillips 
(1985)  and  KHH  dissipation  estimates  are  consequently  of  comparable  order. 
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Figure  5.1.  Comparison  of  the  time  series  of  the  log  of  the  total  dissipation  [kg/s3] 
computed  using  die  formulas  given  in  Table  5.1.  a)  Djj  (solid  line)  and  (short 
dash),  b)  Dp  (solid  line)  and  Ep  ( short  dash). 
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Figure  5.1.  Continued  from  last  page,  c)  Dp  (solid  line)  and  (short  dash),  d)  EP 

(solid  line)  and  Eh  ( short  dash). 
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e)  273  f  (Hz) 


f)  274  f  (Hz) 


Figure  5^L  Hots  showing  the  spectral  dependence  of  the  KHH  dissipation  d^(f)  and 
^(f)  [kg/s3  Hz].  These  spectral  estimates  were  computed  from  the  wave  spectra 
shown  in  Figure  4.2.  Numbers  below  the  plots  indicate  the  Julian  day  from  which  the 
plots  were  taken. 


in  the  discussion  in  Section  2.4,  we  showed  how  Eh  was  derived  by  substituting  an 
irissBrert  or*  ipcctrum  into  d^.  The  mean  frequency  T  and  the  variance  e2  of  the 
ktsaBnrd  cpectram  are  equal  to  thoae  of  die  wind  wave  portion  of  the  measured  wave 
spectrum,  bi  Figure  we  show  the  spectral  dependence  of  the  dissipation  d^f)  and 

CH<f)  for  the  twelve  wave  height  spectra  shown  in  Figure  4.2.  The  solid  line  plot 
indicates  dissipation  computed  usfrig  the  measured  wave  spectrum  while  the  dashed 
line  indicates  the  spectral  dissipation  using  the  idealized  wave  height  spectrum 

♦*(«)■  3fl£o>Jar*.  (2.24) 

The  spectral  estimates  are  very  similar  both  in  magnitude  and  spectral  slope  in  much 
of  the  wind  wave  region  (f  <  1Hz).  Above  1  Hz,  the  slopes  of  the  measured  wave  and 
dissipation  spectra  are  steeper.  While  die  influence  of  die  discrepancy  in  the  higher 
frequencies  on  die  total  dissipation  is  small,  it  does  lead  to  an  underestimate  of  the 
peak  frequency  flOp  since  we  assumed  that  the  spectral  slope  is  constant  throughout 
when  in  fact,  the  measured  wave  spectrum  has  an  or4  slope  from  the  spectral  peak  to 
approximately  1  Hz  and  a  steeper  or7  slope  above  1  Hz- 

Much  of  the  difference  between  and  ^  arises  because  we  chose  not  to  model  the 

dissipation  below  the  wind  wave  spectral  peak.  We  see  from  Figure  5.2  that  dH(f) 

gives  values  of  dissipation  due  to  swell  that  are  significant  Although  the  agreement 
in  the  wind  wave  spectral  range  is  very  good,  die  reduced  formula  EH  completely 
neglects  die  dissipation  in  die  swell  region  computed  in  D^-  While  the  total  energy  in 
die  swell  region  is  a  significant  proportion  of  $(co)  and  in  some  cases  contains  more 
energy  than  wind  waves  (see  sample  spectra  in  Figure  4.2),  its  contribution  to  the  total 
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dfcaipatkm  is  comparably  less  because  the  contribution  of  die  wave  frequency 
component*  i*  weighted  by  a?  (see  Equation*  2.18  or  2.19)  which  favors  the  higher 


frequencies. 

It  is  also  important  to  point  out  that  KHH  implicitly  assumed  a  single  peak  spectrum 
in  deriving  the  wave  evolution  equations.  Their  approach  did  not  factor  in  the 
influence  of  swell  on  the  dissipation  equation.  Since  their  model  assumes  that 
dissipation  in  deep  water  is  due  to  wave  breaking  and  since  swell  does  not  break 
(although  sweD  may  induce  breaking  of  wind  waves),  extending  the  KHH-based 
formula*  to  indude  dissipation  in  the  swell  frequency  range  may  overestimate  the 
total  dissipation.  Nonetheless,  the  two  time  series  in  Figure  5.1a  show  that  this 
difference  is  small  compared  to  the  dynamic  range  of  the  total  dissipation. 

5.13  Dp  and  Ep. 

dearly,  the  difference  in  the  Dp  and  %  estimates  is  more  substantial  than  that 
between  Djj  and  Eh-  This  is  because  Bp  uses  an  equilibrium  spectrum  (Equation  2.8) 
assumption  to  model  Are  measured  wave  spectrum  while  the  Eh  model  spectrum 
«m(cd)  is  computed  from  the  integral  characteristics  of  the  measured  spectrum  itself. 
The  coefficients  of  were  chosen  such  that  the  measured  and  the  idealized  spectra 
have  die  same  energy.  Hence,  the  validity  of  our  wave  spectrum  assumptions  in 
computing  Eh  from  the  integral  wind  wave  spectrum  is  based  only  on  whether  or  not 
the  slope  of  the  wind  wave  portion  of  the  measured  spectrum  <J»  is  close  to  ay*.  The 
variance  a2  and  mean  frequency  f  of  are  identical  to  those  of  O.  On  the  other 

hand,  using  $£  (Equation  2.8)  to  compute  Ep  is  based  primarily  cm  the  validity  of  two 
assumptions:  that  the  spectral  slope  of  die  equilibrium  spectrum  is  indeed  or4. 


160 


a) 


Julian  day 
b) 


Figure  5.4.  a)  Plot  showing  the  time  series  of  the  ^  and  $£  coefficients.  0.89a*co£ 
(solid  line)  and  a  tug  (short  dash)  [mVs3].  b)  Plot  showing  the  relative  magnitude  of 

the  difference  between  0.89a£a£  and  aiug. 
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and  that  the  coefficient  of  or*  in  *  is  cojug.  The  second  assumption  is  also  equivalent 
to  die  assumption  that  the  wind  input,  dissipation  and  nonlinear  spectral  flux  are  in 
equilibrium  balance  under  all  wind  and  wave  conditions. 

In  Figure  4.21,  we  showed  a  normalized  plot  comparing  the  twelve  measured  spectra 
in  Figure  42  with  Phillip*’  (1985)  equilibrium  spectrum  (Equation  2.8).  In  that 
figure,  we  saw  that  the  magnitude  and  slope  of  the  measured  wave  spectrum  Q(co)  are 
comparable  to  those  of  $g(oo)  between  die  wind  wave  peak  and  approximately  1  Hz. 
We  can  explore  die  relationship  further  here.  The  equilibrium  spectrum  is 

$E  *aiugaT4  (2.8) 

where  a  -  0.11  (see  abo  Section  2.2).  In  Section  2.4,  we  showed  that  die  idealized 
spectrum 


*K(a>)  *  3a^D^Q>~*  (2.24) 

gives  a  spectrum  whose  integral  characteristics  are  equal  to  those  of  the  measured 
wave  spectrum.  We  modify  Equation  2.24  by  substituting  Equation  2.5  to  express  it  in 
terms  of  mw  and  get 


«>M((0)»0.89a^fl)'4.  (5.1) 

Note  that  Figure  5.2  shows  that  this  equation  does  well  in  modeling  the  measured 
wind  wave  spectrum. 
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The  variable*  in  the  coefficient  of  Equation  5.1  arc  the  sea  surface  variance  and  the 
ween  wave  frequency  5*  while  in  the  coefficient  of  or4  in  Equation  2.8,  u*  is  a 
variable  while  both  a  and  g  are  constants.  We  can  examine  how  effective  <X>£  is  in 
modeling  the  measured  wave  spectrum  by  comparing  the  coefficients  of  and 
Figure  5.3  shows  a  plot  of  0. 89<£ia*  versus  a  lug  using  die  data  from  the  NOBS 
observation  period.  The  figure  shows  that  the  two  parameters  are  well  correlated 
(conriation  coefficient  ■  0.92).  The  dashed  line  indicates  a  plot  of  slope  « 1.  The 
equation  of  die  orthogonal  linear  regression  of  die  data  is 

low*]-  [0.89<^]+0.0017, 

[on*]-  10.89<#£].  (5.2) 


This  shows  that  $£is,  on  average,  a  reasonable  approximation  to  d>.  It  is  important  to 
point  out  however  that  die  scatter  of  the  data  around  Equation  5.2  (Figure  5.3) 
suggests  that  we  incorporate  a  small  error  when  we  use  this  approximation.  The 
nature  of  die  difference  is  seen  in  Figure  5.4  which  shows  a  plot  of  the  time  series  of 
0.89<£a£  and  [a vug].  Although  the  plot  shows  good  agreement  between  the  two 
parameters,  there  dearly  are  differences  between  the  two.  This  is  demonstrated  in 
Figure  5.4b  which  plots  the  difference  between  0.89<£m*  and  [otu*g]  normalized  by 
0.89<£o£.  From  this  plot,  we  see  that  [au*g]  significantly  overestimates  0.89<£c5* 
during  die  early  stages  of  wind  and  wave  growth  (Julian  day  271).  When  the  wind  is 
decaying  (Julian  day  279-280),  [aiug]  significantly  underestimates  0.89<£io^.  The  time 
series  also  show  that  even  under  fully  developed  conditions,  the  wave  field  0.89a£co* 
may  not  respond  quickly  to  rapid  fluctuations  in  the  wind  speed  [au*g] . 


e)  273  f  (Hz) 


f)  274  f  (Hz) 


Figure  5.5.  Plots  showing  the  spectral  dependence  of  die  Phillips  dissipation  dp(f)  and 
Ep(f)  [kg/s3  Hz].  The  spectral  estimates  were  computed  from  the  wave  spectra  shown 
in  Figure  4.2.  Numbers  below  the  plots  indicate  the  Julian  day  from  which  the  plots 
were  taken. 
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The  plot*  in  Figure  5.4  suggest  that  under  fully  developed  or  slowly  varying  wind 
and  wave  conditions,  Phillips’  (1985)  equilibrium  spectrum  (Equation  2.8)  adequately 
describes  the  wave  field.  Under  conditions  of  rapid  wind  change,  the  assumption  of 
local  spatial  and  temporal  equilibrium  may  not  hold  and  the  balance  among  the 
source  terms  (Equation  2.6)  may  not  be  strictly  valid. 

hi  Figure  5.5,  logarithmic  scale  plots  of  the  spectral  dependence  of  the  Phillips  (1985) 
dissipation  are  shown.  Similar  to  Figure  5.2,  the  plots  compare  the  dissipation 
estimates  computed  from  the  measured  spectra  0  and  the  model  spectra  Og.  Here  we 
note  that  the  difference  between  die  two  spectral  estimates  is  larger  in  Figure  5.5  than 
in  Figure  5.2.  This  could  appear  surprising  given  drat  Figures  4.21  and  5.3  show  that 
0  and  0g  agree  well. 


The  difference  is  mainly  due  to  die  cubic  dependence  of  dp  and  ep  cm  0  and  0g 
respectively.  This  makes  diem  more  sensitive  to  differences  between  die  values  of  the 
of  the  two  spectra.  For  example,  if  0(oo^  *  1.0  and  OgC^o)  =  1-1,  then  A0(co^)  =  0.1  or  a 
difference  of  10%.  However,  03(cqo)  =  1.0  and  0g3(ato)  =  1.33  gives  a  much  larger 
difference  of  33%.  Hence,  although  O(co)  and  0e(cd)  are  close  to  each  other  over  a 
significant  portion  of  die  wind  wave  frequency  range,  die  cubic  wave  spectrum 
dependence  magnifies  the  difference  between  dp  and  ep 

The  spectral  behavior  of  die  Phillips  (1985)  dissipation  also  explains  why  Ep  is  initially 
significantly  larger  than  Dp  during  wind  wave  growth.  The  directional  wave  spectra 
plots  (Figure  4.7  and  4.8)  suggest  that  during  the  initial  rapid  increase  in  wind  speed, 
die  higher  frequency  waves  energy  achieve  equilibrium  much  earlier  than  the  energy 
containing  portion  of  the  spectra.  The  development  of  the  energy  containing  portion 
lags  the  wind  growth.  Since  the  value  of  Dp  is  controlled  by  the  total  wave  energy 
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while  the  value  of  Bp  it  controlled  by  u»,  we  therefore  expect  that  Dp  will  lag  Ep  as 
weO. 

U  Relationship  to  ambient  noise. 

Laboratory  experiments  (eg.,  Melville  et  al  (1988),  Loewen  8c  Melville  (1991a)]  show 
that  file  sound  generated  by  breaking  waves  is  proportional  to  the  total  energy 
dissipated  by  the  event  over  a  limited  range  of  parameters.  Melville  (1992)  suggested 
that  this  proportionality  is  due  to  die  fact  that  the  dissipation  and  ambient  noise 
generation  are  coupled  by  the  air  entrainment  process.  The  correlation  between  sound 
generation  and  dissipation  by  breaking  waves  in  the  laboratory  suggests  that  sound 
may  be  used  to  scale  the  frequency  and  intensity  of  wave  breaking  in  the  field. 
Although  models  (eg..  Fanner  &  Vagle,  1988)  have  suggested  that  the  frequency  and 
intensity  of  wave  breaking  in  the  vicinity  of  a  hydrophone  is  proportional  to  the 
ambient  noise  level,  it  has  yet  to  be  demonstrated  by  field  measurements  that 
dissipation  by  breaking  waves  correlates  with  the  ambient  noise  in  the  field. 

Figure  5.6  shows  plots  of  Dj*  versus  the  ambient  noise  at  4.3, 8.0, 14.0  and  2.5-14  kHz 
(£).  Similar  plots  showing  die  relationship  between  ambient  noise  and  the  other 
dissipation  estimates  are  given  in  the  succeeding  figures:  Figure  5.7  (Eh),  Figure  5.8 
(Dp)  and  Figure  5.9  (Ep).  The  plots  show  a  good  correlation  between  ambient  noise 
and  the  dissipation  except  for  N(14  kHz).  As  we  pointed  out  earlier,  field 
measurements  [eg..  Farmer  &  Lemon  (1984),  Vagle  et  al  (1990)]  suggest  that  sound 
absorption  by  bubbles  generated  by  breaking  waves  can  significantly  alter  noise  levels 
above  8  kHz.  For  this  reason,  we  will  confine  our  discussion  in  this  section  to  the  three 
other  noise  frequencies  given  in  the  figures. 
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log(OH) 


log(DH) 


Figure  5.6.  Plots  of  the  KHH  dissipation  versus  ambient  noise  at  selected 
frequencies.  The  dashed  line  indicates  die  regression  computed  from  die  data  (see 
Table  5.1  for  die  coefficients). 
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Figure  5.8.  Plots  of  the  Phillips  (1985)  dissipation  Dp  versus  ambient  noise  at  selected 
frequencies.  The  dashed  line  indicates  the  regression  computed  from  the  data  (see 
Table  5.1  for  the  coefficients). 
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Figure  5.9.  Plots  of  the  Phillips  (1985)  dissipation  Ep  versus  ambient  noise  at  selected 
frequencies.  The  dashed  line  indicates  die  regression  computed  from  the  data  (see 
Table  5.1  for  the  coefficients). 
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Table  &2.  Coefficients  for  the  ambient  noise  and  dissipation  expicadoii  and  the 


Chapter  5 

We  also  see  from  Pigure  5.5  -  5.7  that  the  distribution  of  the  data  around  the 
regression  lines  appear  to  be  qualitatively  similar  for  and  Dp.  This  should  be 

evident  from  the  dissipation  time  series  given  in  Figures  5.1a  and  5.1c  In  Figure  5.1a, 
we  see  that  the  and  time  aeries  are  nearly  indistinguishable.  In  Figure  5.1c  the 

Dh  and  time  series  also  look  very  similar  with  the  difference  between  the  values  of 

log  Dp  and  tog  nearly  constant 

We  also  see  that  die  dissipation  data  where  tog  D  >  -1.5  show  considerably  less  scatter 
than  data  where  log  D  <  *1.5.  This  is  partly  due  to  the  distribution  of  the  available 
data.  In  Section  4.63,  we  saw  that  the  wind  and  wave  conditions  were  quasi-steady 
with  the  wind  speed  and  dissipation  constant  for  approximately  half  of  die 
observation  period.  As  a  result,  die  tow  N  data  is  sparse  and  die  data  distribution  is 
biased  towards  die  higher  winds  speeds  and  dissipation  values. 

It  is  important  to  point  out  that  except  for  Bp,  the  dissipation  estimates  were  all 
computed  directly  from  wave  field  parameters.  was  computed  from  the  variance 

of  die  wind  wave  field  aj*  and  die  mean  frequency  while  Djj  and  Dp  were 

computed  from  the  measured  spectra.  In  Section  12,  we  described  how  previous  field 
measurements  have  found  that  die  ambient  noise  N  and  die  log  of  the  root  mean 
square  wave  amplitude  a  are  poorly  correlated.  The  correlation  coefficient  between  die 
two  parameters  for  die  NOBS  data  is  0.67  which  is  comparable  to  die  value  obtained 
by  Perrone  (1969)  and  by  Farmer  &  Lemon  (1984). 

Our  analysis  in  Section  4.6.2  showed  that  by  filtering  out  swell  and  comparing  N  with 
die  RMS  wave  height  of  the  wind  wave  components,  we  can  improve  the  relationship 
between  the  two  (Figure  4.13).  By  neglecting  the  contribution  of  swell,  we  increased 
the  correlation  coefficient  between  N  and  log  aw  to  0.87.  A  similar  value  for  the 


174 


CktpttrS 

conriatkm  coefficient  was  obtained  between  N  and  log  s  (the  RMS  wave  slope).  While 
correlation  coefficients  for  N  vs.  log  a*,  is  significantly  higher  than  the  values  for  the 
correlation  coefficient  for  N  vs.  log  a,  it  is  still  significantly  less  than  the  correlation 
coefficient  of  the  N  vs.  log  U  relationship  (correlation  coefficient  ■  0.94). 

Section  4.6.3  showed  that  comparisons  between  wind  and  wave  parameters  and  N 
during  steady  conditions  are  not  meaningful  since  the  dynamic  range  of  die  values 
during  this  period  is  small.  Hence,  we  omit  comparisons  between  N  and  log 
dissipation  during  steady  conditions.  Section  463  also  showed  that  the  results  for  the 
unsteady  cases  are  similar  to  those  for  the  combined  steady  and  unsteady  cases.  Table 
53  gives  the  correlation  coefficients  between  N  and  the  different  dissipation  estimates. 
We  see  from  this  table  that  the  relationship  between  log  dissipation  and  N  give 
correlation  coefficients  (0.93-0.95)  that  are  comparable  to  that  between  log  U  and  N. 

Table  5-2  also  compiles  die  different  regression  equations  for  the  various  dissipation 
versus  ambient  noise  plots.  For  purposes  of  this  work,  we  adopted  the  form  of  the 
WOT  AN  equation  to  relate  ffie  dissipation  D  to  ambient  noise. 

10gD«  m(f)N(f)f  n(f).  (5.3) 


The  coefficients  m  and  n  in  Table  5.2  were  computed  using  the  orthogonal  linear 
regression  formula  (Casdla,  1990,  p.  586).  The  regression  lines  are  shown  as  dashed 
lines  in  each  plot  in  Figures  5.6  -  5.9. 

Equation  5.3  is  equivalent  to  a  power  law  expression  between  the  sound  pressure  and 
surface  wave  dissipation.  We  can  therefore  express  it  in  terms  of  the  equation 
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where  m  ■  l(W1Qna  and  n*l/(10m).  Our  computations  suggest  that  the  power  law 
dependence  between  the  ambient  noise  pressure  level  and  dissipation  is 

p*~I f***  (5.5) 

lor  the  estimates  computed  from  wave  field  characteristics  (Du,  %  and  Dp),  and  a 
slightly  steeper  power  law  dependence 


p’-DP* 


(5.6) 


when  die  dissipation  is  computed  from  u*  (Bp). 

Despite  these  differences,  our  results  show  that  the  different  models  give  dissipation 
estimates  that  are  comparable.  They  also  show  that  dissipation  correlates  well  with 
ambient  noise  in  the  field.  This  is  consistent  with  die  results  of  die  laboratory 
experiments  we  reviewed  in  Section  1.4  which  suggest  Out  the  wave  energy 
dissipated  and  the  acoustic  energy  generated  by  a  breaking  wave  are  both 
proportional  to  the  strength  of  breaking.  Since  surface  wave  dissipation  and  die 
generation  of  underwater  sound  are  coupled  together  by  the  air  entrainment  process 
associated  with  breaking  waves,  then  die  two  correlated  well  in  the  laboratory.  Our 
results  show  that  die  correlation  between  die  two  are  observed  in  the  field  as  well 


t  DISCUSSION  AND  CONCLUSIONS 


In  the  review  ki  Chapter  1,  we  discussed  how  several  decades  of  research  have  shown 
that  the  wind  speed  shows  a  very  strong  correlation  with  ambient  noise  N  in  the 
Knudsen  range.  The  research  in  the  last  decade  or  so  have  taken  advantage  of  this  in 
developing  Weather  Obeervetkms  Through  Ambient  Noise  (WOTAN),  a  method  which 
uses  die  lewd  of  ambient  noiae  at  preselected  frequencies  to  estimate  the  wind  speed 
U  using  the  linear  relationship  between  N  and  log  U. 

In  contrast,  research  relating  the  ambient  noise  level  to  parameters  of  the  wave  field 
have  been  few  and  largely  unsuccessful  The  available  measurements  have  shown 
that  N  and  the  significant  wave  height  (or  the  root  mean  square  wave  amplitude  a) 
are  poorly  correlated.  No  attempt  has  been  made  to  explore  this  subject  further  nor 
have  attempts  been  made  to  establish  a  mathematical  relationship  between  die  two 
parameters. 


Measurement  have  shown  that  breaking  waves  in  the  field  generate  sound  as  low  as 
30  Hz  and  up  to  20  kHz.  Although  acoustic  field  measurements  below  500  Hz  have 
been  known  to  correlate  with  wind  speed,  the  influence  of  shipping,  industrial  and 
other  man-made  noise  can  contaminate  measurements  at  these  frequencies. 
Laboratory  measurements  have  shown  that  die  major  sources  of  sound  from  breaking 
waves  are  due  to  bubble  oscillations  during  creation,  coalescence  and  breakup.  These 
mechanisms  dominate  die  sound  frequencies  above  1  kHz  and  perhaps  as  low  as  a 
few  hundred  Hertz.  At  the  lowest  frequencies  of  wind  dependent  noise  of  the  order  of 
tens  of  Hertz,  there  is  evidence  that  the  sound  is  due  to  the  collective  oscillations  of 
bubble  clouds. 
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If  it  is  indeed  wave  breaking  that  is  responsible  for  the  noise  in  the  wind  dependent 
frequency  tinge,  it  is  then  surprising  that  the  characteristics  of  the  wave  field  and  the 
ambient  underwater  sound  have  appeared  to  be  poorly  correlated.  Although  wave 
breaking  is  ultimately  driven  by  the  wind,  the  wave  breaking  process  is  a  direct 
product  of  die  characteristics  of  the  wave  field.  This  apparent  contradiction  has  yet  to 
be  resolved. 

We  believe  that  die  key  to  understanding  this  relationship  is  to  relate  N  to  wave 
parameters  that  are  either  directly  responsible  for  the  occurrence  of  wave  breaking  or 
are  direct  manifestations  of  die  wave  breaking  process.  The  various  laboratory 
measurements  of  wave  breaking  by  Melville  and  his  coworkers  have  provided  us 
with  some  dues. 


Laboratory  experiments  of  mechanically  generated  breaking  waves  suggest  that  a 
packet  consisting  of  focused  short  and  long  waves  can  produce  breaking  waves 
whose  gross  characteristics  can  be  consistently  reproduced.  The  packets  can  be 
characterized  by  a  slope  parameter  sp  whose  value  is  a  measure  of  die  aggregate 


steepness  of  the  components  of  the  wave  packet  These  experiments  found  that  the 


magnitude  of  Sp  correlates  with  die  magnitude  of  several  important  characteristics  of 


breaking  wave  including,  loss  of  excess  momentum  flux,  mixing,  initial  volume  of  air 
entrained,  wave  energy  dissipation  and  the  acoustic  energy  generated  by  the  event. 


In  particular,  we  focus  on  the  recent  studies  by  Melville  et  al  (1988)  and  Loewen  & 
Melville  (1991a)  which  show  that  the  energy  dissipated  by  a  breaking  wave  is 
proportional  to  the  volume  of  air  entrained  and  die  sound  generated  by  the  event.  The 
results  of  the  work  of  Lamarre  &  Melville  (1991)  on  the  evolution  of  the  bubble  cloud 
generated  by  the  breaking  wave  suggest  that  a  large  fraction  (up  to  50%)  of  the  energy 


dissipated  is  due  to  the  work  done  against  buoyancy  in  entraining  die  bubble  plume 
At  die  same  time/  creating  air  bubbles,  forcing  them  and  die  bubble  cloud  itself  leads 
to  the  conversion  of  mechanical  energy  into  acoustic  energy.  This  process  is  clearly 
not  very  efficient  since  only  a  very  small  fraction  of  die  order  of  10**  (Loewen  & 
Melville,  1991at)  of  the  mechanical  energy  is  actually  converted  to  sound. 
Nevertheless,  both  dissipation  and  sound  generation  in  breaking  waves  are  coupled 
together  by  die  air  entrainment  process  and  we  believe  that  this  is  the  reason  why  the 
two  correlate  well. 

The  laboratory  experiments  clearly  suggest  that  level  of  dissipation  could  be  used  as 
an  effective  variable  in  parameterizing  the  occurrence  of  wave  breaking  in  die  field.  In 
principle,  one  can  compute  the  total  dissipation  in  die  wave  field  by  summing  up  the 
total  energy  dissipated  by  all  the  individual  events.  In  practice,  this  is  difficult  o 
undertake  since  not  only  is  it  hard  to  estimate  how  much  energy  each  individual 
event  dissipates,  it  is  also  difficult  to  came  up  with  a  reliable  and  effective  method  for 
performing  areal  surveys  of  the  occurrence  and  intensity  of  breaking  waves  in  die 
field  using  commonly  available  instruments  and  techniques. 

Alternatively,  we  can  indirectly  estimate  the  total  dissipation  by  examining  how 
energy  is  added,  dissipated  and  redistributed  in  the  wave  field.  T  process  is 
described  mathematically  by  the  energy  transfer  equation  (Equation  1.11)  which 
expresses  die  changes  in  die  wave  energy  spectrum  as  a  sum  of  the  wind  energy 
input  w(k),  the  nonlinear  flux  of  energy  from  the  higher  frequencies  to  die  lower 
frequencies  Vk-T(k),  and  the  dissipation  of  the  wave  energy  d(k) 

The  characteristics  of  the  spectral  dissipation  d(k)  are  at  present  not  very  well 
understood.  In  Chapter  2,  we  described  two  theories  on  how  d(k)  can  be  computed 
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from  the  wave  spectrum.  Using  heuristic  arguments  on  the  influence  of  nonlinear  Lut 
localized  wave  breaking  on  the  wave  field,  Hasaebnann  (1974),  proposed  that 

dHOO-tt^OO. 

The  coefficient  of  this  relationship  was  later  obtained  by  Kotnen  et  al  (1984)  from 
numerical  simulations  of  the  equilibrium  balance  between  d(k),  w(k)  and  Vk-T(k). 
Although  subsequent  numerical  wave  evolution  models  (eg.,  SWAMP  (1985), 
WAMDIG  (1988)]  have  attainted  to  incorporate  more  parameters  and  to  enhance  die 
stability  of  dfj,  the  bask  assumptions  Hasaebnann  (1974)  made  in  deriving  die  form  of 
df|  are  still  used  in  many  of  the  current  wave  evolution  models. 

To  date  the  most  rational  approach  to  the  estimate  of  spectral  dissipation  is  due  to 
Phillips  (1985)  who  derived  an  analytical  expression  for  dissipation  based  essentially 
an  two  major  assumptions  on  die  characteristics  of  die  equilibrium  range  above  the 
wind  wave  spectral  peak:  (a)  that  the  magnitude  of  d(k),  w(k)  and  Vk-T(k)  are 
comparable,  and  (b)  that  compared  to  die  time  and  spatial  scales  of  evolution  of  the 
wave  spectrum,  die  equilibrium  range  wave  components  are  quasi-steady  [or  that  the 
sum  of  d(k),  w(k)  and  Vk-T(k)  is  zero].  From  these  assumptions,  Phillips  (1985) 
derived  an  expression  for  d(k)  and  showed  that  in  the  equilibrium  range,  the  wind 
wave  spectrum 
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While  dH  and  dp  are  both  based  cm  assumptions  of  equilibrium  among  the  source 
terms,  their  major  differences  arise  from  the  different  assumptions  on  the  influence  of 
the  wave  spectrum  on  d(k)  and  the  differences  on  the  assumed  forms  of  the  other 
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source  terms  in  foe  energy  transfer  equation.  Hence,  foe  use  of  either  model  to 
estimate  dissipation  requires  that  we  first  examine  the  differences  between  the  two 
and  how  their  values  compare  using  real  wind  and  wave  data. 


In  Chapter  3,  we  described  an  experiment  in  the  open  ocean  130  km  off  the  coast  of 
Oregon  where  we  made  simultaneous  measurements  of  the  wind  speed,  die  ambient 
noise  and  the  wave  height  spectrum  using  an  array  of  wire  wave  gauges  from  die 
research  platform  Flip.  The  hydrophone  data  was  sampled  at  51.2  kHz  using  a 
directional  hydrophone.  The  power  spectrum  of  die  data  was  computed  in  real  time 
by  a  DSP  board  in  a  PC  and  stored  on  optical  disks.  By  processing  the  data  in  situ,  we 
significantly  reduced  foe  storage  and  post-experiment  processing  requirements  to 
manageable  levels.  All  other  instruments  used  in  our  experiment  were  sampled  by 
Robert  Pinkel  and  Jerome  Smith  of  foe  Upper  Ocean  Physics  Group,  Scripps 
Institution  of  Oceanography.  From  this  data,  directional  wave  spectra  were  computed 
using  the  maximum  likelihood  method  from  the  wave  height  measurements  of  the 
four  wire  wave  gauges.  The  time  series  of  the  root  mean  square  wave  slope  s  was  also 
computed. 

We  found  (Chapter  4)  that  due  to  the  various  noise  sources,  the  higher  frequency 
portion  of  tiie  ambient  noise  data  was  unusable  However,  most  of  the  acoustic  data  at 
frequencies  below  14  Hz  appear  to  be  reasonably  free  of  noise  once  the  sources  of 
noise  were  characterized  and  eliminated.  For  our  purposes,  the  comparisons  with 
ambient  noise  of  different  frequencies  with  wind  and  wave  parameters  give 
qualitatively  comparable  results  so  unless  indicated,  this  section  will  discuss  the 
N(4.3  kHz)  results.  In  general,  we  found  that  the  wind  versus  ambient  noise 
relationships  from  our  measurements  are  consistent  with  those  reported  earlier.  Our 
measurements  give  (Figures  4.9  and  4.10) 
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The  mare  novel  aspects  of  this  study  however  are  in  the  comparisons  between  the 
wave  field  and  die  ambient  nose.  As  we  mentioned  earlier,  we  found  it  surprising 
that  there  seemed  to  be  very  little  in  die  literature  that  sought  to  establish  and  define 
relationships  between  these  two.  The  few  experiments  that  were  reported  show  the 
root  mean  square  wave  amplitude  a  and  N  to  be  poorly  correlated. 

PenhaBow  &  Dietz  (1964)  suggested  that  die  changes  in  the  wind  speed,  a  and 
ambient  noise  are  better  correlated  under  steady  wind  conditions.  When  the  wind  was 
increasing  or  decreasing,  N  and  U  were  still  reasonably  well-correlated  but  N  and  a 
were  not  In  their  paper,  they  showed  a  correlation  coefficient  of  0.90  between  N  and  a 
under  steady  winds  and  0.42  when  the  wind  speed  is  increasing.  We  found  results 
(Section  4.6.3)  that  were  die  opposite  of  those  they  reported.  Our  results  showed  that 
the  data  correlated  better  during  unsteady  cases  (Figure  4.5)  than  during  steady  wind 
conditions.  With  our  data,  the  reason  is  dear  from  Figure  4.18.  In  that  figure,  we  see 
that  the  scatter  in  the  data  is  relatively  constant  under  both  categories  but  the  range  in 
Uio,*  and  N  values  is  larger  for  unsteady  conditions.  For  steady  winds,  the  data  is 
clustered  together  thus  making  the  data  appear  uncorrelated. 

It  is  difficult  to  make  meaningful  comparisons  between  the  Penhallow  &  Dietz  (1964) 
data  and  ours  since  they  showed  no  plots  of  die  data  in  their  paper.  However,  die  fact 
that  the  correlation  coefficients  for  the  unsteady  data  set  and  the  complete  data  set  do 
not  differ  significantly  (less  than  ±  0.02)  [see  Table  4.5]  suggests  that  our  analysis  of 
the  data  is  robust  and  is  not  sensitive  to  the  quantity  of  the  data  we  analyzed. 
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We  found  that  while  a  and  N  do  not  cordate  well,  we  can  improve  die  correlation 
between  the  two  if  we  filter  out  die  swell  components  from  die  wave  elevation  time 
aeries.  When  this  is  done,  die  relationship  between  N  and  a*  is  improved  significantly 
(Figure  413).  Our  results  show  that 


The  effect  of  swell  on  the  sea  surface  variance  a2  is  disproportionate  to  its  effect  on  the 
frequency  of  wave  breaking  and  consequently,  to  N.  Since  swell  is  generated  outside 
die  region,  its  properties  are  not  a  direct  result  of  die  local  wind-wave  conditions.  On 
acahn  day,  we  can  get  high  a  values  if  there  is  significant  swell.  By  filtering  this  out 
and  considering  only  the  wave  fidd  components  that  are  generated  locally,  then  the 
s^  versus  N  relation  is  improved. 

In  addition,  we  note  that  die  spectral  energy  in  the  wave  height  frequency  spectrum 
GKO  is  maximum  at  die  wind  wave  spectral  peak  and  decays  very  rapidly  at  higher 
frequencies  (frequency  slope  f4  up  to  approximately  1  Hz  and  f*7  above  that).  Hence, 
die  change  in  die  value  of  a,,  is  dominated  by  die  changes  in  die  energy  of  the  wave 
frequencies  near  die  spectral  peak.  However,  we  found  that,  as  the  Snyder  et  al  (1981) 
and  the  Plant  (1982)  wind  input  formulas  suggest,  the  higher  frequencies  are  more 
responsive  to  changes  in  the  wind  speed  and  direction  and  consequently,  to  changes 
inN. 

This  is  reflected  in  how  die  directional  wave  spectrum  changes.  In  Figures  4.7  and  4.8, 
we  saw  that  when  die  wind  began  to  pick  up  after  it  changed  direction  by  180*,  the 
wave  energy  in  the  new  direction  was  first  seen  at  the  higher  frequency  wave 
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components  (f  ~  1  Hz).  Eventually,  the  energy  in  the  lower  frequency  components 
began  increasing  until  the  directional  distribution  of  the  wave  energy  spectrum 
regained  its  original  form  before  the  change  in  the  wind  direction. 

In  Figure  114,  we  also  observed  that  although  die  relationship  between  die  log  of  the 
spectral  level  of  the  wind  wave  frequencies  above  the  peak  is  linear  with  N,  the  higher 
frequency  wave  components  do  not  exhibit  the  increased  scatter  in  the  data  at  low  N 
and  U)0  values.  Clearly,  die  key  to  understanding  die  relationship  between  die 
ambient  noise  levels  and  the  wave  spectrum  is  in  the  higher  frequencies.  The 
mechanisms  for  sound  generation  are  more  closely  related  to  these  than  to  a  whose 
value  is  determined  more  by  the  lower  frequency  components. 

As  we  discussed  earlier  in  this  chapter,  experiments  on  laboratory  generated  breaking 
waves  [(Rapp  it  Melville,  1990),  (Melville  et  al,  1988),  (Loewen  8c  Melville,  1991a)] 
showed  that  die  steepness  of  die  wave  packet  correlate  with  the  magnitude  of  many 
of  the  variables  associated  with  wave  breaking.  If  we  argue  that  the  laboratory  packet 
steepness  can  be  considered  as  the  analog  of  die  steepness  of  the  wave  field,  then  we 
can  argue  that  the  RMS  wave  slope  s  (or  some  measure  of  wave  field  steepness) 
correlates  Wuh  the  characteristics  of  waw  breaking  in  the  field  as  well.  We  can 
therefore  expect  that,  similar  to  die  laboratory  results,  wave  slope  correlates  with 
ambient  noise  as  well.  When  we  compared  s  with  N  (Section  4.6.2),  we  found  that  log 
s  correlates  with  N  and  that  die  relationship  between  the  two  can  be  described  by  a 
power  law  where 
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As  we  mentioned  earlier,  one  erf  the  key  objectives  of  this  work  is  to  examine  the 
spectral  dissipation  models  and  see  how  they  relate  to  wave  field  parameters  and  N. 
While  the  dissipation  models  based  on  Hawsctsn aim's  (1974)  assumptions  are 
currently  more  widely  used  in  the  development  of  numerical  wave  evolution  models, 
die  dissipation  estimates  based  on  Phillips  (1985)  appear  to  be  grounded  on  a  more 
rational  basis.  However,  although  there  are  dearly  differences  between  the  two,  as  we 
saw  hr>  Chapter  5,  they  both  give  estimates  that  are  comparable  to  one  another.  In 
general,  die  Phillips  (1985)  estimate  is  approximately  twice  that  given  by  Komen  et  al 
(1964).  Oven  that  die  dynamic  range  of  D  is  at  least  4  orders  of  magnitude  over  this 
data  set,  this  is  reasonably  good  agreement 

The  coefficients  of  die  KHH  dissipation  expression  d^(k)  were  obtained  by  tuning 
djiflc)  with  die  wind  input  term  w(lc)  and  die  nonlinear  spectral  flux  term  Vk-T(k)  in 
die  energy  transfer  equation  (Equation  1.11)  to  reproduce  observed  fetch-limited  wave 
growth-  On  die  other  hand,  Phillips'  (1985)  dissipation  expression  dp(lc)  was  obtained 
by  assuming  mi  equilibrium  balance  between  dp(k),  Vk  T(k),  and  w(k).  The  form  of 
die  equilibrium  wave  spectrum  G^co)  -  otu»ga r4  arises  as  a  consequence  of  Phillips' 
(1985)  choices  for  die  analytical  expressions  of  Vk-T(lc),  w(k)  and  d(k).  Comparisons 
between  die  model  O^co)  expression  and  measured  wave  spectra  «(o>).  (Figure  4.21] 
show  that  $£(<s)  gives  magnitude  and  slope  values  similar  to  those  of  $(co). 

Phillips  (1985)  also  argued  that  the  presence  of  wind  drift  c^-u.  inhibits  the 
formation  of  high  frequency  waves  and  therefore  sets  an  upper  limit  to  4>e(<d).  He 
speculated  that  this  upper  frequency  limit  ^  «  rg/u*  -  0.4g/u».  While  die  observed 
wave  height  spectrum  4(co)  does  show  an  upper  limit  to  the  to*4  frequency  slope,  its 
magnitude  is  between  0.1-0.2  g/u»,  approximately  one  third  of  the  value  predicted  by 
Phillips  (1985).  It  is  possible  that  either  die  true  wind  drift  values  are  smaller  than 
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Philttp*  (1985)  estimated  or  that  toone  other  mechanism  inhibits  the  formation  erf 
waves  at  frequencies  lower  than  the  frequencies  wind  drift  would  afreet 

The  major  differences  between  Eh  and  Ep  arise  because  Phillips  (1985)  uses  ^(es) 
instead  of  4(<d)  to  compute  the  total  dissipation.  Although  we  showed  that  the 
differences  between  die  two  spectra  are  generally  small,  these  differences  are 
magnified  because  of  the  cubic  dependence  of  Ep  on  O(to). 

In  addition,  because  Ep  uses  $e(oo),  die  non-constant  variables  which  control  die  value 
of  Ep  are  u»  and  the  phase  speed  of  the  peak  frequency  Cp.  From  Equation  2.17,  we 
note  that 


Ep  -  t&nlftCp/iu)3].  (61) 

The  dependence  on  uu  is  cubic  while  the  influence  of  the  wave  field  through  the 
natural  log  of  the  wave  age  parameter  is  substantially  weaker. 

On  the  other  hand,  the  non-constant  variables  in  the  KHH  estimate  are  the  Pierson- 
Moskowitz  steepness  parameter  a  /**,  the  RMS  wind  wave  amplitude  and 
ttie  mean  frequency  to*, .  Rearranging  Equation  2.27  and  evaluating  a,  it  is  easy  to 
show  that 


(62) 


indicating  that,  using  the  KHH  model,  the  total  dissipation  of  the  wave  held  can  be 
computed  from  integral  wind  wave  field  characteristics. 
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It  is  also  evident  from  Equations  6.1  and  6~2  why,  as  seen  in  Figure  5. Id,  changes  in 
Eh  can  lag  changes  in  Ep.  The  values  of  a*  and  (5*  are  controlled  by  the  enei  &y- 
containing  portion  of  foe  wind  wave  spectrum.  However,  as  we  saw  in  Figures  4.7 
and  48^  «•  more  directly  influences  die  higher  frequencies  which  achieve  full 
development  at  substantially  lower  energy  levels  than  the  spectral  peak. 
Consequently,  these  frequencies  have  a  lesser  influence  on  the  values  of  and  id,,. 

Our  results  show  that,  consistent  with  the  laboratory  measurements  of  breaking 
waves  by  Melville  et  al  (1988)  and  Loewen  k.  Melville  (1991a),  the  ambient  noise  and 
the  different  total  dissipation  estimates  correlate  well  in  the  field  The  relationship 
between  dissipation  and  ambient  noise  was  empirically  shown  to  be 


for  the  Phillips  (1985)  model  and 


P2-# 


for  the  KHH  model 

In  conclusion,  scientists  have  long  recognized  the  importance  of  understanding  the 
physical  processes  in  the  interaction  between  the  atmosphere  and  foe  ocean.  These 
processes  play  an  important  role  in  determining  climate,  the  general  circulation  of  the 
ocean  and  the  atmosphere,  the  formation  of  water  masses,  and  the  transfer  of  gas  and 
moisture  between  the  two  media.  One  of  the  key  issues  in  describing  these  processes 
is  our  ability  to  characterize  the  evolution  of  foe  wave  field  and  how  the  wave  field 


187 


affects  the  airaea  boundary  conditions  and  duxes.  Although  a  lot  of  research  effort 
has  been  made  to  promote  our  understanding  of  these  processes,  one  of  the  major 
obstacles  in  achieving  this  objective  is  our  ability  to  quantify  and  characterize  the 
breaking  wave  dissipation  source  functions.  Our  research  shows  that  the  relationship 
between  sound  and  dissipation  in  breaking  waves  may  prove  to  be  an  important  tool 
in  quantifying  dissipation  in  the  field. 
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APPENDIX 


A.1  Directional  hydrophone  characteristics. 

A. 2.2  The  hydrophone  equation. 

In  this  section,  we  discuss  die  method  used  to  calculate  the  noise  spectrum  level  using 
die  directional  hydrophone  we  used  in  NOBS.  The  ambient  noise  level  N(f)  [dB  re 
pPaVHz]  for  an  equivalent  omnidirectional  hydrophone  can  be  computed  from  the 
raw  power  spectrum  R(f) 


N(f)  *  R(f)  -  SA+  SL(f)+ AG(f )  (3.1) 


where  R(f)  is  in  dB  re  lVoltVlHz,  SA  is  die  signal  amplification,  and  SL  is  the 
sensitivity  level  of  die  hydrophone  in  dB  re  IpPaViVoli2.  SA  was  set  at  50  dB 
diroughout  the  NOBS  observation  period.  The  array  gain  AG  due  to  spatial  filtering 
by  die  directional  hydrophone  is  defined  as  (Dyer,  1989) 


AG  ■  lOlog 


I 


1 


S(Q)dQ 


S(Q)B2(G)dQ 


(3.2) 


where  Q  is  die  solid  angle,  B(Q)  is  the  beam  pattern  of  the  hydrophone  and  S(Q)  is  the 
sound  source  pattern. 
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Figure  A.l.  Plot  of  the  Sensitivity  Level  (SL). 
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Figure  A.2.  a)  Plot  of  the  beam  pattern  of  the  hydrophone  at  1  kHz. 
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Figure  A.2.  c)  Plot  of  the  beam  pattern  of  the  hydrophone  at  4  kHz. 
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Figure  A  .2.  e)  Plot  of  the  beam  pattern  of  the  hydrophone  at  10  kHz. 
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Figure  A.  2.  f)  Plot  of  the  beam  pattern  of  the  hydrophone  at  15  kHz. 
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Figure  JL2.  g)  Plot  of  the  beam  pattern  of  the  hydrophone  at  20  kHz. 
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A.13  CMUbrmtion  mutts. 


bi  this  experiment  we  used  an  ITC-8181A  positioned  at  the  focus  of  a  42"  foam  filled 
parabolic  dish.  The  hydrophone  was  then  calibrated  at  the  Naval  Ocean  Systems 
Center  (NOSQ  in  San  Diego,  California.  The  diagrams  from  the  calibration  report  are 
shown  in  Figures  A.1  and  A.2.  Figure  A.  1  shows  the  plot  of  the  Sensitivity  Level  (SL) 
while  Figure  AJ.  shows  foe  beam  pattern  B^+h)  of  the  hydrophone  at  1, 3, 4, 5, 10, 15 
and  20  kHz  where  is  the  angle  from  foe  axis  of  the  hydrophone.  Radial  symmetry 
with  respect  to  die  hydrophone  axis  was  assumed. 

A.13  Army  gum  computations. 

In  addition  to  tire  beam  pattern  B^Q),  we  also  need  tire  sound  source  pattern  S(ft)  to 
evaluate  AG(f)  (Equation  3.2).  To  estimate  tire  ambient  noise  source  pattern,  we 
followed  the  formalism  of  Urick  (1986,  p.5-lff.)  which  is  reviewed  here.  He  proposed 
that  the  sea  surface  can  be  modeled  as  a  distribution  of  random  and  densely  packed 
sound  sources  each  radiating  sound  with  intensity  (Figure  A.3) 

/(♦Wo  cos>  (A.1) 

where  ♦  is  tire  angle  from  tire  receiver  at  depth  h  and  horizontal  distance  r  to  tire 
source.  For  a  ring  of  sources  of  area  dA,  the  intensity  is 


d  / 


/(♦) 

~A2+r* 


d>4=^~2nrdr. 

I2 


(A.2) 
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horizontal 


.Figure  A. 4.  Urick's  (1986)  cos  9  model  of  the  vertical  distribution  of  ambient  noise 
(dashed  line)  agrees  well  with  measurements  by  Axelrod  et  al  (1965)  at  891  Hz 
(adapted  from  Urick,  1986). 
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From  Figure  A3,  it  is  easy  to  show  that 


r*Atan+, 

f*Asec+. 

dr-Asec^d*. 

Hence 

d/  «  2x/0  cos"  ttan+dt . 


The  intensity  per  unit  solid  angle  Q  is  then 


(A.3) 


(A.4) 


Figure  A.4  shows  a  plot  taken  from  Urick  (1986,  p.  5-10)  comparing  Equation  A.6  with 
the  ambient  noise  directional  pattern  obtained  by  Axelrod,  Schoomer  &  Von  Winkle 
(1965).  The  figure  shows  good  agreement  between  their  data  and  Equation  A.6. 

Integrating  die  denominator  of  Equation  3.2  is  complicated  by  the  different 
orientations  of  the  coordinate  axes  of  S  and  B2.  While  S  is  radially  symmetric  about 
die  z-axis,  die  axis  of  radial  symmetry  of  the  hydrophone  during  NOBS  was  pointed 
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at  an  angle  45°  from  the  z-axis.  Hence  need  to  derive  a  function  that  will  transform 
data  from  the  sound  source  coordinate  system  (r^,0)  to  the  hydrophone  coordinate 
system  (rA^). 

fat  Figure  A.5,  the  position  vector  for  a  point  A  on  the  hydrophone  axis  in  the  sound 
source  coordinate  system  is 


0 

0 

A~  A 

cos45° 

=  A 

* 

0*45“ 

(A.7) 


The  position  vector  for  a  point  Bin  the  sound  source  coordinate  system  is 


B=B\ 


sin  ^  cos  6 
sin^sinO 


cos$ 


(A.8) 


The  coordinate  transformation  for  the  angles  can  be  obtained  from  the  dot  product  of 
AandB 


AB  = 


AB 

H 


(sin^sin0+cos^)  = 


AB 

n 


cos^h 


$h  =cos  *(sin<»sin0+cos4>). 


(A.9) 


Hence,  the  denominator  of  the  equation  for  AG  in  terms  of  the  coordinates  of  the 
hydrophone  axis  is 
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J  S(Q)B2(Q)dQ  »2xJs(4k)B2(4h)d+h 


(A.  10) 


which  can  be  integrated  numerically.  The  numerator  of  AG  can  be  integrated 
analytically  in  the  sound  source  coordinate  system 


J  S(Q)dft  =  2n/0  Jcosfd*  =  2  xl0 . 


(A.11) 


Table  A.1.  Numerical  values  of  AG  and  SL  used  in  Equation  3.1. 


f  (kHz) 

AG 

SL 

2 

4.3 

-162 

3 

A7 

-162 

4 

11.2 

-161 

5 

13.6 

-158 

6 

15.8 

-157 

7 

17.6 

-156 

8 

19.1 

-155 

9 

20.3 

-153 

10 

21.1 

-153 

11 

21.6 

-153 

12 

21.8 

-153 

13 

21.8 

-154 

14 

21.7 

-155 

15 

21.6 

-156 
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Figure  A.6  shows  a  plot  of  AG(f).  The  values  of  AG  computed  from  the  available 
beam  patterns  are  indicated  by  hollow  circles  (O).  We  also  computed  the  value  of  AG 
based  on  an  omnidirectional  sound  source  assumption 

m = 1. 

These  are  shown  as  filled  circles  (•)  in  the  figure.  Figure  A.6  suggests  that,  for  our 
hydrophone,  the  value  of  AG  is  not  sensitive  to  the  differences  the  assumed  source 
level  distribution.  Figure  A.6  also  shows  the  curve  we  used  to  interpolate  the  data  for 
the  frequency  range  we  examined.  The  numerical  values  of  the  curve  are  given  in 
Table  A.l. 
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Figure  A.6.  Plot  showing '  e  array  gain  (AG)  of  the  directional  hydrophone 
confuted  from  the  available  beam  patterns.  Dipole  source  assumed  -  o. 
Omnidirectional  source  assumed  -  •.  The  line  indicates  the  fit  to  the  data  used  in 
computing  N(f)  [Equation  3.1]. 
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